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PAPERS 


FUNCTIONAL DESIGN FLOOD CONTROL 
RESERVOIRS 


The functional design flood control reservoir has its objective the 
determination the relation between the storage space that must provided 
and the corresponding reduction the flood peak. The method incor- 
porating this relation into the study the economic balance flood control 
system has been explained detail Sherman Am. Soc. 
Mr. Woodward’s which heretofore has provided the 
only direct method for the functional design computations, limited its 
application reservoirs such proportions that the maximum outflow 
small compared with the average inflow during the flood, and reservoirs 
with orifice-type outlets and certain type depth-capacity relationship. 
this paper the writers have generalized Mr. Woodward’s method and have 
extended apply reservoirs with either orifice- weir-type outlets, 
valleys wide range morphological configuration. Although the relation- 
ships derived can used the design multiple-purpose reservoirs, the 
present discussion restricted reservoirs for flood control only. 


INTRODUCTION 


necessary preliminary step the estimation the maximum flood run- 
off that may occur the given site. This usually based upon statistical 
study meteorological and hydrological data, and includes determination 
the shape the flood hydrograph, well its volume. Methods making 
these studies have been the subject extended discussion recent 


comments are invited for immediate publication; insure publication the last dis- 
cussion should submitted February 15, 1940. 
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and further treatment not necessary this paper. Suffice say that 
the flood hydrograph used for the design made somewhat greater than the 
maximum probable flood provide margin safety. matter how 
carefully this work done, the design flood will have the kind uncertainty 
that attached any prediction future events. This sets limit accu- 
racy beyond which useless carry subsequent computations. 

order the storage necessary accomplish desired 
tion peak flood flow, necessary route the design flood through the 
reservoir—that is, compute the moderating effect the reservoir storage 
upon the flood hydrograph. 

general, flood-routing methods are two types: 


(1) Step methods, which the flood inflow hydrograph divided into 
steps time intervals such short duration that during each one the inflow 
and outflow rates may assumed constant. The actual depth-capacity 
data for the reservoir site are used and the discharge-depth relation must 
known. irregular inflow hydrograph may approximated, closely 
desired, using large number steps. 

(2) Integration methods, which the given flood hydrograph replaced 
equivalent flood uniform intensity. The depth-capacity data are 
represented empirical formula. necessary know only the ex- 
ponent and approximate zero-depth the discharge-depth relation. 


The first step methods used required trial-and-error solution each 
step. Later, graphical methods were developed which eliminated this cut- 
and-try Results the step computations include the outflow and 
stage the end each step, that step methods are well adapted for opera- 
tion studies. application functional design, however, step methods are 
indirect. necessary first assume the size the outlet. The corre- 
sponding maximum storage and maximum outflow are determined means 
routing computations. Another size outlet then assumed, and the 
process repeated until the desired part the curve representing the relation 
between storage and percentage reduction flood peak has been obtained. 
the other hand, integration methods give the storage required for given 
maximum outflow directly, with considerable saving time and 

Depth-Capacity Exponent.—The relation between depths and capacities 
the reservoir may represented empirical expression the form 


which: the storage capacity the reservoir depth the depth 
water above outlet; constant; and the exponent. 
For the determination logarithmic plotting used. The depths above 
outlet location (zero depth the orifice type outlet, the crest weir 


type outlet) are plotted ordinates and the corresponding capacities 
Rule for Routing Floods through Storage Reservoirs Lakes,” Posey, Engineering 
News-Record, April 25, 1935, Vol. 114, 580. 


the Miami Conservancy Sherman Woodward, Technical Reports 
the Miami Conservancy District, Part VII, Chapter VII. 
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abscissas. straight line best fitting the plotted points drawn. The 
exponent the reciprocal the slope the line, the slope being measured 
directly with scale having uniformly spaced divisions, disregarding the 
logarithmic scale the plot. The constant not needed the analysis. 

similar procedure may followed with the depth-area data, and the 
depth-capacity exponent found integration the empirical relation 
determined. better, however, determine the depth-capacity exponent 
directly, because integration the empirical relation derived logarithmic 
plotting more likely introduce error than numerical summation the 
original survey data. 

studying the depth-capacity curves for more than 250 reservoir sites, 
was found vary from 1.03 following classification, given 
Sutherland, Assoc. Am. Soc. E., was found justified.’ 


Classification Values 


1.0 1.5 
1.5 2.5 
2.5 3.5 
3.5 4.5 


almost all cases, the points lay very nearly along straight line the 
logarithmic plot. The use the exponential formula obtain capacities 
higher levels than those for which survey data are available inadvisable. 
For some sites, the logarithmic plot will follow straight line over con- 
siderable depth, and then deviate sharply one side. For others, the points 
will lie curve. Since volumes near the bottom the reservoir are rela- 
tively insignificant, important that the best straight line should fit the 
data the higher levels most closely. the straight line close approxima- 
tion the data higher levels, discrepancies the lower levels can 
ignored, with inappreciable error the results. Some judgment required 
deciding whether the degree approximation satisfactory. rigid rule 
can formulated. doubtful cases check should made the step 
method. 

Equivalent Uniform Flood.—The equivalent uniform flood corresponding 
given irregular flood hydrograph flood uniform intensity, having 
the same total inflow the given flood, which will produce the same total 
amount storage with the same maximum rate outflow from the reservoir. 
possible determine uniform flood exactly equivalent given irregular 
flood detailed routing computations. has been found, however, that 
even for hydrographs differing considerably from rectangular, not difficult 
select, arbitrarily, equivalent uniform flood that will give results suffi- 
ciently accurate for design purposes. For example, consider the triangular 
flood shown Fig. This flood was routed the step method through 
2.5) with the orifice type outlet. The maximum outflow 
was 5,160 acre-ft per and the maximum storage, 117,100 acre-ft. 
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the duration equivalent uniform flood assumed hr, the 
intensity should equal the total volume the flood divided 40, 
7,320 acre-ft per hr. Routing this uniform flood through the reservoir the 
step method the maximum rate outflow reached was found 5,140 acre-ft 
per hr, which very close the value obtained for the triangular flood. 


12 100 Acre-Ft per Hr 
Design Flood 


= 7320 Acre-Ft per Hr 
| Equivalent Uniform 
Flood 


100 Acre-Ft 


100 Acre-Ft. 
per Hr per Hr 


1.—TRIANGULAR FLOOD AND UNIFORM 


The error the maximum outflow rate and the maximum storage, corre- 
sponding different assumptions for the duration the uniform flood, 
givenin that any estimate for the duration the equivalent 


TABLE THE OUTFLOW RATE 
AND STORAGE 


Maxi 
Duration Percentage error Maximum Percentage error 
assumption, 7’, in acre-feet in In maximum storage, S, in Maximum 
hours per hour per hour outflow rate acre-feet storage 


121,000 +3.3 
7,320 5,140 115,000 


uniform flood between the values and gives satisfactory accuracy, 
this case. Similar results from other comparisons made the writers 
confirm the belief that arbitrary selection equivalent uniform floods not 
difficult. should remembered that there considerable uncertainty 
the shape and volume the hydrograph. 

Outflow Ratio—The outflow ratio, the ratio the maximum outflow 
from the reservoir the uniform inflow the rectangular flood. the 
same the “‘protection defined Mr. The writers recom- 
mend that the name “protection abandoned misleading. The 


higher the value this ratio, the greater the maximum outflow and the less 
the protection. 


THE DETENTION 


The detention ratio the ratio the maximum reservoir storage reached 
during the flood the total volume the uniform flood. definition, then, 


4 | 
40 =" 
12 Hr 36 Hr 
(2) (3) (4) (5) (6) 
4 
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which the detention ratio, the maximum storage, and the 
total volume the uniform flood. The detention ratio function the 
depth-capacity exponent, the outflow ratio, and the exponent governing 
the discharge the outlets. The writers have computed values over 
the practical range the outflow ratio and depth-capacity exponent, for 
reservoirs with orifice-type outlets (Table 2(a)) and for reservoirs with weir- 


TABLE 2.—VALUES THE DETENTION RESERVOIRS 


Values 


(a) 


0.647 0.489 
0.611 0.441 
0.590 0.415 
0.575 0.397 
0.565 0.384 
0.557 0.374 
0.551 0.366 
0.542 0.355 
0.536 


to 


a 
Ww 


(b) 


Noon 


nw 


De 


WWD orOorwe 


100 


1.0 
1.2 
1.6 
2.0 
2.2 
2.4 
2.6 
2.8 
3.0 
3.5 
4.0 
4.5 
5.0 
5.5 
6.0 


type outlets (Table 2(b)). The values were computed with the aid the 
infinite series developed the Appendix. They can checked (indirectly) 
any the standard step methods. 


worth noting that for manually controlled outlets operated 100% 
efficiency, 


ILLUSTRATIVE EXAMPLES 


The significance and utility the detention ratio the functional design 
flood control reservoirs will illustrated two examples. 

Example 1.—Find the relation between the maximum outflow and the 
storage required for controlling flood 50,000 acre-ft orifice- 
outlet reservoir reservoir site Fig. The depth-capacity exponent for 


(1) (2) (3) (4) (6) (7) (8) (9) 
0.399 
0.348 
0.322 
0.303 
0.290 
0.279 
0.271 
0.258 
0.250 
0.955 0.907 0.536 
0.951 0.900 0.509 
0.948 0.893 0.486 
0.945 0.887 0.467 
0.942 0.882 0.451 
0.940 0.878 0.437 
0.938 0.874 0.425 
0.936 0.870 0.414 
0.934 0.867 0.404 
0.933 0.864 0.395 
0.930 0.858 0.377 
0.927 0.853 0.362 
0.925 0.848 0.349 
0.923 0.844 0.338 
0.921 0.841 0.330 
n 4 
e 
) 
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50,000 
The computations are shown Table the column headings which are 


this site 2.25. The average inflow 2,080 acre-ft per hr, 


TABLE 


Maximum 
per hour 50,000 


(2) (3) (4) 


835 0.668 33,400 
1,040 0.582 29,100 
1,250 0.496 24,800 
1,460 0.406 20,300 
1,670 0.312 15,600 


Detention ratio, Storage required, 


from Table 2(a) 


self-explanatory. The curve showing the relation between the maximum out- 
flow and the storage required can plotted from the data Columns (2) 
and (4). 


2.—Assume that, for the same example, political other non- 
engineering reasons dictate that the maximum storage available shall 
25,000 acre-ft. What the maximum flood reduction which may obtained 
the construction reservoir with orifice type Equation 

0.605. Hence, the maximum outflow 0.605 2,080 1,260 
acre-ft per hr. the weir type outlet used instead the orifice type, 
from Table 2.25, and 0.725. Hence, the maxi- 
mum outflow 0.725 2,080 1,510 acre-ft per hr. 


2. 
t (1) 0! 
0.4 
0.5 
0.6 
0.7 
0.8 
150 
15 
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Manually-Controlled has been much argument the 
relative merits manually-controlled and automatic (fixed) outlet structures 
for flood control reservoirs. Comparisons based upon items difficult 
evaluate fallibility should postponed until the relative economy has been 
investigated. Tabulated values the detention ratio afford easy approach 
such study. For example, assume that reservoir site Example 
used reduce the maximum outflow from flood 50,000 acre-ft 
1,000 acre-ft per hr, and that the average cost storage near the top 
the reservoir $30 per acre-ft. interpolation Table the storage 
required, orifice-type outlets are used, found 29,900 acre-ft. With 
manually-controlled outlets, efficiency 100% could presumably 
reached, which case 25,900 acre-ft storage would required. The 
saving storage 4,000 acre-ft, representing $120,000. manually-con- 
trolled outlets are economical, this $120,000 must exceed the capitalized 
cost the maintenance and operation the mechanical gates, plus minus 
any difference the original cost the outlets. 
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APPENDIX 


Evaluation the Detention Ratio.—The inflow into retarding basin, during 
small interval time, equal the sum the outflow and storage during 
that interval time, evaporation and percolation losses being neglected. 
Thus, 


which: average inflow rate during the time dt; average outflow rate 
during the time dt; small increment time; and increment 
storage capacity during the time dt. The outflow from reservoir can 
expressed terms the depth 


which 0.5 for the orifice type outlet, and 1.5 for the weir type. The 
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storage capacity can approximated monomial exponential function 
the depth: 


that 


Substituting these relations into Equation (4), 
idt= 


For constant inflow, during the entire flood period, Equation (7) 
becomes 


Transposing and collecting like terms, 


Equation (9) can transformed into infinite series direct division, 


Integrating, term term, 
B 2 B 3 


the end the uniform flood, when the elapsed time has reached the 
value the storage water the reservoir will have attained its maximum 
value and the outflow rate its maximum value this instant, Equation 
(12) becomes 


[ 
m 
T 
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> 
4 
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the ratio the maximum outflow the uniform inflow denoted the 
“outflow ratio,” the total inflow during the flood period will 


The total outflow during the reservoir-filling period, equal the total 
inflow minus the maximum storage, 


The average outflow during the period 


x a 


m+n m+3n 


The ratio, the average outflow the maximum outflow given 


x 


factoring and direct division the following expression obtained, which 
gives the limiting value approaches zero: 


The detention ratio can evaluated substituting values into the 
equation 


which follows from the definitions and Equations (18) and (20) 

were used computing the values the detention ratio given Table 
Time time empty reservoir after the cessa- 

tion the uniform flood found letting the inflow equal zero Equation 


(17) 
that 
| 
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(8), which then becomes 


Transposing, 


Integrating, 


H 
and, 


1 Cm m—n 


except when zero. From Equation (24), 


The total time required empty the basin after the uniform flood ceases 
found letting equal zero Equation (25) and therefore 


the 
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PAPERS 


GENERAL WEDGE THEORY EARTH 
PRESSURE 


KARL TERZAGHI,! AM. Soc. 


During the twenty years since the first World War the most important 
development earth pressure theory has been the transition from the original 
theory Coulomb broader concept called general wedge theory.” 

According this theory the shape the surface sliding and the position 
the resultant earth pressure are function the degree and type freedom 
for lateral expansion the supported earth. Coulomb’s theory represents 
special case the general wedge theory. Its validity limited the pressure 


earth supports such retaining walls whose lateral resistance inde- 
pendent their amount yield. every other case the general wedge 
theory should used. The only condition required make this theory 
valid almost complete mobilization the shearing resistance the 
supported earth within continuous zone extending from the foot the bank 
upper surface. This single requirement, however, sufficient exclude 
the general wedge theory (or any other theory involving the condition 
incipient failure) from application the pressure soft clay other highly 
elastic types soil the timbering cuts tunnels. 


GENERAL PRINCIPLES 


The essentials the general wedge theory can described follows: 
the lateral support mass earth yields, wedge-shaped body earth 
descends toward the yielding support. Before this process begins, however, 
the weight, the wedge (see Fig. 1), the resultant reaction, along the 
surface sliding, and the reaction acting along the back the lateral sup- 
port, must equilibrium, and the shearing resistance the backfill must 
almost completely mobilized within continuous zone, extending from the 


paper was presented the meeting the Soil Mechanics and Foundations Division 
Rochester, Y., October 14, 1938. comments are invited for immediate publication; in- 
sure publication, the last discussion should submitted February 15, 1939. 
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foot the backfill its upper surface. This broader concept the 
wedge theory involves assumption concerning the shape the surface 
sliding and none concerning the location the center the resultant earth 
pressure. sliding must satisfy the condition that the correspond- 
ing reaction, required prevent the descent the wedge, shall 
mum and that the yield the lateral support shall consistent with the con- 
ditions stress and strain within the wedge. analysis the manner 
which the lateral support yields, and the influence this yield the conditions 
stress and strain the wedge, was presented the writer This 
analysis led the conclusion that the shape the surface sliding and the 
position the resultant earth pressure are function the degree and type 
freedom for lateral expansion the supported earth. Depending this 
degree freedom, the surface sliding can assume any shape intermediate 
between inclined plane (Fig. and conchoidal surface approaching 
vertical asymptote (Fig. 1(d)). The center the resultant earth pressure 
may anywhere between 0.33 and and the corresponding distribu- 


qe 


tion the lateral pressure ranges between hydrostatic and roughly 

the general wedge theory, the supplementary assumptions are intro- 
duced that the surface sliding plane and that the distribution the pres- 
sures over both the back the wall and the surface sliding hydrostatic, 
the result the well-known earth pressure theory Coulomb. Analogous 
the introduction these assumptions the general wedge theory the 
assumption that plane sections remain plane, which Navier introduced the 
general theory bending obtain the simplified theory bending. this 
simplified theory applied the bending short and very thick beams, 
erroneous results are obtained. Similarly, the theory Coulomb leads un- 
satisfactory results the conditions the problem are incompatible with the 
fundamental assumptions. 


Fundamental Fallacy Earth Pressure Computations,” Karl Terzaghi, Journal, Boston 
E., April, 1936. 
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order define the boundaries between the realm validity Coulomb’s 
original theory and that the general wedge theory, the writer presents herein 
analysis the different conditions that are likely influence the intensity 
and the distribution the active earth pressure. For the sake brevity the 
analysis limited the investigation pressures exerted cohesionless back- 
fills vertical lateral supports. Neither the influence surcharges nor that 
percolating water the lateral pressure will taken into consideration. 


INNER LOWER EDGE 


the back wall perfectly frictionless and if, addition, the wall 
yields tilting around its lower edge, Coulomb’s assumptions are fully satisfied 
and the results obtained his method are correct. This case illustrated 
Fig. 1(a). the back the wall rough, involving angle wall friction 
the forces that act the wedge according Coulomb’s theory are shown 
Fig. 1(b), and the corresponding distribution the lateral pressure over the 
back the wall hydrostatic. The forces shown Fig. 1(b) fail intersect 
one point, which incompatible with the fundamental conditions for the 
equilibrium the wedge. This contradiction inherent Coulomb’s theory 
due the fact that the lower part the real surface sliding curved slightly 
shown Fig. The corresponding distribution the lateral pressure 
also hydrostatic. accurate computation the true shape the surface 
sliding and the true intensity the lateral pressure for the case illustrated 
Fig. 1(c) was accomplished Theodor von Am. Soc. E., 
used different methods, but the results are essentially identical. They demon- 
strate that the error involved Coulomb’s theory does not exceed several 
per cent. For practical purposes such error without significance. 

Similar investigations have been conducted concerning passive earth pres- 
sure (resistance sand against lateral displacement). When applied wall 
with smooth and vertical back (Fig. 2(a)), Coulomb’s theory passive earth 
pressure strictly correct. the other hand, the theory applied 
similar wall with rough back, shown Fig. 2(6), Coulomb’s assumptions 
concerning the forces that act the wedge are conspicuously contradiction 
the conditions for the equilibrium the wedge. Both experiment and the 
advanced theories revealed von and von demonstrate that 
the surface sliding curved shown Fig. 2(c). They also demonstrate 
that the error due Coulomb’s assumptions may exceed 30%, provided the 
value the angle wall friction approaches that the angle internal 
friction, Hence, exceeds several degrees, necessary modify the 
theory assuming that the sliding surface not plane. For all practical 
purposes, one can approximate the shape the sliding surface considering 
the lower portion segment logarithmic spiral and the upper portion 
straight line rising angle (45 0.5 the horizontal, shown 


elastische Theodor von Proceedings, Second International 
Cong. for Applied Zurich, 1926. 


klassische Erddrucktheorie. Veziigyi von Jaky, Budapest, 1931. 
Theorie des Ohde, Die Bautechnik, 1938. 
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Fig. The position the center the spiral must determined 
such manner that the corresponding value the passive earth pressure 
minimum. The most convenient method for accomplishing this result 
trial. accurate method determining the shape the surface sliding 
was introduced 1938 However, the application this method 
requires considerable time and labor, and the difference between the results thus 
obtained and those obtained means the aforementioned approximate 
method does not exceed few per cent. The assumption hydrostatic 
pressure distribution reasonably correct, provided the wall advances toward 
the fill tilting around the lower edge. 

Far more important than the errors due the approximation regarding the 
shape the surface sliding are those due the uncertainties concerning the 
angle wall friction. Thus far the results pertinent investigations seem 
indicate that the angle between the resultant passive earth pressure and the 
normal the wall considerably smaller than the angle friction between 


backfill and wall unless the state stress the mass sand subject the 
later pressure approaches the point failure. Due the important bearing 
the subject passive earth pressure the stability bulkheads, further 
experimental investigations concerning the angle wall friction are urgently 
needed. 

the cases illustrated Fig. 1(a) 1(c) the distribution the active 
earth pressure over the back the wall hydrostatic. The theories which 
apply these cases may also used for computing the active earth pressure 
retaining walls other lateral supports which are capable yielding inde- 
finitely without being weakened. Large-scale earth pressure tests made the 
writer 1929 demonstrated conclusively that the distribution the lateral 
pressure sand always hydrostatic after the lateral support has yielded 
beyond certain critical distance which depends the height the fill and 
the character the 


Retaining Wall Tests. Pressure Dry Karl Terzaghi, Engineering News- 
Record, February 1934. 
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the other hand, the capacity the timbering cut yield under 
pressure very limited and case can expected reach the value 
the yield required induce hydrostatic pressure distribution. For this 
reason necessary consider the influence small yield lateral support 
under the conditions associated with excavation and timbering the intensity 
and the distribution the lateral pressure. 


TIMBERING EXCAVATIONS SAND 


early 1908, the late J.C. Meem, Am. Soc. E., the basis his 
experiences subway construction New York City, maintained that the 
distribution the earth pressure over the timbering cuts decidedly non- 
hydrostatic, and that the surface sliding intersects the surface the bank 
right His statements were corroborated similar practical ex- 
making his large-scale earth pressure tests 1929 the writer observed that the 
surface sliding the backfill sliding wall has abnormally small radius 
curvature and that the center the earth pressure acts elevation 
0.45 above the base backfill with height 1935 the writer demon- 
strated theory,® the assumption plane surface sliding, that the 
center the active earth pressure may any elevation between 0.33 and 
0.66 above the base the backfill, the exact position this point depending 
the degree freedom the supported earth expand horizontal 
also arrived the conclusion that the type and degree 
freedom for lateral expansion mass sand during the process excavation 
and timbering practically identical with that backfill behind wall 
which yields slight tilting around its upper edge. 1938 Ohde published 
method for computing the lateral pressure which corresponds that type 
ayield. for this case that the center pressure should height 
0.55 above the base the fill and that the lateral pressure appreciably 
greater than the Coulomb pressure. The last mentioned fact can expressed 
equation 


which horizontal component the total lateral pressure; unit 
weight the backfill; depth the cut; hydrostatic pressure ratio 
for Coulomb conditions; and the wedge factor, factor that expresses the 
influence the state strain the wedge the total intensity the lateral 
pressure. 

The value can computed means Coulomb’s theory. 
depends only the angle internal friction and the angle wall friction, 
The wedge factor depends the state strain imposed the wedge 


Bracing Trenches and Tunhels, with Practical Formulas for Earth Pressures,” the late 
J.C. Meem, Transactions, Am. Soc. E., Vol. LX, June, 1908, 


Eng and Rock Pressure,’’ Moulton, Transactions, Am. Inst. Mining and Metallurgical 


the Lateral Pressure Sand the Timbering Cuts,” Karl Terzaghi, Pro- 
ceedings, International Conference Soil Mechanics, Vol. Cambridge, Mass., 1936. 
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the lateral boundary conditions. can computed means 
theory, which based the assumption that the surface sliding can 
replaced However, the method far too complicated for 
universal practical application. For this reason the writer developed sim- 
plified method computation. based the empirical facts that the 
surface sliding intersects the horizontal surface the bank right angle 
and that the center the lateral pressure the timbering cut with depth 
bottom the cut. The agreement between Coulomb’s and von 
theory the active earth pressure and between Ohde’s rigorous and simplified 
theory the passive earth pressure prove that moderate departure the 


(W) 


assumed shape the surface sliding from the real one has little influence 
the results. the case under consideration the writer found that the most 
convenient substitute for the real surface sliding logarithmic spiral with 
the polar equation 


which length vector through the center the spiral; length 
arbitrarily selected reference vector; base natural logarithms; 
angle between and ro, radians; and, angle internal friction 
the backfill. Fig. 3(a) shows such spiral corresponding Equation (2). 
has the property that vector passing through the center the spiral and 
any arbitrary point the spiral intersects the tangent the spiral 
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and the normal the tangent point equal the angle internal 
friction. 

order satisfy the boundary conditions the problem, the spiral 
oriented with reference the face the bank Fig. such manner 
that the spiral passes through the foot the bank Fig. and that 
intersects the horizontal surface the bank right angle, arbitrary 
distance (d) from the upper rim The angle wall friction, and the eleva- 
tion the center pressure are assumed known. The general wedge 
theory requires that the frictional resistance the surface sliding should 
almost fully mobilized over the entire surface sliding, Fig. 3(b). Asa 
consequence the reaction, each element the sliding surface Fig. 
intersects the tangent the surface sliding angle According 
Equation (2) every vector also intersects the corresponding tangent the 
same angle, All the vectors pass through the center the spiral 
(point Fig. Hence all the reactions (dF) and the resultant reaction, 
(F), pass through point Force due the weight the wedge 
passes through the center gravity area Taking moments around 
the center the spiral, the following formula obtained condition for the 
equilibrium wedge bc: (W) (dw); 
or— 


(dw) 

The remainder the investigation practically identical with the graphical 
solutions Coulomb’s theorem. Several spirals are assumed intersect the 
horizontal surface the bank different distances (d) from the rim, the 
order facilitate tracing the spirals the writer plots Equation (2) 
piece cardboard and cuts out used pattern for tracing the 
spirals the graph. For each one these spirals determines the corre- 
sponding value and plots these values ordinates above (c) Fig. 
Thus obtains curve According the well-known procedure for graphic 
earth pressure computations, the active earth pressure equal the ordinate 
the highest point curve and the position the corresponding surface 
sliding determined point beneath point 

One would expect the greatest error involved the graphical method il- 
lustrated Fig. occur when the method applied the computation 
the active earth pressure wall, Fig. which yields tilting around its 
lower edge, because this extreme case the surface sliding does not resemble 
logarithmic spiral. The corresponding numerical values are and the 
wedge factor introducing into his computation the values tan 
0.6 and the writer obtained 0.94 instead the correct value 
1.0. Curve Fig. 3(c) represents the curved sliding surface for this case and 
line the corresponding surface sliding obtained means Coulemb’s 
theory. spite the conspicuous difference between these two lines, the 
error does not exceed 6%. the value high, the surface sliding 
actually has the shape assumed the general wedge theory. Therefore the 
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error can expected still smaller. For 37°; 
the computation furnished 1.27 and for 0.50, 
the value increases with increasing values and 


1.18. 


FIELD MEASUREMENTS 


order check the writer’s theory applies the pressure 
exerted the timbering excavations, the Siemens Bau Union, 1936, 
measured the pressures the struts 38-ft subway cut Berlin, 
The sides the cut consisted mostly sand and fine gravel. Fig. 4(a) shows 


(a) CROSS SECTION LONGITUDINAL SECTION 


cross section through the cut, and Fig. 4(b) shows longitudinal section. 
Before excavation was begun, two outer and three inner rows I-beams were 
driven into the ground. the outer rows the spacing was and the 
inner rows was about ft. excavation proceeded, horizontal boards 
were inserted between the flanges the steel beams the outer rows. The 
struts were arranged brace the steel beams the outer rows against 
one another across the cut. After excavation was finished the lateral pressure 


the earth was taken four sets struts. every set the average hori- 
zontal distance between the struts was 6.5 ft. The vertical distance between 
the sets decreased from the top set the bottom set. Each strut consisted 
four pieces bearing against the web horizontal channels and fastened the 
steel I-beams the inner rows. The distance between the webs each pair 
channels was maintained means wooden spacers (a, Fig. 5). 
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order measure the axial pressure strut, two small hydraulic jacks 
were inserted between the channels (Fig. 5(b)) both sides the spacer. 
Then the jack pressure was increased gradually and the corresponding increase, 
Al, the distance between points and was measured. The relation 
between the total jack pressure, and the corresponding values are 
shown Fig. long the pressure exerted the jacks less than the 
pressure the strut, the jacks merely relieve the pressure the spacer without 
changing the pressure the long sections the The corresponding 
increase the distance between points and slight, shown the steep 
rise section the (P; Al)-line Fig. However, soon the jack 


Fia. 


pressure exceeds the pressure the struts, the pressure the long sections 
the strut increases and, consequently, the deflectometer readings increase 
rapidly (see flat section the (P; Al)-line Fig. 6). The pressure that 
acted the strut before the jack pressure was applied determined the 
ordinate the point intersection the vertical axis the graph Fig. 
and the prolongation section the (P; Al)-line. Fig. view the 
test which represents the struts, the spacer, 
the jacks, the extensometer, and the channels. 

The measurements were made the struts six places and 


die Messung der einer Baugrubenaussteifung,” Spilker, Die Bau- 
1937, Heft 
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were limited the top strut and profile the second strut from the 
Fig. shows the results some the measurements. 

order discover whether the frictional resistance the supported banks 
had been mobilized, comparison was made the writer between the theoreti- 
cal value and the measured value the total pressure the struts the dif- 


ferent sections. For this purpose, Equation (1) was used. The wnit weight 
the was 1.7 tons per the angle internal friction was 37°, and the 
angle wall friction was values the angles and undoubtedly 
were estimated and not measured. Due the absence fine material, the 
cohesion was negligible. 

With the numerical values mentioned, Coulomb’s theory yields the value 
0.204; and, from the general wedge theory, 1.22. The computa- 
tions, according the general wedge theory, yield theoretical value for the 
total lateral pressure which 
average 25% greater than the 
value obtained from the measure- 
ments. These findings seem in- 
dicate that the angle internal 
friction the supported material 
was appreciably greater than the 
estimated value 37° (at least 
41°), and that the shearing resis- 
tance the supported banks was 
almost fully mobilized, satisfying 
the only fundamental requirement 
the general wedge theory. 

According the general wedge 
theory for the pressure sand 
timbering excavations, the dis- 
tribution the pressure over the 
lateral support should more 
less parabolic, having maximum 
pressure approximately half the 
depth the cut. This conclusion fully substantiated the results the 
measurements shown Fig. 

The data shown Fig. also indicate the magnitude the deviations the 
true pressure distribution from the statistical and theoretical average. 
ing the geological profile the cut, the supported banks are fairly uniform. 
Hence, seems most likely that the observed variations are essentially due 
accidental variations the procedure followed constructing the timber work. 
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Foremost among these variations undoubtedly the variation the force 
with which the wedges were driven. Another cause may the local variations 
the density the supported material, because such variations are likely 
induce local arching. This might account for the unusual shape pressure 
4in Fig. the other hand, the location the center the lateral 
earth pressure was quite remarkably consistent. varied between the rela- 
tively narrow limits and above the bottom the cut. The 
most frequent values were between 0.53 and 0.55 

These variations must considered inevitable. For this reason the writer 
has always insisted that ample margin provided for possible scattering. 
Fig. this margin indicated the trapezoid indicates the 
depth the cut and the horizontal component the total lateral pressure 


an 


the timbering, the width the trapezoid equal 1.4 Fig. shows 


that the observed variations are least large those which are taken care 
the suggested margin. 


FLEXIBLE BULKHEADS 


Another problem outstanding practical importance concerns the earth 
pressure exerted flexible bulkheads. The total earth pressure can ex- 
pected equal that which acts the back yielding, non-flexible 
retaining wall; but the distribution the pressure quite different, because the 
lateral bulging induces concentration the pressures the vicinity the 
upper and the lower support the expense the pressure the middle, 
yielding part (see Fig. 1(e)). This fact seems have been noticed for the first 
time Danish engineers. When examining some their century-old bulk- 
heads they found out computation that these bulkheads could not possibly 
stand the lateral pressure the supported earth unless the distribution the 
pressure similar that shown Fig. This observation led radical 
revision the Danish specifications for the design case 
sheet-pile computations which came the attention the writer, the value 
for the maximum bending moment obtained means the Danish method 
was about one half the value obtained the basis Coulomb’s theory. The 
economy inherent the application this procedure obvious. 

1936 Stroyer published the results experiments with flexible 
steel bulkheads supporting dry sand backfills The results these 
tests can also explained only the assumption that the distribution the 
lateral pressure over the bulkheads was such that shown Fig. 1(e), which 
turn requires intense within the sand behind the bulkhead. Some 
engineers still seem believe that the arching effect temporary phenomenon, 
destined disappear time. The Danish experience mentioned previously 


des murs Proceedings, World Eng. Cong., Rept. No. 416, 
1929. 

Regler for Beregning Jaernbetonkonstruktioner published the 
Dansk (no date); see also ‘‘Ztschr. des Verb. intern. Schiffahrtskongr.,” No. 
January, 1929. 
Pressure Flexible Stroyer, Paper No. 5024, Journal, Inst. E., 
November, 1935 (discussion this paper the writer, January 1936). 
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represents one many empirical facts which invalidate this pessimistic opinion, 
Some the Danish bulkheads were exposed for centuries the pounding 
waves, and for decades heavy traffic, without the distribution the lateral 
pressure becoming 

The arching between the upper and the lower support flexible bulkheads 
eliminates Coulomb’s theory from application bulkhead problems. The 
influence the lateral deflection the location the center pressure 
not yet known. Undoubtedly, depends not only the type backfill 
and the flexibility the wall, but also the distance through which the 
upper and the lower support the bulkhead are capable yielding and 
the density the backfill. 

Because the outstanding practical importance the problem, the writer 
advised the Soil Mechanics and Foundations Division the Society include 
its research program. The variables considered such experi- 
mental investigation are: (a) Degree fixation the lower part the bulk- 
head, (b) position the upper support (anchorage) with reference the top 
the backfill, (c) yield the anchorage, (d) deflection the flexible part 
the bulkhead with reference plane through the upper and lower supports, 
and (e) the relative density the backfill. Investigations this kind could 
made well-equipped earth-pressure laboratory small expense for 
additional special equipment. would also advisable include the 
program the measurement the deflection full-sized bulkheads, compute 
the corresponding bending moments supplement the results the 
measurements determining the physical properties the backfill. 


PRESSURE EXERTED COHESIVE EARTH 


For cohesive soils the general wedge theory based Coulomb’s equation 
for the shearing resistance, cohesive soils. this empirical equation, 


Coulomb’s cohesion value, the unit pressure normal the surface 
sliding, and angle internal friction. However, there are many 
conditions under which Equation (4) does not give even approximate values. 
Hence, the very best, should considered crude approximation 
reality. used without careful discrimination, can seriously 
For cohesive soils, the influence the type yield the lateral support 
the location the center pressure, and the distribution the lateral 
pressure, still more important than for cohesionless ones. the lateral 
support yields tilting around its lower edge, one must anticipate hydro- 
static pressure distribution, beginning some depth below the upper edge- 
Above this depth the soil can stand without any lateral support. the other 
hand, the lateral expansion the wedge increases with depth (as does dur- 
ing the excavation timbered cut) the heavy pressure acts the upper part 
the timbering, whereas the lowest part the sides can left without any 


des Porenwasserdruckes auf den Scherwiderstand der Tone,” Karl Terzaghi, Deutsche 
Wasserwirtschaft, 1938, Heft 8/9. 
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support. These facts can demonstrated easily means the theory 
analogy with reference the theory cohesionless soils, one 
may also conclude, accordance with experience, that the section through the 
surface sliding should start with vertical tangent the surface the bank 
and descend toward the foot the bank with very marked curvature. Hence, 
there contradiction between the general wedge theory and experience, 
provided the yield the lateral support sufficient mobilize the major 
part the shearing resistance the supported banks. 

This mobilization the shearing resistance, however, can only antici- 
pated comparatively dense strata sandy soils. the soil fairly elastic, 
the case for plastic soft clay, the lateral expansion the sides the cut 
during excavation and timbering can scarcely expected result mobilizing 
more than fraction the shearing resistance the supported material. 
this case the intensity the lateral pressure the earth must greater than 
the value obtained means the general wedge theory. The importance 
the difference depends the elastic properties the supported soil, the dis- 
tance through which the material can yield during the process excavation 
and timbering, the state stress that existed the earth before excavation was 
begun, and the time factor. Because the difficulties connected with 
evaluating these many variables, and the complicated elastic properties 
cohesive soils, the prospects predicting the lateral pressure successfully 
theory are slight. substantial advancement pertinent knowledge can 
expected only from actual measurement the pressure the braces cuts 
through different clay soils, and subsequent correlation the results with the 
physical properties the materials encountered the cuts which the mea- 
surements were made. satisfactory method for making such measurements 
described herein, and illustrated Figs. steel beams are used 
support the sides cut, the pressure the struts can determined means 
strain gage such the Whittemore strain gage. The reference points 
should placed and their distance recorded before the beams are lowered into 
the cut. 

account the appalling scarcity published information concerning 
the pressure clay timbering excavations, the results every reliable 
observation regarding this subject should considered important. Therefore 
the writer has requested the Soil Mechanics and Foundations Division collect 
such data and promote the measurement the pressure the timbering 
excavations wherever there opportunity so. The results thus 
obtained could not fail produce radical revision the specifications for the 
computation earth pressure under such conditions. 
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AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


PAPERS: 


SEWAGE DISPOSAL PROJECT BUFFALO, 
NEW YORK 


The purpose this paper make record, with proper explanation, 
some the engineering data accumulated during the design and construction 
the sewage disposal project Buffalo, During the three years 1936 
1939, the Buffalo Sewer Authority completed the construction intercepting 
sewers, sewage treatment works, and storm-water relief sewers, total proj- 
ect cost $15,000,000. The sewage disposal works were placed operation 

toward the end June, 1938. 


GENERAL STATEMENT 


The Buffalo Sewer Authority was created special act the Senate and 
Assembly the State New York and became law April 1935. The five 
members the Authority were first appointed June, 1935, and later the 
year, December the Authority appointed the organization which was 
serve designing and supervising engineers. Meanwhile, June 18, the 
Authority made application the Federal Emergency Administration 
Public Works for loan and grant, the amount $15,000,000, aid 


and connecting, lateral and outlet sewers, storm water drains, pumping 
and ventilating stations, disposal treatment plants works, and other 
appliances and structures provide effectual and advantageous means 
for relieving Niagara River, Buffalo River and Lake Erie from pollution 
the sewage and waste the City Buffalo and relieving the City 
Buffalo from inadequate sanitary and storm water drainage and for the 
sanitary disposal treatment the sewage the City 


The estimates cost this application included construction work the 
sum $13,200,000, and sum $1,800,000 cover preliminary expenses, 
land and right way, engineering, interest during construction, legal, and 
administrative costs. The proposed storm-water relief sewers were for the 
Bird Avenue Sewer District only. 


comments are invited for immediate publication; insure publication the last 
discussion should submitted February 15, 1940. 
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The first investigation the sewage disposal problem Buffalo was made 
between 1909 and 1918 the International Joint Commission connection 
with investigation the pollution boundary waters. 


The conclusions 
stated their final report are briefly summarized follows: 


“Tt advisable consider the Niagara and Detroit Rivers 
themselves, far remedial and preventive systems are concerned. 
has been stated, the most serious condition existing the bacterial 
pollution these streams. remedy this evil, sewage treatment should 
applied connection with dilution far necessary bring their 
waters the standard mentioned—a mean annual cross-sectional average 
coli not exceeding 500 per 100c.c. This necessary sewage purification 
can effected fine screening sedimentation, and when necessary 
chemical disinfection, cost which will impose unreasonable burden 
upon the urban community responsible for the present pollution.” 


Furthermore, has been clear for number years that the waterways 
within the city limits Buffalo have been grossly polluted, including Buffalo 
River, Buffalo Harbor, and Black Rock Harbor. The function the sewage 
disposal part the project, therefore, threefold, follows: 


(a) relieve Buffalo River and its immediate tributaries from undue 
pollution; 

(b) relieve Buffalo and Black Rock Harbors from undue pollution; and, 

(c) remove contamination from Niagara River and reduce the bac- 
terial content its waters order avoid placing undue pollutional load 


water purification works communities below Buffalo which derive their 
water supply from the Niagara River. 


number the existing sewers had also been overloaded and much need 
relief. this was especially true the existing sewers the Bird Avenue 
District, was decided include storm-water relief sewers for this sewer 
district the project. 

Action for undertaking the project was stimulated the orders the 
New York State Department Health, issued March, 1934, which required 
that during 1935 the City Buffalo must engage engineers, take steps 
acquire land which construct sewage treatment works, submit plans for 
the Bird Avenue storm-relief drains, and make satisfactorily comprehensive 
general plan for interception and treatment the entire sanitary sewage 
the city. Detailed plans and specifications for substantial part the work 
were required completed April 1936, and construction work started 
before July 1936. The order required that the treatment works 
completed and placed operation before October 1937. 

Actual construction work was begun part the Bird Avenue storm- 
water relief drains December, 1935, and contract for borings, which ulti- 
mately required the expenditure about $38,000, was started January, 
1936. The first major construction contract went into actual construction 
June, 1936, under bids which had been advertised April 23, 1936. 

average monthly rate construction was maintained during twenty-two 
consecutive months (from August, 1936, through May, 1938) $547,000. 
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Buffalo situated along the eastern shores Lake Erie and the Niagara 
River (see Fig. southern part the city cut the Buffalo River 
which extends from the harbor southeasterly and easterly direction. Much 
the area tributary Buffalo River occupied industries, and this the 
principal manufacturing district the city. Scajaquada Creek flows east 
and west direction through the northerly part the city and used form 
Park Lake Delaware Park. The northerly part the city much more 
residential character. The population Buffalo the United States 
Census Bureau, and its estimated probable growth, follows: 


Estimated 
Year Population Year population 
352,387 
423,715 
506,775 
573,076 


The total area 1930 was 26,980 acres and the density population was 21.3 
persons per acre. was thought likely that the area the city might 
extended somewhat the north and east the distant future with the larger 
populations, but this does not influence the plans. The length the city ina 
north and south direction 8.6 miles and the width east and west direction 
about miles. 

The average temperature and precipitation Buffalo are indicated 
Table which summary data from the Weather Bureau. The 


Jan- Feb- 
uary ruary 


Sep- No- De- 

Au- 

tem- vem- cem- Total 
gust ber tober ber ber 


Description March| May June July 


Average temperature, degrees Fahrenheit. Average total precipitation, inches. 


Weather Bureau also reports that the prevailing winds are southwesterly 
except during November, December, and January, when westerly winds 
prevail. 

SEWERS 


Practically all the sewers Buffalo are the combined plan, their total 
length being about 710 miles. The first sewers were built prior 1850; 
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during the 5-yr period prior 1855, miles were constructed; and the 
period prior 1895, 139 miles were built. Thus, the existing sewers 
Buffalo have been service, general, for many years. present (1939) 
there are about outlets, all which have now been connected the new 
intercepting sewers. 

The city divided topographically into two parts, one draining naturally 
the Buffalo River and the other directly Black Rock Harbor and the Niagara 
River. With few exceptions, the main drainage system has developed along 
these natural lines. the area north the lower part the Buffalo River, 
system canals formerly connected with the old Erie Canal. Some the 
city sewers discharged into these canals, which have been covered over and 
serve storm overflow sewers. Except the case District VII (see Fig. 1), 
the sewers Buffalo are the so-called combined system, with sewage and 
storm water carried the same sewers. District VII the dry-weather 
flow confined part separate system sewers and conducted the 
Hamburg Street Pumping Station. This station contains steam power 
plant with three horizontal, centrifugal pumps, driven steam engines. 
was formerly operated conjunction with refuse incinerator, which has 
not been used regular operation since before 1929. 

Conditions the harbor, the canals, and the Buffalo River became objec- 
tionable early date, and alleviate this pollution, the so-called Swan 
Trunk Sewer was built 1883. This sewer, with later extension, runs 
roughly parallel the Buffalo River and the harbor from Bailey Avenue 
Albany Street, where passes under Black Rock Harbor the Niagara River. 
diameter for most its length, and was designed intercepting 
sewer receive the dry-weather flow from all the area north Buffalo River 
and south the Bird Avenue sewer district. also receives the discharge 
the Hamburg Street Pumping Station. later years has been much over- 
loaded and, especially near the upper end, overflows have been frequent 
sometimes continuous during dry weather. The capacity this sewer 
utilized the new intercepting sewer system. 

The population and industries adjacent Buffalo River have continued 
discharge their sewage into the sluggish waters the river. Nearly all the 
large industries along Buffalo River have private sewers directly the river. 

There are several small areas, chiefly small separate municipalities, 
outside the present city limits which discharge sewage into the Buffalo sewers, 
amounting the total about 1,000 acres, with population about 13,000. 

The characteristics the major sewer districts Buffalo are shown 
Table 


The waterways Buffalo pertinent the sewage disposal project are (see 
Fig. 1): Niagara River, Buffalo Harbor, Black Rock Harbor, Buffalo River, 
Cazenovia Creek, and Scajaquada Creek. 

Niagara drainage area the Niagara River, the outlet 
Lake Erie, about 175,000 miles. Based records the Lake 
Survey, the average discharge 203,000 discharge 176,000 
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per sec exceeded 90% the time. The minimum monthly average 
discharge record (142,000 per sec) occurred February, 1934; and the 
maximum discharge (251,000 per sec) occurred June, 1876. Thus the 
discharge relatively large and uniform. 

The elevation the water surface the Niagara River opposite the sewage 
treatment plant Bird Island covers quite range. estimated that 
elevation 9.7 (Buffalo City Datum) exceeded the average for 


TABLE 2.—CHARACTERISTICS AND EsTIMATED AVERAGE SEWER 
Districts 


AVERAGE 
District 
No. (see District 
Fig. 1930 1985 


1930 1985 


Ontario-Hertel 


32.9 50.0 860 860 103 108 
31.9 46.0 730 730 160 196 
3.0 40.5 800 800 285 252 
VI Michigan-Emslie............... 55.3 68.7 1,060 1,060 | 10.3} 186 13.8 | 201 
Hamburg Street Pumping 20.3 14.0 760 384 578 
IX 40.2} 150.8 2,145 | 4,830] 10.5] 261 28.6 | 190 
x 43.2} 115.8 2,410; 4,140] 13.3) 308 27.3 | 236 


Other non-sewered areas 


million gallons daily. Ged gallons per capita per day. 


the time and elevation 5.7 exceeded for only the time. 
Extraordinary fluctuations the level Lake Erie cause momentary peaks 
the water level the river Bird Island much higher than the foregoing. 
Thus elevation 2.0 reached, the average, once three years, 
and elevation 5.0 has estimated frequency occurrence once 
fifty years. 

Below Buffalo, the Niagara River used water supply Tonawanda, 
North Tonawanda, Lockport, Niagara Falls, Youngstown, and Fort Niagara, 
New York, with total population about 150,000, and the cities Fort 
Erie and Niagara Canada. There increasing use the water for 
recreation, including boating, fishing, summer homes, and some bathing, mainly 
limited, present, the Canadian shore and the westerly shore Grand 
Island. The average number coli per cubic centimeter the Niagara 
Falls water-works intake was 36.7 1936 and 27.5 1937, with individual 
samples high 240. The analyses the Sewer Authority laboratory show 
great concentration 2,400 coli per the most contaminated 
sections the river downstream from Buffalo. The lowest content dis- 
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solved oxygen found during the river surveys the Buffalo Sewer Authority 
July, 1936, was 6.0 ppm. 

Buffalo Harbor.—Buffalo Harbor made artificially breakwaters Lake 
Erie. Little sewage discharged directly into Buffalo Harbor and has not 
been important factor the design the project. 

Black Rock Harbor.—On the other hand, Black Rock Harbor has received 
considerable quantity sewage and was seriously polluted. fact, 
waterway connecting the New York State Barge Canal with Buffalo Harbor 
and relatively long and narrow, having the shape more canal than 
harbor. Its length 16,500 ft, the area average cross section 8,800 
ft, its surface area 8,100,000 ft, and has liquid volume mean 
depth 140,000,000 ft. Some circulation (about 100 per sec) 
maintained the operation small outlet the lower northerly end 
the harbor and operation the lock. Variations the level Lake Erie 
also have marked effect mixing relatively unpolluted lake water with the 
harbor water, especially the upper southerly end. 

The principal source pollution Black Rock Harbor was the Bird 
Avenue outlet which discharged about 25.0 mgd sewage into the harbor about 
3,000 south its This about 20% the sewage the entire 
city. indication the effect this pollution given the data 
Table 


TABLE Rock WATER, FROM INDIVIDUAL 
SAMPLES 


Date Dissolved oxygen, 5-day B.O.D., Dissolved oxygen, 5-day B.O.D., 
percentage parts per percentage parts per 
saturation million saturation million 
1936 
September 30........ 2.8 
December 3......... 8.1 
1937 
June 2.. 2.3 
September 8......... 2.6 
December 7......... 4.6 


Buffalo River and Cazenovia Creek.—Buffalo River and its main tributary, 
Cazenovia Creek, have drainage area miles. There continuous 
record gagings the quantity water flowing Buffalo River. However, 
the Geological Survey has record routine gagings Cattaraugus 
Creek, the drainage area which adjacent that Buffalo River and 
includes 488 miles. Using the Cattaraugus Creek gagings basis, the 
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discharge Buffalo River has been computed follows: 


Discharge, cubic feet 


Period per second 
Average from May October 31.......... 355 


Median from May October 


discharge 100 per sec exceeded the time for the year 
whole, and for the period from May October 31, discharge per 
sec exceeded 90% the time. 


Mar. 1937; 920 
per Sec 


Dissolved Oxygen 


Biochemical Oxygen 
Demand 


May 1937; 320 
per Sec 


July 1937; 160 


mo 


per Sec 


Dissolved Oxygen (Percentage of Saturation) 


geo 


5-Day Biochemical Oxygen Demand, in Parts per Million 


Nov. 1937; 450 


Distance From Mouth of River, in Miles 


Fig. shows the dissolved oxygen and the 5-day bio-chemical oxygen 
demand for five typical river surveys during 1937. The river discharges indi- 
cated were determined actual gagings the dates the surveys. During 


the surveys May, July, and September the dissolved oxygen was practically 
exhausted. 
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Scajaquada Creek.—F lowing through the northerly section the city into 
Black Rock Harbor, Scajaquada Creek has drainage area miles. 
Where flows through Delaware Park, has been dammed form small 
lake. Due the pollution waters the creek, sludge deposits this lake 
have caused objectionable odor conditions the past. The section the 
creek between the east city limits and Main Street has been enclosed 
large reinforced concrete drain. 


ENGINEERING AND LABORATORY 


Much field engineering work was required done promptly, including 
surveys for design and for the acquisition land and rights way, investiga- 
tions the best location for intercepting sewers, the supervision borings, 
examination and gaging existing sewers, and planning and organizing 
laboratory. Much useful work was accomplished through the laboratory, 
determining sewage quantities and characteristics and the condition the 
water many the waterways. Studies and reports were made most 
the larger industries. Investigation the industrial sewages Buffalo was 
major item the engineering and laboratory work. 

During two periods 1936 (May and August 22), and one 
1937 (July August 28), several the large sewers were gaged and sampled 
determine the dry-weather characteristics the sewage. The 1937 study 
was the most extensive and all the following data refer it. Each sewer 
was sampled for period days. Chlorine demand was determined 


BASED GAGINGS AND ANALYSES 


ANALYTICAL 
PER MILLION 
Aver- 


gallons 
de- 


Districts 


Bailey Avenue 
Babcock Street 
Fillmore 
Emslie Street 

Swan Trunk* 
Albany Street 

Bird Avenue 


Totalst 


Estimate 
entire city....... 573,000 110 


This sewer includes the flow from all the first four. The total includes only the last three items— 
see preceding footnote. part. 


separate samples collected bi-hourly intervals. Bacterial analyses were 
made two samples daily. All other determinations were made 24-hr 
composites, built from samples taken 15-min intervals and proportioned 
according sewage flow. Table gives the results this work. The sewage 


ten- 
tial 
Sewer (see hy- 
Fig. dro- 
gen, 
solids 
400 2,000 4.70 472 375 498 296 8.8 
32,000 6.10 639 186 153 238 180 7.0 
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quantities were determined from measuring the head weirs built the 
sewers. 

During 1936, investigations were made the effect storms the 
quantity suspended solids the sewage, with special attention the charac- 
teristics the first flushings, which would taken the treatment works. 
all, some 150 samples were collected, twelve studies. 

The effect storm runoff suspended solids indicated the study 
the Bird Avenue Sewer, results which are given Table 


Brrp AVENUE SEWER 


Time Sewage flow, solids, Time Sewage flow, solids, Time Sewage flow, solids, 

(p.m.) million parts per (p.m.) million parts per (p.m.) million 


million 


:00 
600 
300 
32.8 252 


Monthly surveys Niagara River have been made since June, 1936. 
During the summer months, when the use boat was possible, samples were 
taken six different cross sections the river between the site the treatment 
plant and the lower end Grand Island, immediately upstream from the 
intake the Niagara Falls water works. About forty stations were sampled 
these times. During winter months, samples were obtained only from the 
two Grand Island bridges, total thirteen stations. Analyses were made 
each sample determine the numbers coli and total bacteria 
agar-agar 37° few occasions, determinations were also made 
dissolved oxygen (D.O.), bio-chemical oxygen demand solids, nitro- 
gen, oxygen consumed, alkalinity, and chlorides. 

Daily bacterial analyses the raw water the Niagara Falls filter plant 
have been made for the Buffalo Sewer Authority since May, 1936. Two 
studies were also made determine the hourly variations bacterial content 
the river water this point. 

Monthly samples from about ten points along the Buffalo River and 
Cazenovia Creek have been collected and analyzed since April, 1936. The 
stream discharge and temperature are determined and the following analyses 
made: Potential hydrogen, dissolved oxygen, bio-chemical oxygen demand, 
total solids, suspended solids, oxygen consumed, total nitrogen, chlorides, 
alkalinity, 37° bacteria, and coli. Samples special locations and for 
special purposes are collected various times. 


iia 


anal 
Ave 
Bric 

dep 

gati 
var 

nat 

5:30 25.3 208 7:30 28.5 
6:30 120.6 964 7:50 204 
7:30 104.0 1,220 8:05 41.9 132 elg 
8:30 56.7 696 46.2 508 
9:30 37.6 364 8:30 46.6 704 

10:30 28.6 172 8:45 43.1 844 

lik 
tic 
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each the Buffalo River surveys samples were also collected and 
analyzed from the following points the harbors: City Ship Canal Michigan 
Avenue, mouth Buffalo River, Black Rock Harbor Ferry Street and 
Bridge Street. 

Occasional studies have been undertaken determine conditions the 
Scajaquada Drain and Park Lake, with particular attention sludge 
deposits and the rate digestion. Many other more less extensive investi- 
gations have been made covering the chlorine demand the sewage, hourly 
variations sewage characteristics, analyses materials used construc- 
tion, 

INDUSTRIES 


Buffalo has relatively large number industries which produce sewage. 
Many these are situated the southerly part the city which drains 
naturally Buffalo River. Some the larger plants discharge directly into 
Buffalo River and not use the sewers. The chemical industries form 
important group, including seven major plants, among which are those the 
General Chemical Company and the National Analine and Chemical Company. 
There are two relatively large coal coking plants, one which has been closed. 
There are twenty-four meat packing plants, eleven breweries, fifteen dairies, 
eight oil handling and refining plants, five rubber reclaiming plants, and five 
tanneries. 

Many these industrial establishments have been surveyed and the 
quantities and characteristics the sewage determined gagings and analyses. 
has been considered that the sewage from some the industries not 
likely readily amenable treatment. has thus been desirable 
undertake cooperative program with the industries for the removal pollu- 
tion from the Buffalo River and from the sewers the city. This program 
has proceeded satisfactory basis. Properly, has not been developed 
until the later steps construction and the early steps operation the treat- 
ment plant. 

RELIEF SEWERS 


Storm-water relief sewers were planned serve the Bird Avenue sewer 
district which has area 5,490 acres. The main sewer for this district 
was inadequate throughout its length, especially the section above Main 
Street, which has capacities, its several sections, only 166 665 (ged), 
based estimated 1985 populations. Overflows had been built from this 
sewer into the Scajaquada Drain, large sewer which, more less, parallel 
the Ferry Street sewer the north. This caused the pollution Park 
Lake Delaware Park. Several lateral sewers this area were also 
insufficient capacity, especially that Michigan Avenue. 

Storm-water relief sewers were planned for rate rainfall estimated 
occur once ten years. runoff coefficient 0.60 was used for the area 
within the city limits and 0.22 for rural areas outside Buffalo. The 
size the necessary relief sewers ranged from 3.5 12.0 diameter. The 
total cost the storm-water relief sewers was $2,330,000. Plans and specifica- 
tions for this part the project were prepared the City Sewer Department. 
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INTERCEPTING SEWERS 


Intercepting sewers were designed serve future population 1,100,000, 
The average quantity sewage each the ten major sewer districts (see 
Fig. was estimated shown Table The unit quantities used com- 
puting these sewage flows are follows: 


Unit rates average 
sewage flow 


Gallons per capita per hr; domestic............ 
Gallons per Acre per Hr: 
Industrial 
Infiltration 


summary the estimated average sewage quantities, million gallons 
daily, follows: 


Description 1930 1985 
Totals, gallons per’ 
capita per day...... 213 186 


For the design sewers and structures serving only part the city, the 
estimated average sewage quantities were increased factors allow for 
different distribution the future growth than that estimated. These factors 
varied with the size the area served and with the estimated density popu- 
lation, having maximum value 1.5. The factors were applied the 
domestic, industrial, and commercial quantities, but not the infiltration. 

For area that will overflow the Buffalo River the intercepting sewer 
capacities were based 3.0 times the estimated 1985 average flow, and for all 
other areas 2.5 times the estimated 1985 average flow. 

Based rainfall records and runoff coefficient 60%, computation 
was made indicate approximately the number and duration overflows 
combined sewage into the waterways. very exact computation was 
attempted, but the approximate results are given Table proportion 
the total sewage flow varying from about one third one per cent, the near 
future, almost one and one third per cent 1985 will overflowed during 
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rain storms. Kutter’s formula, with 0.013, was used determining the 
capacity all the intercepting sewers. 

Careful studies were made determine the economical balance between the 
size and slope the intercepting sewer and pumping costs. order utilize 


TABLE FREQUENCY AND DURATION OVERFLOW FROM 
INTERCEPTING SEWERS FOR THE SuMMER 
May 


Inter- Dry- Dry- Dry- 
cepting FLow; WEATHER WEATHER 
Population 150 GALLONS PER 200 GALLONS PER 250 GALLONS PER 
capacity, density PER Day PER Day Day 
gallons per acre 
per capita 
overflows overflow overflows overflow overflows overflow 
400 230 104 257 
600 168 181 195 
800 132 140 147 
1,000 106 111 116 


all the head available and the same time keep the sewer below existing and 
future utilities, much the main intercepting sewer was built flatter grade 
than the slope the hydraulic gradient maximum flow, that times 
storm the 11-ft and 9-ft portions the sewer will operate with the hydraulic 
gradient much 2.5 above the crown the sewer. 

condensed profile the main intercepting sewer shown Fig. and 
typical sections the intercépting sewer are shown the 
profile, the capacity the sewer varies from 367 1,047 ged, with the higher 
capacities along the Buffalo River. 

For sewers between 3.5 and diameter, bidders were allowed 
choose any one four sections: Monolithic concrete, reinforced concrete pipe, 
brick, and segmental tile. Reinforced concrete pipe (types and Fig. 
was the choice bidders every case, except for monolithic concrete (type 
for earth tunnel division Fig. 

For rock tunnels the horseshoe shape (type was used, with plain 
concrete lining having minimum thickness For the larger sewers 
open cut, the semi-elliptical section was chosen after extensive investigations 
the relative economy several sections. Drawings types and 
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Fig. show the sizes and dimensions used. sections were used for 
shallow cuts where extraordinary loading was anticipated. 
sections were used for deep cuts where other heavy loads were anticipated, 
such from railroad tracks. The entire length the 11.5-ft, type section 
the abandoned Erie Canal was designed for heavy loads provide for 
proposed relocation main-line tracks the New York Central Railroad. 


Reinforced 


Concrete Pipe 
N 


RTH TUNNEL ROCK TUNNEL 

Thicknesses Shown 
are Minimum 


ROCK EARTH ROCK EARTH 
FOUNDATION FOUNDATION FOUNDATION FOUNDATION 
STANDARD SECTION HEAVY SECTION 


TUNNEL 


Many special structures were required. Near Bailey Avenue, the 
southeastern part the city, the South Buffalo pumping station, equipped 
with four vertical centrifugal pumps having total rated capacity 100 mgd. 
tunnel crossing under the Buffalo River leads this pumping station. 
comprises two conduits, one in. diameter and the other, rectangular 
section, 6.5 ft. The main crossing under Black Rock Harbor the 


o 4" 4" 3" o" 
TYPE 
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Sewage Treatment Plant comprises two conduits, each diameter. Some 
intercepting and regulating chambers were built, well large number 
special structures including creek and sewer crossings, junction chambers, and 
several lateral separate sewers. 

The estimated final construction cost the intercepting sewers (exclusive 
unsettled claims) may summarized follows: 


Main intercepting sewer............. $5,254,000 
South Buffalo pumping station....... 223,000 


The cost the Main Intercepting Sewer made $4,703,000 for the 
sewer line proper, including manholes, sheeting left place, rock excavation, 
paving, handling utilities, etc.; and $551,000 for special structures, such 
intercepting chambers, junction chambers, and river and creek crossings. The 
sum $4,703,000 itemized, shown Table 


CLUDING MANHOLES, SHEETING, Rock, PAVEMENT, 


Con- Cost, Con- Cost, 
tract Type Size Length, dollars, tract Type Size Length, dollars, 
Division (see linear per Division (see linear per 
(see Fig. feet linear (see Fig. feet linear 
Fig. foot Fig. foot 
(1) (2) (3) (4) (5) (1) (2) (3) (4) 
6.5 2,059 51.70 11.0 5,074 145.30 
7.0 4,578 41.60 11.0 1,273 148.75 
6.5 1,049 132.10 11.0 1,176 120.40 
6.5 1,407 93.00 11.0 490 252.20 
8.0 1,877 104.20 11.0 1,410 260.00 
11.0 147.10 11.0 1,989 92.70 
9.0 5,720 97.90 11.0 2,864 73.60 
11.5 8,582 80.50 11.0 105 105.80 
6.5 495 58.60 11.0 135.90 
3.5 9,680 28.90 11.0 370 136.60 
6.5 1,375 57.90 


Standard section, rock foundation. section, earth foundation. Heavy section, rock 
foundation. Heavy section, earth foundation. 


STATION AND SEWAGE TREATMENT PLANT 


special permission the City Buffalo, the State New York, and the 
United States, the sewage treatment plant situated Bird Island, along the 
west boundary the city and somewhat north the east and west center 
line. The major bases design are follows: 


Bases Design 
Population: 
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Sewage Quantities (All Million Gallons Daily): 
Capacity intercepting sewer (1985) 
Capacity main pumping station 
1950 basis (population 750,000) 
Average 
Maximum: 
Clarification and disinfection units........... 
Conduits 
4-hr maximum 
24-hr maximum 


Minimum (for less than one hour) 
Near future basis (population 640 000): 
Average 
Minimum (for less than one hour) 
Sewage Characteristics: 
Suspended Solids Raw Sewage Average (1950): 
Parts per million 
Pounds per capita per 
Suspended Solids Removed and Disposed of: 


Period Percentage removal Pounds per 


Average 108,000 
Maximum 


448,000 

351,000 

278,000 

135,000 
(Assume that, yearly average, the volatile matter 70% 
the total suspended solids and that times 

maximum flows 50%.) 
Quantities Grit: 

Cubic yards per Cubic yards 


Period million gallons per 
Average 22.5 


Maximum 


Grit Chambers: 

Basis design, million gallons daily 

Number units channels 

Average depth sewage, feet 

Length, feet 

Width, feet 

Velocity, feet per second: 
Maximum 
Minimum 


1357 
563 
150 
570 
570 
570 
ARN 
300 
200 
’ 
128 
180 
0.3 
200.0 
50.0 
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Coarse trash racks with bar-spacing 5.5 in. 
Bar screens mechanically cleaned with clear spacing bars 
0.75 in. 
Sedimentation Tanks: 


Displacement period average flow of: 


Raw Sludge Quantities: 
Other quantities: 


Gallons wet Pounds 
Period sludge per solids per 

Maximum 

1,053,000 351,000 

278,000 

Maximum month....... 405,000 135,000 


(Assumed weight wet sludge, 62.5 per ft) 


Chemicals for conditioning raw sludge, pounds per 100 
dry solids: 


Average Maximum 
Ferric chloride (100%)........ 2.0 3.0 


general plan the sewage treatment plant shown Fig. and profile 
showing the high-water and low-water hydraulic grade lines shown Fig. 
The elevations Fig. are referred Buffalo City Datum, which 575.89 
above mean tide New York, (1903 level adjustment, Army 
Engineers). The gradients show that, 


(a) With river level 5.7 and with sewage flow 570 mgd, 330 mgd 
will discharged through the submerged outlet and 240 mgd through the 
overflow weir; 

(b) With river level 9.7 and with sewage flow 570 mgd, 500 mgd 
will discharged through the submerged outlet and mgd through the 
overflow weir; 

(c) With river level 5.7, sewage flows much 225 mgd will 
discharged through the submerged outlet, without overflow; and, 

(d) With river level 9.7, sewage flows much 470 mgd will 
discharged through the submerged outlet, without 


The treatment plant occupies somewhat triangular shaped site, having 
area 23.0 acres. the southerly end the main pumping station which 
lifts the sewage about average sewage flow. The sewage then flows 
gravity through screen and grit chamber, venturi meters, aerated channels, 
four sedimentation tanks, and submerged outfall the Niagara River. 
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The pumping building contains six vertical, close-coupled, centrifugal 
pumps, each with rated capacity 100 mgd and low-head capacity 120 
mgd. The total capacity the pumping station thus somewhat more than 
that the intercepting sewer system. Two the pumps are fitted with two- 
speed motors assist uniform regulation flow into the sedimentation tanks. 

There are twelve grit chambers. pair operates unit with one bar 
screen and pump, and there are cross connections between these units. 
Sereenings are discharged into grinders and returned thence the sewage. 
The grit chambers proper are long and wide; and they have 
average liquid depth ft. 

After leaving the grit chambers, sewage flows through two conduits under 
the main building, each conduit having venturi tube. Following the venturi 
tubes, the two conduits are equipped with plates and piping for aerating the 
sewage. This section the conduits 285 long and traverses part the 
site left open, for the present, provide for future structures. Chlorine 
added the sewage the far end these aerated conduits near the sedi- 
mentation tanks. The chlorination equipment housed the main building. 

The sewage then flows through by-pass and control chamber four 
circular sedimentation tanks which provide displacement period min 
rated capacity. Each these sedimentation tanks 160 diameter, and 
the side-water depth 14.5 ft. Effluent from the sedimentation tanks flows 
into the Niagara River about 200 from shore, depth about ft. 

The suspended solids are pumped from the sedimentation tanks four 
sludge storage tanks, having total capacity 750,000 ft, 1.0 per 
capita. Each tank diameter, with side and center depths 
and ft, respectively. the tanks equipped with floating cover. 

the center the space occupied the sludge storage tanks building 
that houses sludge pumping and control equipment with provision for the 
delivery sludge given rates into one all the four tanks, the transfer 
sludge from one tank another, and the sludge-disposal building (see 
Fig. 5). 

The sludge disposal building contains three de-watering filters the vacuum 
type, three flash-type sludge incinerators, and the boiler plant for heating the 
sludge storage tanks and all buildings. Each drum filter 11.5 diameter 
and long, with filtering area 500 ft. estimated that these 
filters will de-water the sludge moisture content 70%, the average rate 
5.0 dry solids per per when conditioned with about 2.0% ferric 
chloride and 8.0% lime (expressed 

Sludge cake from the filters discharged belt conveyer extending the 
length the building which, turn, discharges cross conveyer where the 
sludge weighed and elevated another longitudinal conveyer which delivers 
the sludge incinerator units. 

Sludge cake, mixed with previously dried sludge, delivered flash dryers 
the sludge incinerators; thence blown into the furnaces for combustion. 
Any deficiency heat units made burning oil gas. The sludge 
incinerators include ash-removal equipment, fly ash cyclones, and other 
appurtenances. Each the three incinerators has rated capacity 100 
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tons sludge cake (containing 70% moisture) per hr. Two the inciner. 
ation units will probably sufficient dispose the average quantity 
sludge. 

The various features the sewage treatment plant are connected 
service tunnel along one side which the main ventilating duct constructed, 
The several parts and buildings the plant where sewage exposed are 
provided with large ventilating fans which withdraw air and discharge into 
this main ventilating duct. The chimney the treatment plant, rising 175 
above the ground, actually two concentric ducts. The inner these for 
the hot gases from the sludge incinerators and the boiler plant, and the outer, 
annular one, for ventilation and receives the air blown through the main 
ventilating duct. 

The main building houses the office and laboratory, garage, machine shop, 
blowers for the aerated conduits, chlorine equipment, and the plant water 
supply, including pumps and elevated tank. flow diagram the sewage 
treatment plant shown Fig. 7(a), and sludge incineration unit 
Fig. 

The approximate total cost the main sewage pumping station and sewage 
treatment plant was $4,378,891, distributed follows: 


Contract division Total cost 
Pumping station and grit chamber............ 914,429 
324,000 
Sludge de-watering equipment................ 74,965 
Incineration equipment...................... 257,500 
Sludge disposal building..................... 298,670 
Boiler plant 69,314 
Main building, service tunnels, sewage conduits, 

Yard water 22,327 
Wall around sedimentation tanks............. 16,670 
Miscellaneous informal bids.................. 25,000 


The part this total sum required for sewage treatment about $3,500,000, 
which equivalent $4.67 per capita $23,300 per mgd design capacity. 
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summary the estimated final costs for the major divisions the Sewage 
treatment plant construction shown For the most part, bids were 
received unit price basis. Table shows the high and the low bids the 
successful bidders and weighted average for each kind construction for 
which unit price bids were received. 


TREATMENT PLANT 


Description Weighted 
average 


excavation, cubic yards 
Rock excavation, yards 
Selected fill, yards 
Embankment, cubic yards 
Sheeting Left Place: 
Wood, thousand feet, board measure 
Steel, pounds 
Sheet Piling: 
Wood, thousand feet, board measure 
Steel, pounds 
Timber piles, linear feet 
Steel H-piles, linear feet 
Steel reinforcement, pounds 
Concrete, Cubic Yards: 
Class A-1, formed floors and footings 
Class A-2, formed floors and footings 
Class fill 
Rubbed finish, square feet 
Expansion Joints: 
Copper plate, linear feet 
Mastic, square feet 
Asphalt, square yards............. 
Common brick masonry, cubic yards 
Cast-Iron Pipe, Pounds: 
Bell-and-spigot, 
Bell-and-spigot, specials 
Flanged, 
Flanged, specials 
Structural Metal, Pounds: 
Structural steel 
Miscellaneous steel and iron 
castings 


au 
to 


bo 


to 


oo 


no 


Pipe railing, linear, 
Galvanized floor grating, square feet 


REVENUES AND EXPENSES 


Under the act creating it, the Buffalo Sewer Authority has taken over the 
maintenance and operation all the sewers the city. Thus, estimates 
annual operating costs include items for the maintenance lateral sewers, for 
sewer cleaning and repairs, and for progressive rehabilitation. round nun- 
bers, the annual budget summarized follows: 

Item Annual sum 
Administration and expense main office 80,000 
Collecting sewers 223,000 
Outlying pumping stations 59,000 
Sewage treatment plant 397,000 
Construction work, insurance, and miscellaneous 141,000 
Debt service 600,000 


$1,500,000 


re’ 

Item 
No. 
13.28 
52.67 
0.0258 
93.34 
0.0556 
5 .0488 

20.71 

0.296 fir 

34.33 

0.1143 
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The Authority has adopted two-part rate schedule for raising the annual 
revenue. This rate schedule provides that 55% the total sum (or $825,000) 
shall derived from charge based water consumption; and that the 
remainder (or $675,000 per yr) shall derived from charge paid all 
taxable property proportion the assessed valuation. The portion the 
rate applicable water use 40% the water rates, and the portion applicable 
assessed valuation about $0.73 per $1,000. 
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IMPROVED METHOD FOR ADJUSTING 
LEVEL AND TRAVERSE SURVEYS 


CLARENCE AND JULIUS SPEERT, 
Assoc. Am. Soc. 


simplified method for adjusting level and traverse surveys presented 
this paper. The results are identical with those obtained least squares, 
but the method simple understandable any one possessing 
elementary knowledge algebra and surveying. The labor required for 
this method adjustment considerably less than that required the 
usual least-squares methods, and greater than that cruder methods. 

The paper divided into three parts: The main body gives detailed in- 
structions for the adjustment, such form that the computations may 
performed clerk familiar with the use computing machine. For 
the benefit those interested, the theory this method developed 


Appendix and Appendix explains simple method solving the adjust- 
ment equations. 


INTRODUCTION 


All surveys, whether for the determination altitudes relative 
horizontal positions, are subject two general types error—systematic and 
accidental. Before any adjustment survey attempted, all known 
systematic errors should removed, leaving only those the accidental 
type disposed the adjustment. 

Frequent azimuth checks are usually made well-conducted traverse 
survey. preliminary angle adjustment can applied, therefore, before 
the computation latitudes and departures undertaken. The angle adjust- 
ment simple and well understood that will not discussed herein. 

When latitudes and departures—or differences latitude and longitude— 
have been computed for each course traverse net, some error will usually 


appear the form imperfect closures the net. This error due two 
comments are invited for immediate publication; insure publication the last dis- 

cussion should submitted February 15, 1940. 

Asst. Topographic Engr., Geological Survey, Washington, 

Asst. Topographic Engr., Geological Survey, 
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general causes—accidental errors the straight-line measurements (including 
other errors treated accidental), and errors that still remain 
after the preliminary adjustment the angles. Because the complex 
interrelationship these errors, and the practical difficulties involved their 
separation, customary, except the most precise surveys, treat the 
traverse closure single accidental error with two components; namely, 
northing and easting (or latitude and longitude). this paper assumed 
that all preliminary computations have been performed and that unadjusted 
north and east coordinates have been determined for both extremities each 
line the net. The resultant closure errors, northing and easting, 
will treated independent items, either which may adjusted without 
regard for the other. 

Similarly, level net, the error will appear imperfect closures, 
differences elevation, the junction points. The method adjustment 
described applicable equally level traverse net; but, because 
the latter presents more general problem, having two errors adjusted, 
special attention will given the adjustment traverse net. 


ADJUSTMENT TRAVERSE NET 


The method adjustment described herein based the so-called 
“Junction Point method” for the adjustment levels developed Mr. 
Baldwin the Geological The theory the new method 
presented Appendix this paper. 

illustrate the method and its application traverse adjustment, 
somewhat involved traverse net will adjusted detail. This net shown 
diagrammatically Fig. with the junctions numbered for identification. 
The stations whose positions are held fixed result earlier adjustment 
more precise survey are indicated Fig. and are listed Table 
with their fixed positions. The results the preliminary computations 
north and east coordinates are listed Table the upper entry representing 


Station 


(see Fig. 
North East 


248,367.22 126,478.29 
328,946.18 129,817.26 
315,406.20 154,926.18 
319,029.42 240,792.22 
250,612.16 232,163.26 
293,936.85 187,391.63 


the beginning, and the lower entry the end, the line indicated Column (1), 
accordance with the direction arrow the line Fig. The items 
For complete description the Baldwin method see Instructions the United States 


Geological Survey,” Birdseye, Am. Soc. E., Bulletin No. 788, Geological Survey, 
pp. 
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Columns (4) and (5), Table are for the computer’s reference the 
computation books. Column (7) identifies the stations the extremities 
the various lines, shown Fig. The positions shown Tables and 
have been computed plane coordinates. However, this method may 
used equally well for adjusting traverse geodetic coordinates, differences 
latitude and longitude being substituted for differences northing and 
easting the adjustment computations. 

Special few special and conditions are included 
herein show how such conditions may treated adjustment. will 
noted that line Fig. crosses line without tying it. The omission 
junction this point merely eliminates one equation from the adjustment 
and presents special problem. Although lines and were surveyed 
independently the field, they are reality single line spite the closure 
between them, and they are treated one line the adjustment, the weight 
being that the two lines combined. The same treatment could applied 
lines and but, illustrate the method, these are treated separate 
lines, closing junction indicated. 

Junction junction closure the distance which the com- 
puted position point one line differs from the computed position the 
same point another line. positive the computed position the 
closing line the greater, and negative the smaller. Stated otherwise, 
the junction closure the distance that must subtracted algebraically from 
the position the closing line obtain the corresponding position the line 
station closed upon. junction closure indicated short arrow 
parallel the closing line, pointing toward the line closed upon, and 
numbers representing the magnitude the break closure. The junction 
closures indicated Fig. are derived from inspection Tables and 
For example, junction the end line and the beginnings lines 
(see Table 2); therefore, zero closures are indicated between line and lines 
and The absence closure arrow the extremity any line implies 
zero closure for that line the junction involved. junction station 
has the same position lines but line its position 5.37 less 
northing and 2.59 less easting. Therefore, line closes junction with 
closures 5.37 and 2.59, asindicated. The position station lines 
4.02 greater northing and 0.03 less easting than the fixed position 
station (Tables and 2), indicated the closures these lines the 
fixed station. Fig. starts station with the correct final position, 
and its closure accordingly, zero. All closures, junctions and 
fixed stations, are determined the same manner. When dealing with items 
northing and easting, desirable always write the northing quantity 
first, thereby avoiding confusion well the necessity labeling each item. 

Circuit the next preliminary step the adjustment, 
desirable investigate the circuit closures. Although these not enter 
directly into the adjustment, they offer key the quality the field work 
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and preliminary computations, indicating places apparent weakness and 
frequently aiding the discovery blunders. 

The circuit closures may computed adding all junction closures 
counterclockwise around each circuit, changing the sign junction closure 
its closure arrow opposed the counterclockwise route around the circuit. 
For this purpose, group lines connecting two fixed stations treated the 
same closed circuit. Fig. the circuit closures are shown within 
the respective circuits, inside curved arrow which indicates the direction 
closure. Each circuit closure, and all the junction closures contained it, 
may checked from Tables and using Fig. guide. Proceeding 
counterclockwise around the circuit before, subtract the position the 
beginning, and add the position the end, each line whose direction arrow 
agrees with the direction progress around the circuit; reverse the operations, 
adding the beginning and subtracting the end, for each line whose direction 
arrow against the direction progress around the circuit. circuit 
connects two fixed stations, assume imaginary line joining them, and treat 
the fixed positions the extremities this imaginary line. The circuit 
closures obtained this method should check exactly those determined 
previously. Any discrepancy might due mistake the determination 
the circuit closure either method, mistake one more junction 
closures, combination these mistakes. 

The magnitude any circuit closure relation the length its circuit 
isa measure the precision the survey around that circuit. excessively 
large closure might clearly indicate the existence gross blunder and lead 
its detection. Large closures opposite signs two adjoining circuits 
indicate weakness, possibly blunder, the line common these circuits. 
general, the algebraic difference between the circuit closures opposite 
sides any line may used form rough preliminary estimate the 
correction that the adjustment may expected force into that line. Thus, 
any line that appears below the standard the remainder the work 
may reasonably omitted from the adjustment, order that its weakness 
shall not detract from the strength the remainder the survey. 

the adjustment level net, the lines are weighted inversely 
their relative lengths. For traverse, however, believed that the 
number transit stations line offers better index the strength that 
line. For convenience solving the adjustment equations, avoid coeffi- 
cients either too large too small, desirable that the magnitude the 
diagonal terms the equations between and 10. This condition can 
produced most readily using weights for the lines, such that their average 
value will about Therefore, the effective length of, number stations 
in, each line should divided the approximate average all the lengths, 
numbers stations, before weighting. this adjustment the lines will 
weighted inverse ratio the number stations the weight 
each line this example will taken times the reciprocal the number 
stations it. The number stations each line listed Table 
and the weights are marked the lines Fig. 
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the preliminary positions and the closures have been computed two 
decimal places, desirable that the weights used one place more, 
three places, order that there shall loss accuracy due uncertainty 


TABLE EQUATIONS 


—3.513 |+1.042 +9.134 


—2.737 +0.833 +0.367 


3.741 
(—0.749) 


the last decimal figure. any line has been surveyed with greater less 
care than the others, entitle special consideration, its weight factor 
may changed accordingly. this example all lines are considered equal 
strength. 

Adjustment Form.—After the preliminary examination and preparation are 
completed, the may begun. All junctions between lines the 
survey should numbered, except where lines meet fixed station. Thus, 
there should number for each junction for which correction 
determined. The form for the table equations (Table should prepared, 
allowing one column for each numbered junction, two absolute columns, and 
additional column each side the table for later use. Four lines should 
allowed for each equation except the first, which requires only three. These 
groups lines should numbered the same order the column headings, 
and the first three lines each group should designated and respec- 
tively. The fourth line left blank separation between the equations. 


c 


METHOD FOR ADJUSTING SURVEYS 1373 


the first equation lines and are combined. further guide, hori- 
line may ruled above each b-line. 
adjustment equations may now formed inspection 


TABLE ADJUSTMENT EQUATIONS 


4.865 |+0.367 4.498 


+0.862 3.206 |+3.076 0.863 
+0.3550 0.3554 


—3.805 0.000 0.000 2.174 


—3.499 |+0.769 1.138 |+1.092 2.480 


+1.250 8.867 |—6.322 


+1.315 |+0.595 9.117 |—6.082 


+15.199 |+6.563 |+19.323 


(—1.0) |+0.0776 3.1336 


+0.625 0.893 |—16.506 |+9.963 7.237 


+0.625 0.893 |—17.689 |+9.047 9.782 
3.6788 


1.104 |—8.959 |—10.798 


Back solution 


3.741 


from For this purpose, well for the solution, computing machine 
almost indispensable. The original equations are written the a-lines, 
the other lines being reserved for their solution. The diagonal term for each 
equation entered the line and column bearing the number the corre- 
sponding junction, and always negative. The coefficient each the 
remaining terms equation entered the right the diagonal terms, 
the line and the column corresponding the extremities the corre- 
sponding line the survey, will illustrated. 

Diagonal Terms.—The diagonal term for any equation the negative sum 
the weights all lines entering the correspondingly numbered junction. 
For example, junction Fig. itis (1.042 1.562 0.909) 3.513, 
entered Column (1), line la, Table For junction (0.862 
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1.042 0.833) 2.737, entered Column (2), line 2a. Similarly, 
all diagonal terms are determined from Fig. and entered Table the 
a-line and column corresponding the number the junction involved. 

remaining coefficients Equation (1), Table are the 
weights the individual lines entering junction from the other numbered 
junctions. For example, line enters junction from junction and 1.042 
entered Column (2), line lines and enter junction from fixed 
stations and have coefficients the table. Similarly, line enters junction 
from junction and 0.862 entered Column (3), line line enters 
junction from junction and 0.833 entered Column (5), line 2a; 
line enters junction from junction and 1.042 found already 
entered Column (2), line la. not entered Column (1), line 2a, 
because that space the left the diagonal term, and only those spaces 
the right the diagonal terms are used. like manner, the coefficients 
for all the equations are entered Table 

Absolute absolute terms for northing and easting each 
equation may now computed and entered Columns and respectively. 
Each these terms the sum the weighted closures all lines entering 
its junction; the sum the products the weight each line entering 
the junction the sum the closures that line computed from the con- 
tiguous junction. essential that the sign and direction each closure 
observed carefully, the closure sign being held its direction toward the 
junction whose equation being determined, and reversed its direction 
from the junction. For example, for Equation (1) 1.042 (0) 1.562 
0.909 0.03) 9.134. Similarly, for Equation (7), 0.625 7.67 
3.68) 0.397 0.03 7.69) 0.746 (0) 7.754. The absolute terms 
should computed the same number decimal places the coefficients. 

When all the absolute terms have been computed and entered Table 
the adjustment equations are complete and may checked. The absolute 
terms should checked first, follows: Compute the weighted closures 
all lines into all the fixed stations the same manner which the absolute 
terms were computed; write the sum all these northing closures 3.741) 
the bottom the (N)-column and the sum all these easting closures 
2.992) the bottom the (E)-column. Including these items, the 
sum all entries each these columns should zero, very nearly 
small discrepancy the last place may due the rounding off residual 
decimals beyond the last significant figure. preferred, the (N)- and (E)- 
columns may checked simultaneously recording only the sum the 
northing and easting closures into the fixed stations 0.749), and adding 
this both columns equal zero. 

computing all sums products, the individual products need not 
recorded, but should accumulated the computing machine. For this 
purpose the decimal points should placed properly all the dials, and 
should observed carefully all operations. Only the keyboard and multi- 
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plier dial are cleared after each multiplication, the product being left the 
machine. Negative products are obtained with the reverse multiplier. The 
product dial cleared only when the desired sum has been recorded. 

(S)-Column.—The (S)-column the right Table offers further 
check the formation the equations, well step-by-step check 
the solution. The (S)-item for each equation determined adding the 
negative sum the weights all lines entering junction from fixed stations 
the entries the (N)-column and (E)-column. For example, junction 
lines and enter from fixed stations; therefore equals minus the sum 
the weights these lines plus plus (1.562 0.909) 5.593 9.134 
12.256. junction only line enters from fixed station; hence 
2.174 0.000 0.000 2.174. There are fixed stations con- 
tiguous junction therefore, equals the sum and 4.865 

When the (S)-column completed, the sum all entries line through 
and should equal exactly. The sum all entries Column (2) 
and line through and should equal Similarly, the sum 
Column (3) and line should equal 83. Thus, each equation checked 
adding its column down the diagonal term, then across through 
equal its (S)-entry. mistake any coefficient will appear discrepancy 
the check the two equations containing that term. discrepancy 
peculiar particular equation indicates mistake either its diagonal term 
its (S)-term. 

Solution equations Table may now solved 
the Doolittle method‘ shown Table The detailed solution especially 
adapted and simplified for use with computing machine explained with 
complete instructions Appendix II. The (N)-, (E)-, and (S)-columns are 
treated identically the forward solution, the check the (S)-column being 
applied the usual manner after each step. the back solution, inde- 
pendent computation made for CN, CE, and for each equation, using 
only the appropriate items for each, and omitting the other two quantities 
from each computation. The back solution checked, step step, 
adding CN, CE, and 1.0000 for each equation; the sum should equal 
for that equation, very nearly. 

(CN)- and (CE)-items the back solution are the 
required corrections applied the corresponding junctions Fig. 
These are rounded off two decimal places and are shown small rectangles 
with arrows pointing the line the extremity which the correction 
applied. should borne mind that these corrections apply the 
position closed upon each junction. The correction the extremity 
the closing line may determined subtracting the junction closure from 
the computed correction. Thus, junction the computed correction 
applies equally lines and junction the computed correction 


instructions for the solution are given several publications the Coast and Geodetic 
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applies lines and the northing correction the end line 
The corrections lines fixed stations are equal and opposite the closures 
these stations. The difference between the corrections the ends 
line represents the correction distributed prorated through the line, 

Final this point the adjustment the junctions complete, 
corrections having been determined for each extremity each line. One 
final check may applied, however, test the efficacy the adjustment. 
According the theory least squares, the most probable values for the 
corrections should satisfy the condition that the sum the weighted cor- 
rections applied all lines around given junction will zero. This test 
may applied around each numbered junction multiplying the weight 
each line entering the difference the corrections applied its ex- 
tremities. this purpose, the and E-corrections may tested together, 
the corrections the far end each line being applied with their proper signs, 
and those the near end with reversed signs. For example, junction 
0.909 4.02 0.03 1.30 3.13) view the number and 
magnitude the items included this test, the discrepancy 0.01 reason- 
able, and indicates satisfactory adjustment. The same test may applied 
all numbered junctions the net. using lines and check junction 
the break station should considered, thus: 0.397 4.02 0.03 
1.46 7.69 0.97 4.13) 0.746 (0.00 0.00 0.97 4.13) 0.625 
0.48 3.15 3.88 0.45) 0.00. this example, the discrepancies 
the eight numbered junctions vary from 0.02, with average value 
0.007. 

SuMMARY 

brief summary the operations required for complete adjustment may 
value the formation mental picture the entire Fora 
traverse net, the preliminary computations should completed and the 
and E-coordinates, the latitude and longitude, should determined for 
each junction. For level net the orthometric elevation each junction 
point should determined. Then the operations the adjustment should 
performed the following order: 


1.—Prepare diagram the survey, showing junction closures, circuit 
closures, and weights all lines. Number all junctions, except those fixed 
stations. 

2.—Prepare form similar Table for the adjustment equations. 

3.—Compute the diagonal terms the negative sum the weights all 
lines entering each numbered junction. 

the weight each line the space numbered according the 
junction numbers the ends that line, using only those spaces the right 
the diagonal terms. 

5.—Compute and enter the (N)- and (E)-columns, respectively, the 
weighted closures northing and easting all lines entering each numbered 
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junction. level adjustment, course, there would only one absolute 
column. 

6.—Check the (N)- and (E)-columns summing them zero with the 
weighted closures all lines entering all fixed stations. 

the for each equation the negative sum the 
weights all lines entering its junction from fixed stations, plus its (N)- and 
(E)-entries. Add each equation down its column its diagonal term, then 
across through its absolute terms equal its (S)-term, exactly. 

8.—Solve the equations the Doolittle method, forward and back solutions. 

9—On the adjustment diagram, enter the results the back solution 
corrections the lines closed upon the numbered junctions. 

10.—Subtract the junction closures from the computed corrections 
determine the corrections the ends the closing lines. 

the adjustment comparing with zero the sum the weighted 
corrections all the lines around each numbered junction. 

12.—Distribute prorate the difference between the end corrections 
each line the stations that line, complete the adjustment. 


believed that the results any survey, even poor one, can im- 
proved appropriate adjustment, that all parts the survey shall fit 
together with mathematical nicety after the computations are completed. 
other words, any survey worth making worth improving small additional 
cost careful adjustment. The relatively small cost adequately adjusting 
survey undoubtedly improves the results greater degree than could 
obtained similar expenditure for improved field technique. Obviously, 
poor work cannot made good adjustment, but the data resulting from 
any survey can improved value and usefulness all mathematical in- 
consistencies are removed. When the funds available for survey are limited, 
believed desirable reserve portion the money for the purpose 
adjustment. 

traverses, except for the most precise geodetic surveys, customary 
combine all accidental errors the closures latitude and longitude, 
latitudes and departures. this respect, practically all methods adjustment 
agree. Beyond this, variations exist, from adjustment “estimation” 
judgment adjustment least squares. For any adjustment, regardless 
method, desirable that diagram prepared show the various closures 
and the relative weights the respective lines. The diagram, although not 
essential part the adjustment, invaluable aid the determination 
and distribution the adjustment corrections. 

generally recognized that the method least squares offers the best 
adjustment possible accord with the theory probability. Although the 
method described herein not derived through consideration the squares 
residuals errors, based the principle that the correction any 
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junction shall weighted mean the corrections brought from the con. 
tiguous junctions (see Appendix I). This principle perfect harmony with 
the principle least squares, and the two methods give identical results, 

may emphasized, furthermore, that the adjustment equations that 
are set inspection from the adjustment diagram are identical every 
respect with the normal equations that would obtained the use the 
observation-equations method least squares. This new method 
definite advantage over the latter that the equations may written 
inspection, requiring auxiliary tables. has similar advantages over the 
condition-equations method, that table correlate equations not needed, 
and the solution the adjustment equations gives the desired corrections 
directly, without the necessity applying the correlates the correlate 
equations and computing the line corrections before the junction corrections 
can determined. 

believed that the method adjustment described herein simple 
easily understood any one familiar with the elements surveying 
and algebra. Although the Doolittle method suggested the shortest 
known method solving the type simultaneous equations obtained, these 
equations may solved equally well the usual algebraic methods. (For 
algebraic solution, the equations should written complete form; see 
Appendix I.) Although the method described identical with least squares, 
not expected have the psychological effect that the name squares” 
has upon most surveyors—the suggestion something too complicated 
understood except expert mathematician. Some technical understanding 
helpful, course, preparing and studying the closure diagram, but the 
mechanics the adjustment are such that, the instructions are observed 
carefully, complete adjustment can made clerk familiar with the use 
the computing machine. 

compared with cruder methods adjustment, this method requires 
little additional work for the theoretically most probable distribution the 
errors. The only work peculiar this method the writing and solving 
the simultaneous equations. this respect there compensating saving 
labor, particularly traverse adjustment, that corrections for both 
north and east coordinates are obtained simultaneously, whereas most 
methods these must adjusted separately. compared with the Baldwin 
method for level adjustment, this method gives more balanced adjustment, 
eliminates the need for exercise judgment the selection junctions, and 
permits adjustment entire network single set-up. The time required 
for adjusting ordinary level net this method compares favorably with 
that required the Baldwin method. 

Because this method suitable for the adjustment any type line 
survey, another possibility becomes apparent. many surveys happens 
that traverse and levels are run over the same routes for the purpose de- 
termining horizontal and vertical control for all monuments. the lines are 
run that the weights used for adjusting the levels are the same ratio 
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the weights used for adjusting the traverse, can made so, then 
both traverse and levels can adjusted the same time using three 
absolute columns the adjustment equations instead one two. With 
three absolute columns, the (S)-column should include the effect all three, 
and the solution may completed exactly the same manner. 
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APPENDIX 


THEORY 


necessary preliminary the development the new theory, the 
Baldwin method will reviewed 

the simple level net shown Fig. line runs from bench mark 
(the elevation which fixed) another fixed bench mark, closing 


Fixed 


Fixed 
Benchmarks 


ADJUSTMENT DIAGRAMS 


bench mark with error 0.03 indicated the small arrow 
point and the number 0.03 (this means that the elevation bench mark 
bench mark Similarly, line begins fixed bench mark and closes 
bench mark with error 0.05 (or the elevation bench mark 
determined line 0.05 less than the unadjusted elevation deter- 
are 10, and miles, respectively. 

For the adjustment, the lines are weighted inversely their lengths 
the ratio 2.00 1.00 1.25, indicated Fig. This adjustment 
based the principle that the adjusted elevation bench mark 
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weighted mean the elevations brought through the various lines tying it, 
This principle strict accord with the principle least squares and may 
readily derived from it. The weighted mean most conveniently obtained 
considering corrections and closures instead actual elevations. The 
tion line point (as shown small rectangle Fig. zero; the 
correction line point (to correct the closure 0.03 ft) must 
0.03; the correction line point Bis there are breaks 
continuity between bench mark and the indicated corrections points 
and the corrections these points may used for weighting point 
order refer common elevation point the closure 0.05 must 
added the correction point before the latter may weighted 
bench mark Then the weighted correction point may computed 
the sum the products the various line corrections the weights their 
respective lines, divided the sum the weights; that is, 


2.00 1.00 1.25 0.02 ft. 


express this principle general equation (see Fig. 3), let and 
represent the corrections points and respectively; ad, bd, and 
and respectively; and Paa, and represent the weights the 


Transposing all terms the right side the equation, 


Cee doses eon 


Note that represents the closure breaks the line whereas 
represents the same closure breaks the reverse direction; thus, 


Equation (10) represents general equation for the correction any 

tion terms the corrections other points and the mathematical breaks 
continuity the lines connecting the junction with the other points. 
narily, the Baldwin method has been used determine the correction 
particular junction when the corrections the other points have been known. 
Extension Baldwin Method.—The same principle may applied more 
complicated level net, indicated Fig. Denote the correction any 
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junction the letter used designate that junction; the weight the 
letter and the subscript letters the ends the line; and the mathematical 
continuity any line the lower-case letters the ends the line, 

and are fixed points, their corrections are known from the 
closures these points, and need not computed. Equations similar 
Equation (10) may written for each the remaining junctions. Thus, 


Poa (the sum the weights all lines entering junction 
designated and (ab) (db) (fb) (the sum the products 
the breaks the lines into the junction the weights these lines) desig- 
nated Equation (12) becomes 


Similarly, equations for the other junctions adjusted may written 
follows: 


and, 


there are six unknown junction corrections and six equations, rigid 
solution possible. 

Equations (13) (18) may written most conveniently inspection 
the form table, the unknown corrections being used designate the 
columns, which only the numerical coefficients need placed (see Table 5). 

the equations are written the same sequence the columns, will 
seen that the negative terms form diagonal, about which the table sym- 
except for the Q-column; that is, all coefficients the left (or the 
right) any negative term appear again above (or below) the same relative 
Thus, except for the Q-column, the table may read across 
down with the same results. This symmetry makes possible condense the 
table omitting all coefficients the left the terms. Then each 
equation read down its diagonal term and across through its absolute term. 
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this form the equations may solved readily the Doolittle 
give the required junction corrections (see Appendix II). 


TABLE EQuaTIONS 


(A) 
(E) 
+Pok 


This method adjustment not only suitable for levels, but even more 
advantageous its application traverse net. Instead single closure 
error each junction, traverse survey has two closures—in latitude and 
longitude, northing and easting. These two closures are considered inde- 
pendently and, other methods adjustment, require individual treatment. 
However, the method presented herein, the equations for the corrections 
any particular junction are identical except for the absolute Q-terms. 
Thus, both sets equations may combined single table, similar 
Table but having two absolute columns—one for northing and one for easting. 
The effort required for the solution both sets only slightly greater than 
that required for one alone. 

provide check the formation and solution the equations, extra 
column added the table. The items this column are dimensionless 
numbers, each being the sum the coefficients and absolute terms its equa- 
tion. This sum column not essential part the adjustment, but makes 
possible check each step before proceeding the next. 


APPENDIX 


SIMPLIFIED SOLUTION ADJUSTMENT EQUATIONS THE DOOLITTLE METHOD 


The Doolittle method semi-mechanical method for solving simultaneous 
equations the symmetrical type formed this method adjustment. 
much quicker than the regular algebraic solution, especially when there are 
many equations solved simultaneously, and particularly adaptable 
the use computing machine. this method, others, the equations 
are combined successively, one the unknowns being eliminated each step. 
When all the equations have been combined, the resulting equation contains 
only one unknown, which may evaluated readily. the back solution, the 
evaluated unknowns are substituted successively back into the combined equa- 
tions for the evaluation the preceding unknowns. 


tio 
ide 
cal 
tio 
of 
» 
Ce 
res 


1939 IMPROVED METHOD FOR ADJUSTING SURVEYS 1383 


Tables and the original equations have been written the a-lines, 
and each equation read down its column the diagonal term, then across 
through the absolute terms. The b-lines, Table contain the combined equa- 
tions, each which includes the influence those preceding it. Because the 
elimination one additional unknown, each combined equation contains one 
term less than the preceding one. The combined equations, the b-lines, are 
read straight across. The equations the c-lines are mathematically identical 
with the corresponding b-equations, but they express the principal unknown 
explicit function the remaining unknowns. The items the c-lines are 
the multipliers used the formation the combined equations the 
succeeding b-lines. obvious that the first equation the and b-lines are 
identical and need written only once. 

The procedure for combining and solving the equations follows few definite 
rules; when once mastered, becomes almost mechanical and requires practi- 
cally mental effort. The method best illustrated completing the solu- 
tion the equations Table The complete solution shown Table 

The first step the solution the formation Equation (1c) from Equa- 
tion (1b). Equation (1) represents the correction junction terms the 
all other junctions; therefore the entry Column (1) (the diag- 
onal term Equation (1a)) will called the principal term Equation (1b). 
Similarly, the b-item under each diagonal term will called the principal term 
the other corrections, its coefficient should made 1.0. This done 
multiplying all items line the negative reciprocal the principal term. 
The negative reciprocal (0.284657) 3.513 entered, for convenience, 
Column the left the table and multiplied the items line 
give the items shown line Line may now checked adding 
across through the absolute columns, including the first term, equal 
the (S)-item, very nearly. small difference may due the dropped 
residual decimals beyond the last significant figure. general, the c-lines 
should recorded one decimal place more than the b-lines, order retain 
the same number significant figures. 

The next step the solution the formation Equation (26) combining 
Equation (1b) with Equation (2a). this step, everything the left 
Column (2) ignored and each item Equation (1b) multiplied the 
Column (2); then added the (2a)-item its own column. 
Thus, Equation (2b) formed follows: 


1b-2(+ 1.042) 0.2966) 2a-2(— 2.737) 2b-2(— 2.428) 
1b-N(+ 5.593) 0.2966) 2a-N(— 4.865) 2b-N(— 3.206) 
1b-E(+ 0.2966) 2a-E(+ 0.367) 2b-E(+ 3.076) 
1b-S(+12.256) 0.2966) 2a-S(— 4.498) 0.863) 
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Each item designated its line and column, thus: Item line 2a, 
Column (3). Equation (2b) may now checked adding across equal 
its (S)-term, very nearly. Equation (2c) formed from Equation (2b) 
exactly the same manner Equation (1c) was formed from Equation 
checked adding across from 1.0 through the absolute columns 
equal its before. 

Equations (2b), and (3a) are now combined form Equation 
this step, everything the left Column ignored, and the c-items 
Column (3) are used multipliers for their respective equations. there 
item (giving zero multiplier), Equation (1b) drops out this step. 
Each item Equation multiplied item 2c-3 and then added the 
(3a)-item its own column form Equation 

Equation added across check its Equation (3c) 
formed from Equation (3b), and checked. This same routine followed 
throughout the forward solution. The first four equations are all short ones 
and are combined easily, explained. 

Equation formed combining Equations (1b), (2b), (3b), (4b), and 
(5a). this step, everything the left Column (5), Table ignored, 
and the c-items Column (5) are used multipliers for their respective 
b-equations. there item Equation (1b) drops out this step. 
compute item any column line multiply each b-item the 
column the c-item its equation Column (5), and add the products the 
(5a)-item above the term being computed, thus: 


always using the c-items Column (5) the multipliers for the corresponding 
b-items the other columns. Equation (5b) now checked adding across 
5b-S. Equation (5c) formed multiplying Equation (5b) the 
tive reciprocal its principal term, and checked adding across, including 
the 1.0, equal 5c-S, very nearly. Equation formed combining 
the previous with Equation (6a) already explained, and 
checked 

Back same routine followed the formation and check 
Equation (8c), which expresses the correction junction explicit 
function other unknown. other words, the corrections junction 
appear the absolute columns line 8c, and are 5.1507 northing and 
5.1115 easting. the first unknowns evaluated, these are entered 
lines and CE, respectively, Column (8) the back solution. provide 
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check, 8c-S entered CS-8. Then, Column (8), 1.0000 
should equal very nearly. (This check was already applied line 8c.) 

The corrections junction may evaluated next. Equation (7c) ex- 
presses each correction junction explicit function the corresponding 
correction junction which already known. The corrections junction 
(which will called, collectively, C8) are therefore substituted into Equation 
toevaluate C7. these corrections, reality, are independent variables, 
will necessary, the back solution, substitute them separately, and 
evaluate but one correction Thus: 


CE-8 7c-8 7c-E CE-7 
and, 


The (CS)-item not correction but used check the others and 
computed the same the check, before, CN-7 CE-7 
Equation (6c), and checked adding CN, CE, and 1.0000 equal CS, 
very nearly. Thus: 


CE-7 6c-7 CE-8 6c-8 6c-E CE-6 
and, 


like manner, all evaluated corrections are substituted successively into 
the next preceding c-equation for the evaluation the corrections the next 
preceding junction. The same check applied all; namely, 
1.0000 CS. Note that Equation (5c) there are items Columns 
(7) and (8); and C8, therefore, not enter into the computation C5, 
their multipliers are zero. Similar conditions are found throughout the 
remainder the back solution, materially reducing the labor required for this 
part the problem. When all the corrections have been computed and 
checked, back they are entered the diagram Fig. already 
explained. 


ep. 

the 

ved 

nes 

and 

ep. 

the 

ling 

ross 

ling 

licit 

and 

vide 


Py 
=? 
p 
0 
e 
4 
af nN 


AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


PAPERS 


RELATION THE STATISTICAL THEORY 
TURBULENCE HYDRAULICS 


The basic principles the statistical theory fluid turbulence are outlined 
this paper. Only those ideas and conceptions importance hydraulics 
problems are presented; the more complicated mathematical details are 
omitted. Turbulence can characterized intensity parameters which are 
the root mean square values the fluctuating velocity components, and 
length parameter proportional the most probable eddy size. The specific 
parameters used characterizing the energy dissipation and the diffusion 
power turbulence are derived, and experimental techniques and data are 
presented illustrate how these parameters are obtained. The two major 
problems hydraulics which knowledge the turbulence mechanism 


importance are those involving the increase decrease energy dissipation, 
and the transportation sedimentary materials. The practical application 
the theories turbulence some specific problems indicated. The 
experimental data were obtained injecting color streams and suspended 
immiscible liquid droplets into turbulent water streams and using motion- 
picture photography. 


PURPOSE PAPER 


Since about 1930, both Europe and the United States, considerable 
research fluid flow has been concerned with the fundamental aspects the 
mechanism fluid turbulence. aeronautics, wind-tunnel and atmospheric 
turbulence has been important field research since about 1925. Some 
the findings the aeronautics engineers and physicists have been permeating 
gradually into the hydraulics literature. Certain fundamental investigations 
water turbulence have been conducted, and more and more are being started. 
However, all this talk and writing about turbulence generally leaves the 
practical hydraulic engineer quite uninterested. The practical usefulness, 


comments are invited for immediate publication; insure publication the last dis- 
cussion should submitted February 15, 1940. 


Assistant Prof., Iowa Inst. Hydraulic Research, State Univ. Iowa, City, Iowa. 
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either immediate future, knowledge relating turbulence does not seem 
very apparent. 

the purpose this paper correlate the knowledge fluid turbulence 
with the solution some practical hydraulics problems, and present the 
information fluid turbulence, already existence, more simplified 
Some experimental data obtained the writer and others will presented, 
and various experimental techniques for obtaining pertinent information 
liquid turbulence will discussed. 


TuRBULENCE? 


Turbulence, present wind-tunnel air streams, has been studied for 
some time and has gradually attracted the attention hydraulics engineers, 
because there are many hydraulics problems which would better understood 
more were known the mechanism turbulence. 

When observing turbulence injecting dye some other material into 
water, when observing smoke coming from chimney, most haphazard 
mixing process seen, seemingly unpredictable nature. Fundamentally, 
may said that fluid any point state turbulence the direction 
and magnitude its velocity vector vary with time. This variation 
relatively rapid and general cannot predicted except the probability 
sense. Turbulence may caused result high shearing stresses set 
through screen honeycomb, some means such mechanical stirring. 
Confusion has arisen the minds some regarding whether not there 
definite periodic relationship between velocity component single point 
with respect True turbulence entirely random nature without 
any kind regular periodicity. stream discrete vortices which form 
directly behind body immersed moving fluid, such the von 
vortex street, not true turbulence; referred vortex motion, which 
breaks down eventually into true turbulence. 

Turbulence manifests itself the form large and small eddies, whirling 
about from one layer fluid another. These eddies can transfer momentum, 
mass, heat, etc., from layer layer, and can also cause the high dissipation 
energy that characterizes turbulent flow. apparent then that, charac- 
terize various quantities regarding turbulence (such as, for instance, the size 
the eddies, the velocity point, the frequency velocity fluctuation), 
average quantities must dealt with and the principles statistics employed 
determining and using these quantities. Additional fundamental ideas 
regarding turbulence can found the writings Dryden,? Hunter 
Rouse,? Assoc. Am. Soc. E., and Boris Am. 


Intensity and Scale Wind-Tunnel Turbulence,” Dryden, Report No. 
581, National Advisory Comm. for Aeronautics, 1937. 


Conceptions the Mechanics Fluid Turbulence,’’ Hunter Rouse, 
Am. Soc. E., Vol. 102 (1937), 463. 


Turbulent Flow,” Boris Bakhmeteff, Princeton University Press, 1936. 
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PRINCIPLES STATISTICAL THEORY TURBULENCE 


the outset let noted that complete and mathematically rigorous 
theory turbulence has not yet been developed. However, far the 
hydraulics engineer concerned, such theory, even available, would 
little use. mathematical difficulties are great that the development 
rigorous analysis has taxed some the best mathematical minds. The various 
ideas that have helped greatly understanding the mechanism turbulence 
Am. The statistical principles they have set forth have 
been checked experimentally, and are fundamental importance under- 
standing the turbulence mechanism. 

convenient represent the varying velocity vector any point 
turbulent flow means the velocity components along three rectangular 
axes. The velocity the direction mean flow, which taken along the 
z-axis, any instant represented u), which the mean velocity 
that direction and the fluctuation that direction. Obviously, the 
arithmetic mean zero. The velocity along the other two axes, and 
the x-direction, and are zero. Since and are quantities that 
vary with time, useful represent them certain statistical averages. 
Before this done, the statistical distribution these quantities should 
investigated. what the form, for instance, which function 
such that indicates the probability that the value lies between 
dv? Such function, f(v), called statistical distribution function. 
the value could determined any point for sufficiently long period 
time the general form f(v) could found. 

The value the velocity component with respect time shown 
These data were obtained turbulent water stream point near 
the center 12-in. 12-in. open channel where the mean velocity was 7.8 
in. per sec (see heading Technique”). block diagram indi- 
cating the frequency occurrence various values shown Fig. 
The smooth curve Fig. has the equation, 


Fig. bar over term terms indicates mean taken with respect 
time. 

This the ordinary normal error law. Other investigators* have found the 
same type distribution law, thus indicating the random nature the velocity 
935, Theory Turbulence, ‘‘by Taylor, Proceedings, Royal Soc. London, Vol. 151A, 


von Journal, Royal Aeronautical Soc., Vol. 41, No. 324, December 
1109. 


and Analysis Velocity Variation Turbulent Flow,” Simmons and 
Salter, Proceedings, Royal Soc. London, Vol. 145, 1934, 242. 
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fluctuations. the distribution function the same general form, different 
degrees turbulence can indicated some average value the quantities 
u,v,andw. The most significant average that called the root mean square, 


Velocity, in Inches per Second 


4.0 
Time, in Seconds 


1.—TRANSVERSE VELOCITY FLUCTUATION FOR TURBULENT FLOW CHANNEL 


indicates the spread the distribution curve; also, its square, such 
proportional the energy the turbulence for that velocity component. 

Various types turbulence can designated depending the relative 
magnitude and interrelationship and The simplest type turbu- 
lence, from the mathematical viewpoint, that referred isotropic turbu- 
lence. this type turbulence and also quantities the 
type ware zero. These terms are referred mean products 
statistics. quantities such are different from zero, mean shearing 
stress and mean velocity gradient exist the fluid. ordinary two- 
dimensional flow the mean shear stress, represented as, 


which unit density. simple derivation this given Messrs. 
Rouse* and Isotropic turbulence, therefore, cannot exist near 
boundary. only approached the center large conduits and 
artificially created turbulence, such that beyond screens grids. The 
turbulence hydraulic jump probably approaches isotropy. 

Another fundamental concept that useful discussing turbulence that 
the correlation coefficient. order characterize turbulent condition 
more completely, addition the quantities v?, and which characterize 
the intensity the turbulence, scale factor necessary. Two different 


Values of fv) 
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turbulent flows may have the same intensity, but the average size the eddies 
may entirely different. order identify the scale the turbulence, 
Professor Taylor® has suggested that the instantaneous velocity components, 


Equation (1) 


0.6 


Values 


Velocity Fluctuation, Inches per Second 


DISTRIBUTION CURVE FOR TRANSVERSE VELOCITY FLUCTUATION 


either measured simultaneously two different points, along either 
the z-axis the y-axis, over appreciable period time. Suppose 
measured simultaneously two points distance apart. Let the ve- 
locity one point and that the other point the same instant, and define 


The term known statistics correlation coefficient. Now, the 
two points are close that they are within the average size eddy, and will 
tend approximately equal magnitude and the same sign, and 


J 
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will nearly unity. However, the two points are average eddy diameter 
more apart, will approach zero. Experiments verify this fact. 
fessor Taylor then defined quantity, which characteristic the scale 
the turbulence, thus: 


which the distance when approaches zero. Experimental values 
R,, determined the center 12-in. 12-in. water channel where 
the mean velocity was 7.8 in. per sec, are shown Fig. Investigations 
the variation the correlation coefficient have proved helpful the study 
wind-tunnel turbulence, and should undoubtedly useful identifying 
various conditions hydraulic turbulence. 


TURBULENCE AND ENERGY CONSIDERATIONS 


First, the concept energy turbulence should clearly 
The kinetic energy per unit volume point turbulent fluid any 
instant for the velocity component along the z-axis can represented 0.5 


(The term represents the weight unit volume water.) 


The average kinetic energy then 0.5 since the arithmetic mean 


value wis zero. The term 0.5 the kinetic energy the mean velocity. 
The energy represented the other two velocity components 0.5 and 
0.5 Therefore, the total energy due turbulence per unit volume is, 


the turbulence uniform the direction flow (that is, does not 
change, such uniform pipe), the energy present due the turbulence 
need not taken into account any energy equation. However, the energy 
turbulence varies from point point, such occurs conduit expansions, 
beyond bends, through valves, etc., its omission, any study energy changes, 
results error. important item regarding turbulent energy that 
form degraded energy and will eventually transformed into heat, 
and can never recovered. 

How does turbulence cause the dissipation energy? All energy dissi- 
pated into heat through the action viscosity. there were viscosity, 
there would practically energy loss, and any turbulent energy created, 
such sudden expansion, would tend persist indefinitely. The high 
internal shear stresses set the turbulent eddies cause the large energy 
dissipation turbulent flow. 

pipe flow, for instance, potential energy being dissipated constantly 
into the form heat. For turbulent flow, simplified derivation for the rate 
transference the potential energy producing flow into energy turbulence 
follows: Let the mean tangential shearing stress point circular 
pipe The mean difference velocity between two layers fluid AU. 
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The work done per unit time the force area (Az Az) then equal 
Az). The work done per unit time per unit volume fluid 
or, the limit, which the mean velocity gradient. The 


mean unit shearing stress, point distant from the center pipe, 
radius equal to: 


which: loss pressure head pipe length, and unit fluid 
weight. This shear stress can also represented as, 


which: viscosity; unit density; and mechanical viscosity due 
turbulence (may vary from point point). Except near the wall, fully 


developed turbulent flow, the quantity very small comparison 


into turbulent energy rate represented per unit volume. The 


turbulent energy then changed into heat the action viscosity. The 


Values of Correlation Coefficient, Rz 


4 
Volues 2, Values of n, in Cycles per Second 


FOR TRANSVERSE DIFFERENT FREQUENCIES 


rate which potential energy transferred into turbulent energy not 
necessarily equal the rate energy dissipation into heat, any one point. 
The turbulent energy can diffused other points the action the eddies, 
and then dissipated into heat. 
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The mechanism which turbulent energy dissipated into heat quite 
complicated, even the simple case isotropic turbulence. The rate 
this energy dissipation proportional the viscosity times the mean square 
the various instantaneous turbulent velocity gradients—that is, quantities 


2 
such The partial derivative sign used because the component 


also varies with The general formula for turbulent energy dissipation 
given and discussed Professor For isotropic turbulence Pro- 
fessor Taylor has shown that the mean rate dissipation energy per unit 
volume fluid is: 


Various simple relations exist between the different mean turbulent velocity 


2 
gradients isotropic turbulence. For instance 


2 
One method determining the value quantity such use 


the correlation coefficient concept. Recalling the definition from 


Equations (3), the values v2)? are obtained for small values then 
approximately, 


Thus, 


Therefore, 


istic the turbulence far energy dissipation concerned. The value 


2 
for the data shown Fig. 0.56 in. The value then 3.32, 


since was 0.52. 

analyzing the velocity fluctuation single point, the statistical 
distribution function was found indicate the frequency occurrence 
similar type statistical analysis the distribution the energy the 
turbulence, such 0.5 among the various frequencies. This will useful 
showing what the average magnitude the frequency the velocity 


Theory Turbulence,” Taylor, Proceedings, Royal Soc. London, Vol. 
430. 
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function f(n) determined such that f(n) indicates 
the most probable proportion the total energy, represented 0.5 
lying between the frequencies and dn). More exactly, f(n) does not 
indicate the fractional part the total energy lying between the frequencies 
nand dn), but rather, the probability the total energy lying that 
range frequencies. The method the determination from experi- 
mental velocity fluctuation data, has been indicated the writer previous 
can shown that, 


0 0 
shown Fig. for the velocity data shown partly Fig. The value 


3 2 


checks very well with the calculation the correlation coefficient method 
given Equation (10). That the two methods are mathematically equivalent 
can demonstrated. However, different data are used each method 
obtain the same final result. 


CHARACTERISTICS TURBULENCE 


Turbulence eddies are capable causing the exchange momentum, heat, 
matter, from one layer fluid another. Near solid boundaries the 
existence mean velocity gradient across turbulent fluid indicates that 
momentum being transferred from one layer fluid another. temper- 
ature gradient indicates the transfer heat; and suspended material gradient 
turbulent stream indicates the transfer matter. The characterization 
turbulent flow its diffusing exchange power important item 
any investigation turbulence. The process mixing turbulent 
fluid stream often referred ‘‘eddy 

The laws molecular diffusion have been applied mass transfer processes 
turbulent fluid However, the complete validity this procedure 
open question, since the motion fluid masses not comparable that 
gas molecules. common and convenient method studying the diffusing 
power turbulence observe the spread heat, matter, color from 
point source the heat, matter, color injected the turbulent fluid stream. 
Experiments were made injecting constant stream color through 
fine tube, and also injecting colored droplets mixture carbon tetra- 
chloride and benzol, having the same specific gravity water, and photo- 
graphing the transverse spread the color, the droplets, with motion- 
picture camera. The important quantity determine from such experiments 


Van Driest, Jun. Am. Soc. E., Proceedings, Fifth International Cong. Applied Mechanics, 


Eddy Diffusion,” Towle, Sherwood, and Seder, loc. cit., 1938. 
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the mean square cross-stream transverse travel the color various 
distances downstream from the point injection. Fig. shows some sketches 
typical motion-picture frames the color and the droplets. The transverse 


Water Surface 


Color Stream 


AND COLORED IMMISCIBLE 


travel any distance downstream called and its mean square value 
Professor has developed theory relating and His general 
premise that the statistical correlation between the velocity fluid mass 
one instant and another later instant high when the time interval 
small, and this correlation approaches zero when the interval becomes large. 
This theory shows that varies for small values and gradually 


approaches direct variation with increases. Professor 
showed that 


which: mean square transverse velocity fluctuation; correlation 
coefficient between velocity fluid mass one instant and instant, 
t 


0 


Fig. shows experimental data obtained the writer the relation 
and turbulent water stream near the center 12-in. 12-in. 


open channel point where the mean velocity was 7.8 in. per sec and 
was 0.78 in. per sec. Note the initial linear relationship and and, 


Continuous Movements,” Taylor, Proceedings, London Math. Soc., Vol. 20, 
1921, 196. 
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larger values the linear relationship and Note also that the 


initial slope the determines the value Vv? according Equa- 
tion (14). 


> 


Values of and /V? 
w 


5 
Distance, x, Downstream, in Inches 


Data DIFFUSION TURBULENT WATER STREAM 
(Bars Over Squared Terms Indicate Mean the Squares and not the Square the Mean) 


mathematical relationship between and not very convenient 
means characterizing the diffusing power turbulence. diffusion coeffi- 
cient would the logical thing determine. Using the ideas molecular 
diffusion theory, such coefficient would defined as: 


4 


The value for the data Fig. 0.059. The dimensions 


2 ‘ 
are velocity times length. The velocity the root mean square 


transverse velocity, and the length term the length associated with 
mixing processes. This length term proportional the size the turbulent 
eddies, and Prandtl’s so-called length.” 

should added that the relative concentration the color droplets 


transverse the turbulent stream varied according the normal error curve; 
that is, 


| | | aJ/y? 0 10- Ly 
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The probability that particle would lie between and given 

knowledge the diffusing power turbulence most important 
studies suspended material transportation. The diffusion coefficient 
proportional the so-called coefficient eddy viscosity, mechanical 
viscosity, usually designated and occurring the formula for shear 
stress caused turbulent exchange momentum: 


The proportionality between and must determined experimentally, 


EXPERIMENTAL TECHNIQUES 


The different quantities which desirable determine order 
characterize turbulence have been discussed. more detailed explanation 
the experimental methods that have been used and are available for the 
determination these various quantities will outlined briefly. The in- 


practically impossible determine all three these fluctuating compo- 
nents the same instant. aeronautics and wind-tunnel studies the hot- 
wire anemometer has become the standard instrument for measuring 
cross components, and can determined with special type hot-wire 
anemometer. The hot-wire anemometer has not been used very successfully 
water, due primarily difficulties such collection air bubbles and 
other foreign matter the wire, the extremé magnification the fluctuating 
current that necessary high-velocity water, and various other troubles. 
Nevertheless, has been used water lower velocities and undoubtedly 
will find further use. 

One the most convenient methods determining simultaneous values 
pended the flowing water. illuminated, such particles make streaks 
motion-picture film, and from the length and direction the streaks two 
components the velocity can determined. Particles made from mixture 
carbon tetrachloride and benzine, having specific gravity equal that 
water, have been used. The particles are illuminated narrow plane 
light, that only those particles the light are visible. 

Another method used the writer obtain values the cross velocity, 
for very short intervals time, involves thin color stream injected into 
the stream through fine needle. motion picture was taken this stream 
close the needle. The velocity for very short intervals time obtained 
measuring the transverse travel for short distances, downstream. 
Distance short enough that the direction the color stream practically 


National Advisory Comm. for Aeronautics, Technical Report No. 581. 
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values enable the calculation Such color stream shown 
Fig. The writer has found that good measurements simultaneous 
values and can made photographing, with constant speed motion- 
picture camera, the travel the immiscible droplets for distance about 
from the point injection. magnifying type close-focusing lens 
used the camera that picture obtained area about in. 
Following the movement the particles from frame frame the film 
permits the determination and and sufficient data are obtained the 

Another method determining from diffusion experiments, 
finding the initial slope the shown Fig. and using Equation 
(14). The energy dissipative characteristics are determined finding 
factor proportional the mean square the turbulent velocity gradients 
such which was defined Equation (11), the factor defined 
Equation The term can determined measuring the correlation 
and the term calculated analyzing the velocity fluctuation 
single point. air flow determined use two hot-wire ane- 
mometers placed distance apart. The writer was able measure 
injecting color two points, distance apart, simultaneously, and photo- 
graphing the two streams together. 

The determination the diffusion coefficient diffusion experiments has 
been discussed. Undoubtedly, this the most direct means finding the 
diffusion characteristic turbulence. Evidence exists that the diffusion 


coefficient related the intensity the turbulence given and 
the length factor The term proportional the size the eddies, 


and the greater is, the greater the diffusion for given value 

Unquestionably, the photographic method accurate and ideal method 
for making fundamental studies turbulence the laboratory. time- 
consuming and tedious, particularly the analyzing the data recorded the 
motion-picture film. instrument that would give accurate recordings 
the instantaneous velocity direction either longitudinal transverse 
the flow would useful. course, such instrument would have 
compared with data obtained the photographic method order insure 
that following the velocity fluctuations accurately and not introducing 
turbulent effects its own. forward step the development such 
instrument has been made who has developed so-called 
“angularity” indicator which measures the mean angle flow. This angle, 
course, proportional the mean transverse velocity. 

The cup type current meter very sensitive turbulence. However, 
because its inertia does not follow all the velocity fluctuations exactly. 
does give some information regarding the magnitude the velocity fluctua- 


Meeting, 1938. 


ten 
) 
‘ 
as 
4 
f 
‘ 
‘ 
3 
the 


1400 THEORY TURBULENCE Papers 


tions and the size the eddies, the speed rotation the vane obtained, 
some automatic recording arrangement, for short intervals time over 
extended period. 


RELATION TURBULENCE RESEARCH PROBLEMS 


Energy Dissipation—There are many problems aeronautics, physics, 
chemical engineering, and meteorology which knowledge the turbulence 
mechanism quite important. However, this paper will confined entirely 
outlining some specific problems hydraulic engineering which could 
solved more accurately and more completely with information 

Considering the general problem energy dissipation, there are various 
individual problems which require knowledge the manner which energy 
dissipated. are some design problems which desired obtain 
the maximum energy dissipation possible, and some which the energy dissi- 
pation must minimum. 

conduit expanding cross section, the additional losses energy asso- 
ciated with such expansion are directly attributable the creation 
turbulence energy beyond that ordinarily present uniform flow. such 
expansion, whether gradual sudden, there first increase the energy 
turbulence and then decrease. 

order study the energy transformations completely, the energy 
turbulence given Equation (5) must found, addition the ordinary 
kinetic energy, computed from the average velocity distribution. Fig. 
shows some typical velocity data obtained the method streak motion 
pictures 3-in. circular pipe. Fig. shows similar data for section 
3-in. 5-in.expansion. The rate flow both cases was 0.0818 per sec. 
computing the turbulent energy the transverse velocity component, was 
assumed equal this quite true for flow circular conduit, except very 
near the wall. The total mean kinetic energy obtained from the mean 
velocity distribution given by, 


0 


The total energy turbulence 


0 


The value the 3-in. pipe 0.255 ft-lb, and 0.00793 this 
case 3.1% Em. The value the expansion the 3.75 diameter 
section 0.157 ft-lb, and 0.0135 ft-lb, 8.6% Em. Note the relatively 
higher value the turbulent energy the expanding portion. Data such 
that given Fig. enable the accurate study the energy changes taking 
place from point point. 

Although turbulent energy can created very quickly, for instance, 
sudden expansion, the dissipation that energy does not occur immediately. 
The expression loss” used many times describe the energy loss 


Octe 


| 
é 
> 
| | 
2 | 
ig 


October, 1939 THEORY TURBULENCE 1401 


0.10 


Mean Velocity, in Ft per Sec 


Mean Square Velocity Fluctuation, (Ft per Sec)? 
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sudden expansion. This misnomer and gives wrong impression 
how the energy loss really occurs. The sudden change velocity causes 
treme turbulence; the energy represented such turbulence then gradually 
dissipated downstream the action viscosity. the final analysis, vis- 
cosity, friction, the cause the dissipation energy into heat. 

desired dissipate considerable energy, such below dams 
spillways, reduce the mean velocity flow steep chutes conduits, 
much can learned regard such problems from the theory the turbu- 
lence mechanism. The rate which potential energy per unit length 
circular conduit converted into turbulent energy is: 


which: loss head per foot conduit; total weight water flow- 


ing unit time; unit shear stress distance from the center; and 
piperadius. For fully developed turbulent flow the shear stress equal 


which proportional the diffusion coefficient, more specifically, 


the quantity which the intensity the transverse turbulent 
velocity, and length factor determined the size the eddies. 
fore, the rate potential energy transformation can increased increasing 
which can done introducing large-scale roughness produce large, 
intense eddies. 

For given intensity turbulence the smaller the eddies the higher the 
rate energy dissipation will be, seen from Equation (11). The 
will smaller when the eddies are smaller, and therefore the average rate 
energy dissipation will greater. Thus, appears that once the poten- 
tial energy mean kinetic energy flow the form turbulent energy, the 
dissipation this turbulent energy can hastened the breaking down 
the large eddies into smaller ones. This principle can readily taken ad- 
vantage the destruction kinetic energy below dams and spillways. 
Thus, converting potential energy into turbulent energy, such for flow 
pipe channel, the mean velocity can reduced the creation intense 
large-scale eddies. destroy the kinetic energy high-velocity flow below 
dam, desirable produce intense, small-scale eddies. 


TRANSPORTATION SUSPENDED MATERIAL 


The relation the distribution suspended material the exchange 
characteristics turbulence has been quite well For equi- 
librium conditions the general statistical equation is: 


which: instantaneous vertical velocity component; instantaneous 


Theory Turbulent Flow and Sediment Transportation,” O’Brien, Am. 
Soc. E., Transactions, Am. Geophysical Union, 1933, 487. 


Suspended Material Hunter Rouse, Proceedings, Fifth International 
Cong. Applied 1938. 
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sediment concentration per unit volume; average velocity fall single 
sediment particles still water; and mean sediment concentration any 
point above the stream bed. The left side Equation (22) indicates the net 
transport material upward the turbulence, and the right side the quan- 
tity that settles due gravity. Under equilibrium conditions they are equal. 
The value any instant can represented Then the 
upward, will positive, and when down will negative, since 
the concentration sediment general increases from top bottom 
stream. The mean value then AN, since the mean value 
zero, because the arithmetic mean value zero. The question now 
arises how large is. This depends the distance from which the 
suspended material brought down or, other words, the scale 
the turbulence size the eddies. Introducing length factor, character- 
izing the diffusing scale the turbulence: and Equation (22) be- 


comes 


The term which may called proportional the diffusion coefficient 
defined Equation (16). Whether identical the momentum exchange 
coefficient, occurring the shear formula, must investigated experimen- 
tally. they have been assumed equal, and preliminary investigations 
indicate that they are identical. Integration Equation (23) leads rela- 
tion between the unknown sediment concentration, any point above the 
stream bed, and known concentration, point distance above the 
bottom stream canal: 


Log. 


This expression for silt distribution has been checked both artificially created 
turbulence and also natural streams. The term varies vertically 
natural stream, tending increase away from the bottom, reach maximum 
about the mid-depth, and then decrease the surface approached. 
However, the exact variation has not yet been determined. 

Another important problem regarding turbulence and suspended material 
the predication the actual concentration some point (for instance, just 
above the bottom), knowing the composition the bottom, the various hy- 
draulic factors, such mean velocity and depth flow, and the intensity and 
diffusing power the turbulence. Such relation necessary the expression 
for sediment distribution (Equation (24)) great practical usefulness. 
The mechanism the picking sediment from the stream bottom the 
turbulent eddies merits intensive experimental investigation. However, such 
investigation will little fundamental value the intensity and the 
diffusing power the turbulence not determined the same time. 

the problems silt sedimentation, the diffusion characteristics turbu- 
lence are also extremely important. illustrate how the theory turbulent 
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diffusion can applied the sedimentation silt, two general problems 
analyzed. 

Problem 1.—Consider the problem which sediment settling turby. 
lent stream, and desired know what the silt distribution various 
distances downstream from the point which settling starts, and which 
point the silt distribution known. The assumption will made that the 
turbulence such magnitude that most the silt that reaches the bottom 


Surface 


not picked again. Referring Fig. the concentration section (1) 
point above the bottom units per unit volume for particular size 
silt, what will the distribution these units some distance down- 
stream? the particles were weightless water they would diffused ac- 
cording the laws turbulent diffusion, and their distribution would ac- 
cording the normal error curve given Equation (17) (with 
Fig. the point would the height Then the average concen- 
tration point the section downstream will be: 


sion distance measured from the initial elevation, the concentrations 
and are per unit volume, can taken unity. words, then, 
the area under the distribution curve between 0.5) and 0.5) 
times the concentration per unit volume The term definite 
Since the silt settles still water velocity apparent that the 
highest concentration (the point the distribution curve where will 
found distance below or, other words, the distribution curve will 
shifted downward distance equal For simplification the velocity 


used the average velocity the entire stream. The concentration 
point will then 
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2 


topography. obtain the total average concentration per unit volume 
values for all values hi, from must summed: 

Note that function and that and are constants for each summa- 
tion any one point section (27) gives the relation between 
the initial silt distribution and the distribution another section downstream 
when settling occurring. 

Problem 2.—The other part the sedimentation problem that finding 
how much silt settles distance course, the difference between the 
total concentrations two points distance apart what has settled the 
time required for the water flow distance However, this roundabout 
process finding the quantity that has settled the length therefore an- 
other method used finding this quantity. first pre- 
sented the method outlined herein. The probable quantity silt that will 
settle distance from all the silt that passes point given time 
given the shaded part the area under the distribution curve, shown 
Fig. 


(AM), 
which: total silt passing unit area point given period 


The total silt that will settle unit time, the 


time; and 


reach then given 


\ 2 To 0 


Equation (29) can simplified, especially approximate expression 
The probability integral can evaluated readily reference 
books probability and statistics. However, for the present deemed 
sufficient present the basic equations only. The next important step 
conduct fundamental experiments order check these basic relations. 


MISCELLANEOUS PROBLEMS 


addition the general problems energy dissipation and suspended 
material transportation, which practical application made knowledge 
relating turbulence, few others may stated, follows: 


(1) Many explanations have been proposed the past the reason for 
the surface velocity open channel stream being less than the maximum. 


Velikanov, Transactions, Scientific Research Inst. Hydrotechnics, Leningrad, 1936, 50. 
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The fundamental explanation related the turbulence mechanism, and 
specifically the process diffusion fluid masses the form eddies 
which are formed the channel boundaries and then diffuse into the 
body the stream. 

(2) The relation between high surface-water velocity and the entrainment 
air dependent the degree turbulence the water, which 
depends the depth flow and the channel roughness. 

(3) problem that has recently been studied detail that the relation 
the turbulence conduit the accuracy the Allen salt-velocity method 
high degree turbulence (that is, for relatively high velocity) the method 
more accurate than for slight turbulence. Errors may also result the measure- 
ments are made reach pipe which the turbulence changing, such 
beyond bend, valve, pipe expansion, other turbulence producing device. 

(4) Instruments for measuring mean velocity and static pressure turbu- 
lent water stream are affected more less the fluctuating velocities. Cur- 
rent meters, especially those the cup type, greatly over-register turbulent 
water. The velocity head read pitot tube greater than that due the 
average velocity depending the degree turbulence. The static pressure 
read static-pressure tube placed turbulent fluid has reading greater 
than the true static pressure amount proportional the mean squares 
the cross-velocity components. However, with intensities turbulence ordi- 
narily present water conduits channels, the error the reading 
designed pitot tube static-pressure tube should less than 1%. 

(5) conclusion, should mentioned that aeronautics the measure- 
ment the intensity, and the scale turbulence, rapidly becoming stand- 
ard laboratory procedure all fundamental experimental work. the near 
future, undoubtedly, such measurements will made part many hydraulic 
laboratory investigations. For example, model studies are made prob- 
lems involving energy dissipation problems dealing with suspended 
material transportation bed-load movement, order that such studies 
fundamentally sound, dynamic similarity must maintained far intensity 
and scale the turbulence are concerned. Turbulence should thought 
terms definite parameters instead just qualitative descriptive term 
relating general condition flow. 
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PAPERS 


EFFECTIVE MOMENT INERTIA 
RIVETED PLATE GIRDER 


The importance the plate girder structural member has warranted 
many tests its behavior under load. approach the solution this 
problem, series tests was planned discover, possible, how much the 
flange stresses plate girders are affected the rivet holes the flanges. 
other words, attempt find what, this report, called the effec- 
tive moment inertia, and, far the writers know, the first attempt 
secure such information. This report describes the methods used Swarth- 
more College, Swarthmore, Pa., for testing two plate girders and gives 
resumé the results. The intent was obtain knowledge the action 
plate girder under working conditions; therefore the girders were not loaded 
until they failed. 

Both sets tests show that the effective moment inertia very close 
the gross moment inertia, and that good correlation measured flange 
stresses with stresses exists the gross moment inertia used. 
Some data are also presented justify the location the neutral axis the 
center gravity the gross cross-sectional area. Finally, the writers present 
formula for effective moment inertia. 


INTRODUCTION 


customary for structural designers consider net moment inertia 
the design plate girders. For the purpose this paper, ‘‘net moment 
will defined the moment inertia computed deducting 
the effect the material removed holes the tension flange. For clarity, 
“gross moment inertia” defined the moment inertia the solid cross 
section the plate girder, allowing deductions for the effects holes 


Prof., Civ. Eng., and Chairman, Div. Eng., Swarthmore Coll., Swarthmore, Pa. 
Asst. Prof., Civ. Eng., Swarthmore Coll., Swarthmore, Pa. 
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the tension flange. The use moment inertia” connection with 
the beam theory implies displacement the neutral axis that can com. 
puted. The net moment inertia assumes that holes the gage 
the tension flange are slotted for the full length the girder. 
say that the usual plate girder, having rivet pitch the flanges from 2.5 
in. in., does not represent girder with slotted holes the full length the 
tension flange. The practice designers, who use the moment 
when computing the deflection plate girder, has been confirmed 
some experimental data. 

secure answer the question whether the net moment inertia 
the gross moment inertia, intermediate value which may called 
the effective moment inertia, should used the determination the flange 
stresses plate girder, the same section the same plate girder was tested 
under pure moment with the flange the following condition: 


Series solid flanges (no holes outstanding legs flange angles 
cover plates) 

Series open holes outstanding legs flange angles and cover plates; 

Series with machine bolts the holes described series and, 

Series with rivets the holes described series 


evident that the same section the same plate girder tested under 


the conditions described series and any displacement the 
neutral axis will disclosed the strains the web and the flanges. 


GIRDERS 
The choice test girders was predicated upon maintaining depth-to-span 
ratio consistent with practice, the same time selecting size that could 
handled readily the laboratory. The test girder chosen has depth-to-span 


| 
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SHEAR SECTION PURE MOMENT SHEAR 


ratio 10.6. Fig. detail drawing the test plate girder, PT-1, and 
represents the girder the state which was received the laboratory. 
The holes marked were drilled the shop after tests Al, A2, and had 


4 

us 

4 

3" 

3"" 25 


October, 1939 MOMENT INERTIA GIRDER 1409 


been completed. Test girder, PT-2, same when received, except that 
the rivets connecting the web plate flange angles the pure moment section 
were spaced 2.5 in. apart. holes were punched the outstanding legs 
the flange angles the cover plates the center section the girders. 
The makeup the girders was the same throughout the full length, namely: 


Web plate in. in. 

Flange angles 2.5 in. 2.5 in. in. 


The depth, back back angles, was 12.5 in.; and, 8-in. cover plates were 
used that strains could measured the side these plates next 
the 6-in. cover plate. This procedure gave useful information and was 
abandoned. 


‘ 


2—View Test Set-Up, THE FLANGE PLACE 


The “pure moment section” the test beam was the only one which the 
four comparative conditions girder section could secured. The cover 
plates were developed under the conditions series and bolts the 
shear section only. Fig. illustrates the typical installation the flange 
clamps necessary prevent the flange plates the compression side from 
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buckling during series and The clamps were designed that two 
set-screws came into point contact with the flange plates and two 
screws came into contact with the inner side the outstanding legs the 
flange angles along the gage lines. The plates were pulled into contact and 
held there Their construction was such that they added nothing 
the flange area. 


The method applying the loads was follows: One end the test girder 
rested upon roller upon the platform 100,000-lb testing machine; the 
other upon fixed cylindrical reaction point, supported steel trestle and 
adjusted the same height the platform the ‘testing machine. The 
load was applied the test girder equalizer beam (see Fig. through 


East End West End 


5) | 


2P 
Equalizer Beam 


Loading Beam 


Movable Head of 
Testing Machine 


Tie Rod 


Reaction on 
Platform of 
Testing Machine 


Test Section 


5'0" 


ANALYSIS FOR PLATE-GIRDER TEST 


rollers situated points 33.5 in. from the supports each end the test 
girder. (Test loads are noted all cases, this paper, and not 2P. 
Therefore the bending moment the pure moment section always 33.5 
in-lb.) The load was brought the middle the equalizer beam 
loading beam. One end this loading beam was connected the movable 
head the testing machine clevis and pin; the other end was attached 
the corresponding end the plate girder two bolts, 1.75 in. diameter, 
passing through steel slabs, top and bottom, contact being made through 
cylindrical reaction points. Contact between the loading beam and the 
equalizer beam was made means roller. this arrangement the 
concentrated load applied each load point exactly that which weighed 
the testing machine, making possible total load twice the rated capacity 
the testing machine. the presentation data the testing machine load 
the recorded value. all cases, the total load the girder twice this 
value. The east end the girder (see Fig. rested the testing machine. 

Pure Moment Section—The method measuring the deflections the pure 
moment section was follows: Two pins (0.25 in. diameter), spaced 25} 
in. each side the mid-point the girder, were passed through the web plate 
such manner that the tops the pins were level with the center line the 
web plate. These pins extended through the web plate sufficiently far 
each side that they could support steel side bars in. in. cross 
section. The holes the ends these bars were in. diameter, and the 
distance between the supporting pins was the middle the plate 
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girder third pin was placed through the web similar manner, and each 


end the pin was mounted Ames dial, reading in., with the contact 


10,000 
points resting upon the side bars. (These dials are identified the north and 
south dials.) Under this condition loading, the middle the test girder 
was subjected uniform moment, the moment due the weight the 
girder itself neglected. evident that the reading the dials gave 
the deflection the center line the web below two points that line which 
were spaced 50.25 in. apart (see Fig. 1). 

measure the stresses the girder itself, gage points with 10-in. centers 
were drilled and readings taken with 10-in. strain gage having least reading 

least reading was equivalent stress persqin. the test procedure, 
readings were checked independent observers until they agreed within six 
divisions the dial, insuring readings with maximum deviation 300 
per in. 


Assuming modulus elasticity 30,000,000 per in., the 


Fig. shows the cross sections the plate girder for each series tests. 
each case initial load was applied and then the load was increased until 


TEST SECTIONS A-1, B-1, TEST SECTIONS B-2, TEST SECTIONS A-3, B-3, D-3 


the increment deflection the pure moment section was 0.0300 in., 
slightly greater. After each series had been run, with loads give the afore- 
mentioned deflection, the load was increased 63,000 for series C-3 and D-3 
for specimen PT-1, and readings were taken the pure moment section. 
detailed statement the test program follows: 


Test Section Description (see Fig. 

Web plate with two pairs flange angles 
riveted on, but with holes the out- 
standing leg. 

A-2 6-in. cover plate added top and bottom. 
There were holes the plates the 
pure moment section. The plate was held 
contact with the flange angles the 
clamps shown Fig. 
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Test Section Description (see Fig. 

A-3 8-in. cover plate added the 6-in. cover 
plate, top and bottom. this case, also, 
the pure moment section had holes 
and the plates were held contact the 


clamps. 

B-1 Web plate and flange angles were tested 
test section 

B-2 6-in. cover plates added, clamped, and 
tested test section 

B-3 8-in. cover plate added the assembly 


described test section B-2, the two plates 
clamped and tested test section A-3. 

C-1 machine bolts placed the open holes 
in. diameter) the outstanding legs, 
tightened with proper washers, and tested 
test specimen B-1. 


C-2 6-in. cover plates added, bolted on, and 
tested test specimen B-2. 
C-3 8-in. cover plate added the assembly 


described under test specimen C-2 and 
tested test specimen B-3. 
D-3 The two cover plates were riveted the 
flange angles and tested test speci- 
men C-3. 


series (the section with open holes), the plate girder with the 6-in. 
and 8-in. plates attached was taken drill press and holes were drilled 
through the outstanding leg the flange angles and the cover plates, both 
compression and tension flanges. specimen PT-1 the pitch was in.; and 
specimen PT-2 the pitch was 2.5 in. series the holes were open 
the pure moment section. (The pitch the rivets specimen PT-2 was first 
made in. specimen PT-1, and tests were run which checked specimen 
PT-1 closely that they are not reported herein. Then the girder was re- 
turned the drill press and the additional holes drilled make the pitch 


Beginning with the assumptions the common theory flexure, can 
proved that the neutral axis must pass through the center gravity the 
cross section beam under flexure. This substantially correct for homo- 
geneous beams the usual kind. However, the design riveted plate 
girders, the engineer must decide whether will locate the neutral axis 


the center gravity the gross cross section, the center gravity the 


net cross section, considering that the material that filled the holes has been 
removed from the tension flange, and that the holes the compression flange 
are completely filled. the net cross section existed each section, 
would were obtained slotted holes running the full length the girder, 
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and the compression flange were solid, the center gravity should un- 
doubtedly computed from the net cross section. The net section not 
continuous, however, and questionable whether the rivets fill the holes 
the compression flange sufficiently consider the compression flange solid. 

That the neutral axis remains near the center line the web plate shown 
Table where data from tests the two plate girders, PT-1 and PT-2, 


TABLE OBSERVED AND COMPUTED DISTANCES THE 
Axis ABOVE THE CENTER LINE, INCHES; 


(a) PLATE (b) 6-In. 
AND FLANGE Cover PLATES 
ADDED 


Cover ADDED 
Description 


Observed, Web Data: 
Specimen 

Observed, Flange Data: 

Specimen PT-1............ 

Computed* 


worm 


Computed the basis net section. 


are presented. observed results were obtained measuring strains 
10-in. gage lines both sides the web plate. The gage lines were dis- 
tances in., in., and in. above and below the center line the girder. 
Astraight line was passed through the plotted data and the deviations obtained 
directly. The deviation the neutral axis was also computed using the 
measured strains the outside fibers the girder, and these results are com- 
pared with those obtained from the web data. The experimental results all 
indicate that the neutral axis deviates only slightly from the center line the 
web plate, and not has been considered deviate the designers using 
the net section. This means that the deformation the compression flange 
almost great that the tension flange. 


Tensile tests were made the National Bureau Standards eight test 
coupons from specimen PT-1, obtained from the extra pieces material sent 
the fabricator (see Table 2). The coupons had reduced width 0.5 in. 
and gage length in. The coupons for specimen PT-2 were tested 
Swarthmore College. The results were similar that the data are omitted 
for present purposes. The average value the modulus elasticity was 
30,000,000 per in. 


OBSERVED EFFECTIVE INERTIA 


the pure moment section, the deflection length 50.25 in. was read 
the nearest 0.0001 in. Ames dials. observe values the effective 
moments inertia, the moments inertia the gross section, (Column (2), 


A-1 B-1 C-1 A-2 B-2 C-2 A-3 B-3 C-3 D-3 
0.19 0.1 0.14 0.16 0.08 0.07 0.09 0.04 0.05 0.10 
0.09 0.0 0.01 0.04 0.09 0.05 0.17 0.07 0.07 0.10 
0.11 0.1 0.07 0.11 0.05 0.11 0.09 0.05 0.03 0.14 
0.34 0.63 0.25 0.68 0.68 0.18 0.66 0.66 0.66 
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Table 3), were assumed correct for series and load, was applied 
until the deflection the pure moment section was 0.0300 in. series 
and loads were applied each test section the correct value cause 


TABLE ror TENSILE STRENGTH; SPECIMEN PT-] 


Thickness, 


ulus tional Tensile in. 


elasticity limit strength 


Coupon 


Coupon cut from: 


Web plate 

Web plate 

6-in. cover plate 

6-in. cover plate 

8-in. cover plate 

8-in. cover plate 

Flange angles 

Flange angles 27,900 


deflection 0.0300 in. Then was assumed that since the deflection, 


the plate girder was the same each these tests, and since 


which and are the load and moment inertia for any particular section. 
For example, the load that caused deflection 0.0300 in. for section A-2 was 
26,600 and the gross moment inertia was 315 the load that caused 
the same deflection for section B-2 was 25,200 lb. Therefore, the observed 
25,200 
this manner, values the observed effective moment inertia for each 
test section were computed and tabulated Column (2), Table com- 
paring the values Table 3(a) and Table should recalled that the 
pitch this pure moment section was in. for specimen PT-1 and 2.5 in. for 
specimen PT-2. Instead tabulating Column (2), Table the value 
the gross moment inertia opposite the test section series the value 
was computed from the value the deflection, 0.0300 in., determined 
the moment-area method, using the value determined from the test 
coupons. These values check the values the gross moment inertia within 
the limits the experimental error and justify the foregoing procedure. 
When the values the observed effective moment inertia are plotted, 
the graphs show that they approach the values the gross moment 
the tests proceed from that giving the greatest deviation series 


effective moment inertia for section B-2 


D-3 43.0 
| ; 
(1) 
and the observed effective moment inertia each case 
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holes), through series (bolted sections), series (riveted sections). The 
deviation very small value series D—namely, 2.1% the gross moment 
inertia for specimen PT-1 and for specimen PT-2. 


Moments INERTIA, 


Computed Observed 


Computed, Using the Method Indicated 


Test 


Net section Net Section 
section Ob 4 
Ob- serve 
served section Neutral Neutral 
section axis axis 


eccentric center 


(6) (13) 


(14) 


A-2 


B-2 298.5 315.0 270.7 19.3 19.7 19.2 19.5 18.4 18.4 19.2 23.7 20.8 20.8 
C-2 306.5 315.0 270.7 18.9 19.6 18.6 19.2 18.4 18.4 19.2 23.7 20.8 20.8 


A-3 494.0 490.8 476.9 477.4| 18.5 19.0 18.7 19.3 19.1 19.1 19.3 20.3 19.7 19.7 
B-3 460.5 490.8 430.0 19.9 20.2 20.2 20.5 19.1 19.1 20.1 24.4 21.5 21.5 
C-3 472.0 490.8 430.0 19.7 19.9 19.8 20.0 19.1 19.1 20.1 24.4 21.5 21.5 


PT-2 


C-1 179.0 192.9 19.0 19.2 18.5 18.6 17.2 17.2 18.2 22.3 19.5 19.5 
315.0 315.0 300.9 17.9 18.6 18.1 18.7 18.4 18.4 18.5 20.0 19.2 19.2 
B-2 282.0 315.0 270.7 20.0 20.9 20.2 21.0 18.4 18.4 19.2 23.7 20.8 20.8 
295.0 315.0 2/0.7 275.5| 19.4 20.3 20.0 19.2 20.8 20.8 


Neutral axis eccentric. Neutral axis center. 


Top flange. Bottom flange. 


FLANGE STRESS 


all the test sections, flange stresses were obtained experimentally 
measuring strains the top and bottom flange. 10-in. strain gage was used, 
with the gage lines placed symmetrically with regard the longitudinal center 
line the girder. Four gage lines were used for test sections A-1, B-1, and 
C-1, six gage lines for test sections A-2, B-2, and C-2, and eight gage lines for 
test sections A-3, B-3, C-3, and D-3. The layout the gage lines the 6-in. 
and 8-in. cover plates was such allow two gage lines bisect the open 
holes the flanges and later the bolts and rivets. The observed flange stress 
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was computed from average the readings all gage lines every case, 
using the observed value modulus elasticity. The stresses were computed 
for single load for each series, follows: 


Test sections Load, kips 
15.9 
B-3, C-3, and 41.6 


These particular values were selected because they caused increment 
deflection 0.0300 in. test sections A-1, A-2, and A-3 for specimen 

Using these loads the flange stress has been computed using the gross mo- 
ment inertia; using the net moment inertia with the neutral axis the 
center the web plate; and finally using the net moment inertia with the 
neutral axis the center gravity, considering the compression flange 
solid and the tension flange slotted. 

The stresses computed using the gross moment inertia are compared 
with the observed stresses Table 3(a), Columns (5), (6), (9), and (10). The 
computed flange stresses are somewhat lower than the observed stresses for 
all sections except for the initial solid flanges series For those sections 
good correlation exists. The stress the compression flange greatest for 
open holes, and gradually approaches the computed value the method 
holding the flange together changed from clamps, bolts, and rivets. 
such decrease shown for tension flange. This probably indicates that 
bolts and rivets tend relieve the stress the compression flange. 

Table 3(b), Columns (5), (6), (9), and (10), the observed flange stresses 
are compared with the flange stresses computed using the net moment inertia, 
with shift allowed the location the neutral axis from the center line 
the web plate. The computed flange stresses are about 2,000 per in. 
higher than the observed stresses throughout this comparison. 

comparison between the observed stresses and the flange stresses, com- 
puted using the designer’s net moment inertia (see and 
assuming the neutral axis above the center line the girder, accordance 
with some design procedures, shown Table 3(b), Columns (5), (6), (7), 
and agreement between the compression flange stresses very good. 
However, the computed stress for the tension flange approximately 4,500 
per in. higher than the observed. This would seem indicate that this 
design procedure would give much larger flange areas than are needed 
maintain safe stresses the tension flange. 

addition the foregoing information, Table gives the flange stresses 
computed the observed effective moment inertia, using the observed 
deviation the neutral axis given Table from observed flange data. 
The excellent agreement between this column and the observed flange stress 
independent check the observed results. 


After completing the tests where the loads were held values giving 
flange stress about kips per in., test sections C-3 and D-3 specimen 
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were subjected loads causing flange stresses kips persqin. The 
was increased and decreased increments, and strain readings were 
taken the flange both the loading and the down loading. 

The strain readings gave indication permanent set. breaking 
the scale around the rivet heads through the upstanding legs the flange 
angles was noticed. close observation showed perfect set orthogonal 
breaking the scale was noticed around the rivets through 
the flange plates, although this region was stressed higher. re- 
membered that the holes the cover plates were drilled, whereas the holes 
through the web and flange angles represented shop practice punched holes. 
investigation one these holes showed poor alinement the plates and 
rough perimeter. possible that the condition the holes accounts 
for the observed strain lines? 


CONCENTRATIONS 


Stress concentrations, with their accompanying areas plastic flow around 
the holes the flanges plate girder, present some questions regard 
the validity and the interpretation the strain gage readings explained herein. 
treated this paper strains have been reduced terms stress applying 
Hooke’s law. There may objections this procedure two grounds: 
(1) That this does not represent the 
stress any point; (2) that the aver- 
age stress through the net section rep- 
resents the true stress. answer 
these objections, the writers would 
point Fig. showing the way the 


stress varies across the net section. Net Section 
known that the maximum stress 
reaches the elastic limit, and further 


there nothing that can done 
about it; second, unit stress arrived 
using net area has meaning. 

reasonable assume that 
any strain gage, operating 10-in. 
gage length over the holes, will reg- Uniform Tension 
ister some plastic, well elastic, 
deformation. The writers not 
know how separate these compo- 
nents. Nadai states that the 
areas plastic flow the region 
hole are the shape cross, 
the direction the cross making angle approximately 45° with the 
principal stresses tension also well known that the 
intensity stress highest the boundaries the are 
diameter perpendicular the applied principal stresses. Fig. depicts this 


Nadai, McGraw-Hill Book Co., Inc., New York ,N. Y., 1931, 103. 


Fic. CONCENTRATIONS 


re A 
ns 
or 
m- 
1), 
11s 4 
we 


1418 MOMENT INERTIA GIRDER Papers 


condition graphically. the principal stress increased, the stress the 
boundary the hole increases the yield point value and plastic flow 
the areas defined the cross. 

For any reasonable allowable design stress, the yield point stress will 
reached the hole, and plastic flow will occur. The designer’s choice 
kips per in. instead kips per in. will not prevent this from happening, 
although may limit, some extent, the area over which plastic flow may 
extend. 

Once ductile material strained beyond the yield point and the load 
released, residual stresses are set which will beneficial resisting the 
next applied load. 


EFFECTIVE MOMENT INERTIA 


seems reasonable believe that both the strength and the stiffness 
plate girder are affected the size and spacing the rivet holes the flanges 
plate girder. the pitch were small enough there would slot the full 
length the girder, and the diameter the rivet holes were large enough 
there would material carry the load. The writers have developed 
the following formula which incorporates the results the tests described 
this paper, and which contains terms that include the effect the diameter 
and pitch the rivet holes. 


DEFLECTION CURVE 


DIAGRAM 


diameter rivet holes; the part the rivet hole that acting 
full net section; the deflection point with reference tangent 
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the elastic curve point gross moment inertia; clear 
moment inertia (holes are deducted from the top and bottom flanges); and 
effective moment inertia. The new term, the moment 
has been introduced since these tests show that the neutral axis remains 
the middle the web plate. This fact based the practically equal 
stresses the top and bottom flanges. moment areas: 


2 
general, times constant; and, the present case: 


2 


The quantity Equation (5) must determined experimentally. The 
value probably depends upon the diameter the hole, and the dis- 
tance between rivet lines. The writers not have enough experimental 
results available determine this relation. thought that the value 


must lie between 0.6and 1.0. called Table the 


(a) PT-1 (b) PT-2 


Percent- Percent- periment 
age devi- age devi- 
ation* 


Percent- Percent- 
age devi- age devi- 


0.962 0.96: 0.976 


All deviations are computed from the experimental value. 


value obtained from experimental data compared with the computed 
“m,” using equal and 0.6, respectively. will noted that the per- 
centages deviation the computed value from the observed value change 


ul 
That is, 
Value 
auion* ation 
(1) (2) (3) (4) (5) (1) (2) (3) (4) 
C-1 +0.2 0.928 0.907 0.943 +1.6 
C-2 +0.2 0.936 0.914 —2.3 0.954 +1.9 
C-3 +1.5 0.916 0.927 +1.2 0.956 +4.4 
C-4 0.930 0.927 0.956 +2.8 
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from positive negative with these values that evidence that has 
value between 0.6 and 1.0. 


CONCLUSIONS 


measured flange stresses are smaller than the computed 
using moment inertia based upon deduction rivet holes from the 
tension flange. the case the girder with 5-in. pitch, the difference was 
considerable; with the girder having 2.5-in. pitch there was closer correla- 
tion. This was anticipated and furnishes some evidence that effective 
moment inertia should depend upon the rivet pitch well the rivet 
diameter. The net moment inertia, used design with top and bottom 
flanges alike, should yield conservative stresses. 

2.—The experimental values the effective moments inertia are less 
than the computed gross moment inertia. The effect rivet pitch shown 
clearly. However, the test girder having the 5-in. rivet pitch and riveted 
flanges gave nearly exact agreement between the experimental and computed 
values the moment inertia. 

3.—The tests indicate that the neutral axis should assumed the center 
gravity the gross cross section. 

effective moment inertia can computed using the theo- 
retical expression, Equation (5), which involves the rivet pitch. Further ex- 
perimental evidence needed establish this conclusion, and the value 

5.—Since stress concentrations reaching the yield point must exist soon 
the unit stress the gross section approaches the stress ordinarily specified 
for structural work, and nothing can done prevent them, the writers 
believe that design methods should made simple possible. Any 
method that may adopted must give sections that have been proved 
practice satisfactory. 


believed that the data from these tests point the acceptability 
the gross moment inertia for design purposes. For ordinary girders, under 
usual loadings, requires modification; for cantilever girders, girders 


small pitch, should used with coefficient such indicated 


Equation (5). 

hoped that this paper will stimulate discussion this subject and 
reveal additional data which may confirm disprove these findings. Many 
engineers may have opportunity measure strains and deflections full- 
size girders. this could done widely, and such data carefully studied, 
all uncertainties would soon removed. 


ACKNOWLEDGMENTS 


The investigators wish acknowledge the cooperation The American 
Institute Steel Construction conducting this research project. 
Frankland, Am. Soc. E., and all members the Research Committee 
have contributed many suggestions that have been incorporated this study. 
The writers are particularly indebted Wilson, Am. E., 


: 
0 
his 
git t 
Nee 
' 
‘ ~ 
A 
. 
4 
| 


October, 1939 MOMENT INERTIA GIRDER 1421 


engineer bridges and buildings, the Pennsylvania Railroad Company, for 
his sympathetic and timely discussions during the preparation the paper. 
Mr. Windenburg, physicist the Bureau Construction and Repair 
the Navy Yard Washington, C., has cooperated throughout the 
project. His assistance making several tension tests riveted plate speci- 
mensis acknowledged. Robert DeMoyer, Assoc. Am. Soc. E., graduate 
student Swarthmore and associate professor civil engineering 
Pennsylvania Military College, Chester, was great assistance while 
performing the tests. Mr. Paul Kutz and Mr. Bourdelais, the Swarth- 
more Engineering Division, are responsible for many the details developed 
during the set-up the testing equipment. The Belmont Iron Works, Phila- 
delphia, Pa., performed the final riveting the test girders free charge. 


§ 
- 

‘ 

)- 

dey 


é 
re’ 
=> 
N 
z ta 
‘ 
is 
t 
. a 
‘a - 
P 
Z 
7 


AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


DISCUSSIONS 


ANALYSIS RUN-OFF CHARACTERISTICS 


Discussion 


Orro Assoc. Am. Soc. (by discus- 
sions have indicated gratifying interest this subject, and have directed 
attention the need for further investigation, particularly infiltration and 
seepage. 

Mr. Cochrane attempts derive the shape the falling limb the hydro- 
graph inversion the rising limb. The writer believes that the dis- 
charge, any time after the time concentration and after the cessation 
the rain, must determined function outflow from storage. This 
sufficiently well established numerous studies depletion and 
the study referred Mr. Barnes. The sharp break the top the 
basic hydrograph, questioned Mr. Cochrane, and also Messrs. Zoch 
and Turner, definitely verified the hydrograph Fig. presented 
Mr. Sherman. The writer has criticism offer formulas for maximum 
flood discharge; however, these are outside the scope this paper. 

Mr. Barnes makes timely comment the unreliability the per cent 
run-off, run-off factor. This factor used where nothing better available. 
mentioned Messrs. Sherman and Zoch, better results will follow from 
the deduction infiltration from rainfall, and the consequent use 
However, infiltration and seepage are separate subject, and con- 
sideration them would complicate the subject unduly without affecting the 
characteristics being analyzed. any case, distribution precipitation 
intensity giving constant “‘net not impossible and would fit the 
assumptions which the analysis was based. 

The “subsurface storm base mentioned Mr. 
Barnes item run-off sufficient importance merit considerable 


paper Otto Meyer, Assoc. Am. Soc. E., was published November, 1938, 
Proceedings. this paper has appeared Proceedings, follows: February, 1939, Messrs. 
Vietor Cochrane, and Bertram Barnes; March, 1939, Messrs. LeRoy Sherman, Richmond 
Zoch, and Merrill Bernard; April, 1939, Messrs. Franklin Snyder, and Hoyt; May, 1939, 
Clark, Jun. Am. Soc. E.; and June, 1939, Messrs. Jarvis, and Howard Turner. 

Asst. Engr., Engr. Dept., Engr. Office, Sacramento, 

Received the Secretary August 15, 1939. 

Between Rainfall and Run-Off from Small Urban Areas,” Transactions, Am. Soc. E., 
Vol. 101 (1936), Equation (2), 151. 
Monthly Weather Review, September, 1934. 
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study. can doubt that this type run-off occurs and often 
considerable magnitude. The writer visualizes the source this run-off 
sheet water separated from the true possibly cushion 
air filling the voids the soil. Presumably, the same methods 
used for true surface run-off may applied study subsurface storm 
flow; however, must treated separately from surface run-off. 

Mr. Sherman has apparently misinterpreted the hydrograph Fig. 
termination indicated; the tail intended infinite length, only 
portion being shown. Mr. Sherman states, the use unit hydrographs 
combination, shown his example Table and Fig. 10, practical 
working tool, and the writer also uses this method, where applicable, account 
its simplicity. However, the unit hydrograph used that for storm 
duration close to, equal to, the time concentration, inconsistencies may 
expected develop. For example, with unit hydrograph for one-day 
storm, the tail several days long, and the time concentration one day, 
one-day rain produces peak equal the rainfall rate, but two-day rain 
must necessarily cause the computed peak greater the amount the 
rate discharge the tail the unit hydrograph one day after the peak 
the unit hydrograph. This cannot occur nature. The observed hydro- 
graph Fig. obviously has the characteristics described the writer’s 
analysis, the time concentration being about min. 

Mr. Sherman states that the time concentration and the distribution 
for three-day given Appendix III, should derived from the 
distribution for one-day storm. This could done, or, the example, 
the distribution for one-day storm could derived from that the three- 
day storm, which this case the basic hydrograph, the time concentration 
being three days. The process reversible. 

Mr. Zoch’s distinction between run-off and rate run-off, and between 
run-off and discharge, commended. The writer does not agree with Mr. 
Zoch’s statements infiltration and internal drainage; however, that con- 
sidered outside the scope this paper. 

The writer cannot agree with Mr. Zoch the beginning the depletion 
curve. true that with storms short duration the part the falling 
limb prior the time concentration contains both items accretion and 
items depletion; however, the succeeding portion must pure depletion. 
This has been established empirically, and may reasoned logically. The 
contribution flow from every part the basin has reached its maximum 
when the time concentration reached, and thus can only diminish. For 
the same reason, the maximum discharge can never later than the time 
concentration. The transition curve mentioned Mr. Zoch must lie entirely 
between the end the rain and the time The 
not part the basic hydrograph, occurs only the duration 
the rain exceeds the time concentration. The writer, therefore, adheres 
the statement that the basic hydrograph consists the two parts: (a) The 
concentration curve and the storage, depletion, curve. 
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Mr. Zoch states that the writer has not mentioned the determination 
the time maximum discharge. This time indicated the peaks Fig. 
and adequately described the accompanying text. 

Mr. Bernard calls attention the lack sufficient and suitable basic 
data for hydrologic studies. this very lack which forces hydrologists 
use unsatisfactory expedients such the run-off factor. 

The writer would like call attention the fact that the definition 
time concentration appears the second paragraph the section 
“Analysis the Hydrograph.” The mention time concentration the 
“Synopsis” not intended for definition. 

describing, graphically, the sequence occurrences water-shed 
during storm, Mr. Bernard makes use method which could well 
adopted generally investigations physical phenomena. Too often there 
are only the who into transports pure mathematics, and 
the who generalize from meager data. only visualizing 
actual physical occurrences that engineers can sure that their derivations 
are the right 

Mr. Snyder’s conclusion that the time concentration not the duration 
unit hydrograph correct. The synthetic unit graph Fig. includes 
subsurface storm flow (secondary base flow), and not all true surface run-off; 


the basic hydrograph true surface run-off would reach maximum 245 


per sec. separate hydrograph subsurface storm flow must con- 
structed the writer’s methods are used. The relative time peak the 
various hydrographs Fig. commented upon Mr. Snyder, are subject 
errors drafting. rain would certainly cause slightly earlier 
peak than 7-hr rain. 

The writer concurs thoroughly with Mr. Hoyt his remarks the neces- 
sity for research the behavior melting snow. 

Mr. Clark ties flood routing and hydrographic analysis together 
admirable fashion. complete hydrologic study large basin must always 
involve the construction hydrographs head-waters and tributaries and 
routing the combination. integration methods used these two steps 
real value. 

The basic hydrograph used construct one hydrograph Fig. was 
synthesized from the unit hydrograph used for the other; the differences may 
best understood from careful study the method construction. 

Mr. Clark’s question the apparent lack effect storage the rising 
limb best answered thus: The steepening the slope the water surface 
during rising stages opposes and nullifies the storage effect, whereas the flatten- 
ing the slope during falling stages increases the storage effect and thus makes 
storage the dominant factor. 

cautioning against indiscriminate use the expression for time travel, 
Mr. Clark has emphasized important point. The variation velocity 
along stream has not been given the attention deserves. assumptions 
this subject should accepted without most critical analysis. 

The effect land treatment described Mr. Jarvis must considered 
hydrologic study. The writer has encountered this problem attempt- 


4 
4 
5 5 xy 
y 
> 
1 4 
of, 
‘ 


1426 MEYER ANALYSIS RUN-OFF CHARACTERISTICS 


ing estimate run-off from irrigated lands, where the irrigation works 
almost complete retention. 

Mr. Turner calls attention the difference shape the writer’s hydro- 
graphs for storms the same volume but different intensity, and states that 
actual hydrographs are more similar. shows, Fig. 17, hydrographs 
resulting from storms and the entire time concentration, 
Actual storms lasting through the time concentration are rare large 
streams, and hydrographs resulting from storms shorter duration than 
one-half the time concentration are much more similar than those shown. 
The writer the opinion that the time concentration the Connecticut 
River Sunderland, Mass., about days. Thus the hydrographs 
Fig. are not more similar than might expected, especially the shapes 
were probably determined largely the precipitation occurring the one 
two maximum days the storms. 

adverse criticism the Committee Floods the Boston Society 
Civil Engineers was intended, the writer recognizes well the value its 
work; was indicate the necessity for refinement that the shortcomings 
the original formulas were mentioned. 

conclusion, the writer wishes emphasize that the physical sig- 
nificance the shape the hydrograph, and not the empirical mathematical 
equations describing it, that important. 

Corrections for Transactions: Delete Footnote 17b (Proceedings, March, 
1939, 514) and the proposed change which refers; Proceedings, 
November, 1938, 1773, the sentence beginning Line should read: “It 
may also used synthesize the basic hydrograph from one produced 
Line above Equation (17), delete and change Equa- 
tion (17) read See also corrections published Proceedings, 
February, 1939, 377; March, 1939, 526; and June, 1939, 1052. 
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DISCUSSIONS 


TESTS BUILT-UP COLUMNS STRUCTURAL 
ALUMINUM ALLOYS 


Discussion 


wishes thank Professor Niles for his discussion this paper. stated, 
theory for columns that fail twisting, and reliable equations for computing 
the strength such columns, are not readily available. When satisfactory 
methods become available for predicting the strength members that fail this 
manner, the use such members special locations may become more common. 

The prediction twisting failure for four the specimens Table was 
based the following reasoning: The outstanding plates would fail local 
buckling because the strength for this mode failure less than the strength 
failure were occur lateral deflection the specimen unit (compare 
the computed strength obtained Equation (6) with the strengths obtained 
Equations (4) and (5)). the supported edge the plate were free 
rotate (but not deflect), the free edge the plate would buckle into single 
wave; but the supported edge were fixed against rotation, the free edge would 
deflect into number waves, the number depending the relative propor- 
tions the plate. Since the torsional rigidity these specimens small, 
the degree fixation the supported edge the outstanding plate would 
not cause the development more than one wave the free edge. Thus, 
the center the specimen would merely twist relative the ends. 

Failure under this type action may may not result the collapse 
the specimen but occurs that load which the lateral deflection the 
free edge the outstanding parts increases rapidly increasing rate (see 
Figs. 2(m) and 2(n)). The maximum load sustained the specimen may 
greatly excess this buckling load and, course, will accompanied 
collapse the specimen. There seems method available for predicting 
the spread between the buckling load and the maximum load, but obviously 
depends the relative proportions the cross-section. Considering failure 
this light, the strengths the specimens that failed twisting were pre- 
dicted within reasonable limits Equation (6). 


paper Marshall Holt, Assoc. Am. Soc. E., was published November, 1938, 


Proceedings. Discussion this paper has appeared Proceedings, follows: March, 1939, Alfred 
Niles, Assoc. Am. Soc. 


Research Engr., Aluminum Research Laboratories, New Kensington, Pa. 
Received the Secretary September 13, 1939. 
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DISCUSSIONS 


TRAFFIC PROBLEMS METROPOLITAN AREAS 


Discussion 


Buckley, Mickle, Mitchell, and Coe amplified the original paper very 
effective manner. Although little difference opinion has been expressed, 
some questions have been raised which will treated herein. 

Mr. Herrold quite right his statement that accident-prone driver 
presents problem with which the careful driver cannot There can 
little doubt that the person who, through improper driving habits un- 
willingness conform with the requirements traffic safety, becomes re- 
peatedly involved accidents merace others who want drive walk 
safely. However, accident-proneness, indicated the repetition acci- 
dents, does not appear characteristic very large percentage drivers. 
fact, Forbes, the Bureau for Street Traffic Research Yale Uni- 
versity, concerning study 500 drivers from the State Con- 
necticut for the years that are thus left with the conclusion 
that not more than 1.3 per cent the drivers were accident-prone and that 
only 3.7 per cent the accidents the states studied could attributed 
the accident-prone 

the other hand, many street locations are definitely accident-prone 
indicated the repetition accidents almost the same kind. Analyses 
many hundreds collision diagrams street intersections, and other 
problem locations, have shown definite patterns accident experience 
individual locations. many cases these patterns are not all similar 
the traffic-flow pattern, for they show vehicles from certain directions pre- 
dominating the accident experience and those from other directions having 
little part it. Surely, one cannot conclude that the careless drivers 
come from only one two directions and that the drivers approaching from 


paper Reeder, Esq., was published December, 1938, Proceedings. Discus- 
sion this paper has appeared Proceedings, follows: March, 1939, Messrs. George Herrold, and 
Thomas Buckley, and April, 1939, Grant Mickle, Robert Mitchell, and Cleveland Coe. 

Chf. Traffic Engr., National Safety Council, Inc., Chicago, 

Received the Secretary September 15, 1939. 

Journal General Psychology, April, 1939. 
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other directions are careful. Rather, the only conclusion that can drawn 
such cases that there something unusual and unanticipated about the 
condition which the drivers and pedestrians are subjected from certain 
approaches. Further evidence this found the fact that corrective 
measures often result the prevention these accidents simplifying the 
traffic and making easier drive walk safely. 

Mr. Mitchell’s interpretation the statement the belief 
many, city street not continuous scene human and mechanical 
somewhat different from what was intended, and some explanation should 
made. The purpose was not convey the impression that great many 
people over-estimate the seriousness the traffic-accident problem. That 
certainly not the case. was intended, rather, mean that many believe 
that accidents are rather indiscriminately distributed along the street and are 
likely occur almost anywhere. fact, however, they are concentrated 
certain points and along certain sections streets, whereas other sections are 
relatively free accidents. 

The difficulty with the belief many, the statement was intended 
convey, that would place the responsibility for accidents too largely upon 
abnormal unsocial drivers and pedestrians who are not interested their 
own safety. However, the patterns accidents individual 
locations, and the freedom from accidents long sections street bearing the 
same traffic that becomes involved accidents other places, show quite 
conclusively that substantial part the responsibility for accidents lies with 
the conditions peculiar individual locations. 

The very basis the engineering approach the solution traffic problems 
the belief that the vast majority drivers and pedestrians drive walk 
safely the extent that they can foresee and avoid hazards. The function 
traffic engineering make easier for drivers and pedestrians use the 
streets and highways safely and expeditiously either eliminating the hazards 
providing warnings, safeguards, and aids for avoiding them where they 
must remain. The more effectively this function performed, the easier will 
the functions traffic-safety education and enforcement. 

The field traffic engineering has the very fascinating job determining 
and evaluating driver and pedestrian practices important factors the 
design traffic facilities, control measures, and safeguards. The more 
quantitative data that are gathered concerning traffic the more obvious 
becomes that driver and pedestrian practices, under similar conditions, are 
approaching similarity. Vehicle speeds along given section street will fall 
into quite definite patterns, the dissenters from which are likely relatively 
few. Nearly all drivers will obey sign where the computed “critical 
speed,” determined view obstructions and other conditions, very low. 
However, where the critical speed relatively high, much smaller percentage 
will obey because the others resent restrictions which are obviously unnecessary. 
These and other facts which traffic engineers observe enable this branch 
the profession design measures for the prevention accidents and congestion 
and obtain anticipated results. 
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Another important function the traffic engineer make the 
his efforts reduce accidents congestion existing streets highways 
serve the basis for foreseeing and preventing similar conditions the design 
and construction new ones. the assembly that body information 
concerning the value different remedies for different problems and the 
conditions which differentiate one problem from another which will give his 
proposals for avoiding accidents congestion new street highway plans 
the prestige which comes from successful experience. 

Fortunately, setting new traffic-engineering department city, 
the engineer can start with treatment the many accident congestion 
problems that already confront him and work from them the treatment 
the problems which foresees the new streets highways that are being 
planned. time goes on, his job will change somewhat from the treatment 
existing problems the anticipation those that will develop the proper 
measures are not taken the design and construction improvements new 
streets highways. 
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DISCUSSIONS 


THEORY LIMIT DESIGN 


Discussion 


tional requirements any structure are: (1) That strong enough 
withstand the application the loads specified magnitude; and (2) that 
under the application these loads exhibit certain degree resistance 
deformation and thus maintain its geometric integrity. Strength and stiffness 
within the prescribed limits are the basic considerations any design theory 
whether applied individual members entire structures. Since the 
indices both the strength and the stiffness any structural member are 
expressed terms its dimensions and the physical properties the material 
which made, strength and stiffness structural member are definitely 
related. Parametric equations for strength and for stiffness can written 
each case and, varying the parameters, simultaneous values can obtained 
from which curves can drawn. From these curves the functional relation 
strength and stiffness can ascertained conveniently. the character 
this that determines whether the principle superposition (see head- 
ing applies. 

admit the application the principle superposition structure, the 
following conditions must postulated: 


(1) The character the impressed load must not change the load 
increased 


(2) The geometric configuration the structure must not appreciably 
change with the increase the load; 

(3) the Load changes, the axial (and bending) loads each member 
the structure must change proportionally the change Load 

(4) the loads the members change, with the change Load the 
unit stresses and total deformations must remain proportional the loads; and, 


paper Van den Broek, Am. Soc. E., was published February, 1939, 
Proceedings. Discussion this paper has appeared Proceedings, follows: May, 1939, Messrs. 
ohn Meursinge, Silverman, Edward Godfrey, Basil Sourochnikoff, Mirabelli, Goodrich, 
George Winter, and Francis Simpson; June, 1939, Messrs. Joseph Wise, Alfred Freudenthal, 
Hans Bleich, Alfred Niles, and Floris; and September, 1939, Messrs. Nishkian, and 
Erie Peterson. 
Associate Prof. Mechanics, Univ. Tennessee, Knoxville, Tenn. 
Received the Secretary August 1939. 
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(5) The total deformations the various members—tensile, flexural, and 
compressive—must the same order. 


Conditions (1) and (2) must also apply each member the 


The author’s predication (see heading structure 
will fail only when many members, restraints, the structure re. 
dundants plus one, are loaded their capacity strength” and his third dictum 
(heading are unquestionably valid within the limits the operation 
the principle superposition; his design procedure illustrated the 
application the design tower. the writer understands it, this design 
procedure consists three main steps which can outlined follows: (a) 
Apply over-stress factor, the design load structure; (b) the basis 
this load, draw continuous force diagram, and closes, 
may regarded satisfactory; and (c) assume the dimensions one member, 
conformance the force diagram, and proportion the neighboring members 
relation indicated the force diagram. 

design structure this manner the basis the configuration 
the structure when unloaded, assume that this configuration will not 
change appreciably while the capacity load gradually diminished zero value 


and then gradually increased back the portion its capacity value. 
assume that when this accomplished each member (which was designed for its 
capacity load) will subjected portion its capacity load, postu- 


late that throughout the process unloading and partial loading the loads 
each member, its unit-stresses and its total deformations, are constantly pro- 
portional Load must assumed that 


which the total load the structure; the load member; 
the unit stress the member; and the total deformation the 
member. This assumption equivalent stating that the structure must act 
that is, the structure acts whole and analogous single 
member like character. However, condition this nature can approxi- 
mated only extremely rigid structures, which the members each joint 
can assumed virtually equivalent rigidities, giving the structure 
character homogeneity. Their axial deformations, whether tensile, com- 
pressive, flexural, can assumed very small, and the resultant dis- 
placements the joints are negligible. The application the principle 
superposition can admitted and the design can made the basis 
static equilibrium. However, due attention permissible joint displacements 
and consequent deformations has resulted, the case suspension bridges, 
considerable economy Similar economy can undoubtedly shown 
for other structures. 


Transactions, Am. Soc. E., Vol. (1927), 894. 
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The stiffness structure whole function the ratio the loads 
(dead and live) that the structure must carry its over-all dimensions. 
certain types structures this ratio rather small and, the deflections are 
consequence and are pefmitted large, the design results 
type structure. The behavior these structures under applied loads 
very dissimilar the behavior structures heavy rigid type. This 
latter type will display degree elastic stability under any condition stress 
intensity its members the yield point. (In certain cases elastic stability 
may persist when the stresses some members are beyond the yield point.) 
The light type structure, however, may reach condition neutral equi- 
librium, presenting the danger failure buckling, when the magnitude 
stresses its members considerably below the elastic limit. This dissimilar- 
ity action can explained only the supposition that the ratios stresses 
deformations are not the same the two types. light structure the 
deformations proceed greater rate than the stresses, or, stated briefly, the 
principle superposition does not apply. Such light structures are very sensi- 
tive the effects the eccentricities and nonsymmetries the joints, and 
their light members will often subjected combination axial and lateral 
loads. 

When applied whole structures, the principle superposition must 
postulated from the proportionality total deformations the applied loads. 
the manner indicated the first paragraph this discussion, can shown 
that, general, the relation between the capacity loads various structural 
members and their stiffnesses not straight-linefunction. Itis asymptotic 
curve, similar many respects Euler’s curve, the asymptotic part which 
the range low slenderness ratios (that is, high stiffness) the members. 
Assuming this part the curve straight line can seen that the axial 
deformations heavy, stiff members, whether they are subjected tensile, 
flexural, compressive loads, will the same order. Experiments validate 
this assumption. the curve recedes from its asymptote the direction 
higher slenderness ratios, the rapid change the slope tangents the curve 
indicates that such similarity axial deformations slender members cannot 
assumed. other words, the principle superposition cannot applied 
structures joined slender members heterogeneous character. For 
similar reasons the principle superposition cannot applied slender 
members subjected the simultaneous action heterogeneous loads such 
compressive and bending loads. can shown that under these conditions 
the load-resisting function member, regards lateral loads, approxi- 
mately proportional minus the ratio the applied axial load the ca- 
pacity axial analogy similar proportionality can postulated for 
the structure whole. 

Where the conditions are such that the law superposition cannot ap- 
plied, the distribution loads various members cannot made the basis 
simple statical conditions—namely, the force diagram—but must made 
the basis Maxwell’s equations—that is, the method. 


Transactions, Am. Soc. E., Vol. (1930), 882. 
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Herein, the writer’s mind, lies the anomaly the author’s thesis. Although 
Professor Van den Broek emphasizes the importance deformations, the force 
diagram the key his analysis; although holds the deformations 
structure prime importance, unsound, the very methods 
analysis which are based considerations deflections. 

Just how far theory postulated upon the law superposition applies 
structures which the law does not apply, and within what limits deviation 
from this law, not problem pertinent this discussion limit design. 
theory pushed beyond its rational limitations becomes hypothesis. That the 
standard specifications, formulas and methods used structural design are open 
questioning cannot denied. They have been developed with reference 
certain types structures and meet specific problems. Their application 
the light type structures, particularly the highly indeterminate character, 
may questioned the grounds economy. this end, the author has 
rendered valuable service raising the problem. The theory limit design 
undoubtedly applicable certain cases. With modifications, particularly 
with reference compression members and the factor safety, perhaps its ap- 
plication could extended. The limits its application and the extent the 
necessary modifications can ascertained, however, only considerable 
additional analysis and checks against the methods that experience has 
shown reliable and economical. 
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HYDROLOGY THE GREAT LAKES 
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Discussion 


Symposium represent hydrologic balance sheet for the Great Lakes. The 
method used illustrated its application Lake Superior. The resulting 
values found for the individual items such balance have values which differ 
for some items from those more generally considered applicable similar 
conditions. The authors’ results can reviewed both from the point view 


their detail procedure and from general comparisons their results. 

Colonel Pettis states, rigid solution impossible where there are three 
unknowns and only two equations. meets this condition means 
assumptions which provide the equivalent athird equation. The correctness 
his conclusions thus depends the correctness these assumptions. 

Colonel Pettis assumes that the area supplying water Lake Superior 
limited its topographic drainage basin. This assumption the usual and 
natural one for this area. Its use requires that the sum all items water use 
must equal, and cannot exceed, the total rainfall the topographic basin 
tributary the lake. 

The known items are combined the expression Al- 
though all these items cannot measured completely exactly they are 
determinable within reasonable limits. The greatest uncertainty the 
item Colonel Pettis states the run-off only 49% the tributary drainage 
has been measured. 

shows negative values for April and May. This means that the values used 
for the rainfall the lake, the surface run-off into the lake and the lowering 
the lake are not sufficient supply the outflow from the lake these 
months. This, turn, means that some additional item inflow required 
supply the outflow even evaporation occurs from the lake. 


Symposium was published April, 1939, Proceedings. Discussion this paper has 
appeared Proceedings, follows: June, 1939, Young, Assoc. Am. Soc. E.; and September, 
1939, Messrs. Alfred Cooper, Jr., and Adolph Meyer. 


Prof., Irrig., Univ. California, Berkeley, Calif. 
Received the Secretary June 14, 1939. 
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Although the uncertainty the values for and may account for the 
-negative values found for all experience requires the inclusion 
some material amount evaporation these months. This, turn, 
results requiring additionally larger item inflow from some unmeasured 
source. 

Proceeding with the assumption that there unmeasured ground-water flow 
into the lake, Colonel Pettis derives values for its amount which produce the re- 
sults shown his summary balance for both land and water areas. The curve 
value the unknown evaporation from the lake and any other unknown inflow, 
indicates that either there are major omissions errors 
that other sources inflow must exist appreciable amount. also indi- 
cates that the values for evaporation months must have relative distribu- 
tion through the year that materially different from those that have been 
usually considered. the unknown item inflow assumed have fairly 
constant distribution through the year, the resulting values for evaporation 
from the lake must have winter values that are larger than the summer values. 

Using his assumptions Colonel Pettis completes the analysis secure the 
resulting values for land losses, ground-water flow into the lake and evaporation 
from the lake. Just Colonel Pettis used his results for 
judge the amount distribution months his unknown items, his con- 
clusions can used judge the reasonableness his results. 

The rainfall the drainage area must supply the land losses, the rainfall 
penetration which becomes the source the ground-water and underground 
flow into the streams and the lake well the surface run-off. For Lake 
Superior alone the items water supply must balance similarly. one item 
over-estimated some other item must reduced correspondingly the land 
balance. inflow the lake over-estimated, evaporation from the lake 
will have increased. 

Some the limitations the values that may secured for evaporation 
from Lake Superior can shown varying the assumed value the land 
losses for consumptive use. The total rainfall 28.2 in. may distributed 


Description Pettis For Various Assumed Values 


Consumptive use, inches depth land area. 5.9 10.0 15.0 16.6 
Inches Depth Land Area: 
Underground flow .................. 10.7 6.6 1.6 
Inches Depth Lake Area: 
Surface inflow lake, 17.9 17.9 17.9 
Underground flow into lake, 16.5 10.2 2.5 
Resulting evaporation from lake, 35.6 29.3 21.6 19.1 


various amounts illustrated Table These results illustrate the 
tions the different items. 16.6 in. the rainfall the land are used 
soil-moisture evaporation and plant transpiration, including rainfall intercep- 
tion vegetation, and Colonel Pettis’ estimate stream flow correct, there 
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remaining supply available for ground-water flow into the lake. meet 
Colonel Pettis’ conclusions land losses can only 5.9 in. and rainfall penetra- 
tion would have 10.7 in. plus the ground-water portion the stream flow. 

Although exact values can set for these items, Colonel Pettis’ conclu- 
sions differ materially from most the present values used for 
secures land loss 5.9 in. where water surface evaporation 35.6 in. land 
are one-sixth water losses. The drainage area Lake Superior 
contains many other lakes. Their total area not available; they represent 
one-sixth the drainage area, their evaporation would require all the water 
assigned Colonel Pettis land losses, leaving nothing for loss from the five- 
sixths the areainland. For the vegetative conditions this area, consump- 
tive uses should exceed 5.9 in. The rainfall sufficient support active 
vegetation. Cropped areas would require more than 6-in. loss for usual 
yields; forest areas would consume excess this amount. Colonel Pettis 
maximum monthly land loss less than in. per month. General 
values shown available records are well summarized Joseph Kittridge, 
who places the area surrounding Lake Superior the group annual water 
losses between and in. Mr. Kittridge found such losses less than 
in. only sagebrush and desert shrub areas. Although his method 
deriving results includes unmeasured deep ground-water movement with the 
land losses, his values are derived from areas adjacent Lake Superior where 
such deep ground-water movement not factor. 

The reasonableness the land-area balance can also approached from 
another point view. All the underground flow the lake, plus least 
part the surface run-off, must represent rainfall that has penetrated below 
the plant roots and reached the water table. The portion the surface run-off 
that represents ground-water flow into the streams not discussed numerical 
terms the Symposium and would difficult estimate closely. the 
stream flow fairly well sustained throughout the year such ground-water flow 
appears material part the measured stream flow. 

The quantity the derived underground flow into the lake represents 37% 
the precipitation. With ground-water movement into the surface streams, 
the summary presented Colonel Pettis would require that more than one-half 
the rainfall penetrate the ground-water. This much larger proportion 
rainfall penetration than has been found for other similar areas. The 
general results observations rainfall penetration the Great Plains area 
would indicate little penetration from the 28.2 in. rainfall heavy soils, 
perhaps in. sandy lands conditions are favorable the drainage 
the Loup River Nebraska. The probable amount rainfall penetration 
the drainage area Lake Superior would not appear sufficient account 
for more than the ground-water flow into the streams which included the 
stream item Colonel Pettis’ water supply balance. 

The small value found for consumptive use requires correspondingly high 
value for the total run-off. For both surface and underground run-off, the 


Magnitude and Regional Distribution Water Losses Influenced Joseph 
Jr., Journal Forestry, Vol. XXXVI, No. August, 1938, pp. 775-778. 


Ground-Water from Rainfall,” Harding, Am. Soc. E., Trans- 
actions, Am. Geophysical Union, 18th Annual Meeting, 1937, pp. 368-371. 


{ 
, 
y 
’ = 
l- 
in ; 


1438 HARDING HYDROLOGY THE GREAT LAKES 


total run-off derived 22.3 in. from rainfall 28.2 in. This exceeds the 
quantity run-off obtained from adjacent similar areas. Areas farther west 
the same latitude, with about in. rainfall, have surface run-off only 
in., and have little rainfall penetration. 

Evaporation.—The Symposium presents two determinations evaporation 
from Lake Superior, separately derived. Their close agreement does not neces- 
sarily represent proof the correctness either. 

Mr. Hickman’s method observation evaporation represents com- 
mendable effort eliminate one the factors which causes observations with 
pans differ from the loss from large water areas. The method used, however, 
raises additional questions regarding the interpretation the pan results 
terms lake evaporation. 

The maintenance the temperature the water the pans the same 
value the temperature the surface water the lake does not itself 
assure the same rate evaporation. Although may reduce the extent 
the difference, comparisons between the temperature-controlled pans and 
measured evaporation from large areas are needed before the assumption that 
the controlled pan has coefficient unity all months can justified. 
Although the difference temperature between the water the pans and 
lakes one the important factors their different rates evaporation, other 
factors such rim effects the pan, differences exposure and differences 
pan size still need considered. 

Evaporation dependent the available heat supply. the lake the 
heat supply primarily that received from solar radiation. The evaporation 
will less than the quantity water that the total solar radiation can convert 
into vapor because back radiation and any loss heat surrounding bodies. 
Heat stored the lake during months rising water temperatures and 
withdrawn during those lowering temperatures. The authors state that 
temperature variations extend depth They present the monthly 
mean surface temperatures but not show the monthly temperature profiles 
the lake. 

the heat added the pans during the cooling months, maintain the 
same water temperature the pan the surface water the lake, equal 
the heat obtained the lake surface the cooling the lake water, then 
the loss from the pan would represent the probable evaporation from the lake 
far the heat balance items are involved. the summer months the heat 
removed from the water the pan would need equal the heat stored 
the water the lake. Without records the heat added removed from 
the pan, temperature profiles the lake, comparisons the basis 
heat balance cannot made. 

Mr. Hickman’s results differ from direct observations other lakes which 
not freeze during the winter. some western lakes conditions permit the 
same type analysis used Colonel Pettis where ground-water inflow 
outflow negligible amount. Among such lakes areas fairly low winter 
temperatures are Pyramid and Walker lakes Nevada. These lakes also 
not freeze during the winter; both have sufficient depth comparable, 


are 

exe Al 

De 

i e\ 


October, HYDROLOGY THE GREAT LAKES 1439 


this respect, Lake Superior. Results show much less relative 
evaporation during the winter months than are found the authors. The 
evaporation and air temperatures the three lakes are compared Table 
For the air temperatures the Weather Bureau records the nearest station 
are used, Lahontan being used for Pyramid Lake and Schurz for Walker Lake. 


TABLE 7.—CoMPARISON EVAPORATION AND AIR TEMPERATURES; PYRAMID 
AND WALKER LAKES NEVADA AND LAKE SUPERIOR 


Pyramid Lake Walker Lake Lake Superior Percent the Annual Total 


Mean air Mean air Mean air 
Evapo- temper- Evapo- temper- Evapo- temper- 
ration, ration, ration, ature, in| Pyramid| Walker Lake 
degrees degrees degrees Lake Superior 
feet Fahren- Fahren- inches Fahren- 


January 
February 
March 
April 
May 
June 

July 
August 
September 
October 
November 
December 


Total 


Until the heat-controlled pan can tested adjacent large water surface 
for which the evaporation can determined, the coefficient that needs 
applied the results with the heat-controlled pan must remain unknown. 
Although the coefficient might prove unity, both for annual and for 
monthly periods, the writer will venture the prediction that the loss from the 
heat-controlled pan will found exceed that from the adjacent large water 
area and that the monthly ratios the pan the large area will variable. 
The size the heat-controlled pan also factor its amount evaporation 
performance similar other types pan. Comparisons with different 
sizes heat-controlled pans would not difficult make; size were found 
factor would demonstrate that the assumption value unity for 
the coefficient their pan used the authors was not justified. 

Pettis meets the apparent contradiction negative 
evaporation April and May, the ground-water inflow Lake Superior 
taken zero, assumptions which remove that difficulty but lead im- 
probable values for land losses and rainfall penetration. Mr. Hickman’s heat- 
controlled evaporation pans give values for evaporation which agree with 
Colonel Pettis’ results but, applied Lake Superior without correction, lead 
the same difficulties land losses that result from Colonel Pettis’ analysis. 

Although these two papers represent commendable effort determine the 
hydrologic balance for the drainage area Lake Superior, the problem does 


ng, Transactions, Am. Geophysical Union, 16th Annual Meeting, 1935, pp. 
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not appear Further work needed appears include: (1) more 
direct study land losses determine what part the rainfall may 
able for ground-water movement into the lake; (2) perhaps review and 
ble revision the estimated surface flow from the unmeasured portions the 
drainage area; (3) comparisons the heat-controlled pans with large water 
areas whose evaporation can measured; and (4) monthly temperature pro- 
files the lake determine cooling the late fall and winter months can 
supply the amount heat for the estimated amount evaporation 
these months. Although previous estimates evaporation from Lake Superior 
may have been too small, Mr. Hickman’s estimate, when all factors are con- 
sidered, appears too high. Actually the evaporation from the lake must 
have value that can supplied from the rainfall and leave consistent water 
supply available for land losses. Mr. Hickman’s results not appear meet 
this standard. 


not discuss the geology and ground-water conditions the drainage basin 
tributary Lake Superior. The area known underlain, for the most 
part, pre-Cambrian crystalline rocks which general are impervious 
have only low permeability. most places these hard old rocks are overlain 
thin mantle glacial drift which is, part, water-bearing. 

far known, the rocks this basin not, except locally, have any 
artesian structure other structural features for conveying the ground-water 
laterally toward the lake. The ground-water believed essentially under 
water-table conditions, with the water-table sloping toward the streams and 
branches and not generally toward the lake. these conditions prevail, large 
part the stream flow doubtless derived from ground-water, but the quantity 
water percolating directly into the lake can only small fraction the 
10.7 in. postulated the paper. 

the basis the great mass information that available regard 
the geology and ground-water conditions the United States, can stated 
confidently that, for most the large drainage basins the United States, the 
loss gain water subterranean flow relatively small, and the discharge 
transpiration and evaporation approximates the difference between precipi- 
tation and stream discharge. general, the rate discharge transpiration 
and evaporation changes gradually and systematically from one region 
another with changes climate. 

Colonel Pettis gives the land losses for the drainage basin Lake Superior 
5.9in. Ifa larger rate loss conceded for the extensive swamp areas 
and numerous minor lakes within the basin, the rate loss from the actual 
land areas would still smaller. 10.7 in., the purported value 
added 5.9 in., the purported value the total 16.6 in., which ap- 
proximately the discharge transpiration and evaporation that would 
expected the Lake Superior basin from the known regional loss. 


Geologist charge Division Ground Water, Water Resources Branch, Geological 
Washington, 


Received the Secretary August 18, 1939. 
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The writer has not been able get clear understanding the process 
which the evaporation from the lake has been estimated computed 
this paper. apparently quite conjectural. The value 10.7 in. for 
ground-water inflow (U) has been derived merely subtracting from 
Furthermore, the value 5.9 in. for the land loss (L) has been 
derived merely subtracting Ufrom Apparently effort was made 
study the actual ground-water conditions. 

would seem that, before the results 10.7 in. and 5.9 in. can given 
favorable consideration, systematic ground-water survey the drainage basin 
should made, including quantitative ground-water study least one 
locality that appears relatively favorable for percolation ground-water into 
the lake. The quantitative study should include such items contour map 
the water-table and determination the extent, thickness, and coefficient 
permeability the water-bearing beds. 


Colonel Pettis reaches the conclusion that, for the drainage area the Great 
Lakes, the evaporation and underground flow are greater, and the land losses 
are less, than the values that have been generally accepted. 

mean annual evaporation from Lake Superior has been estimated 
notable investigations recent years follows: in. Mr. Freeman, in. 
Messrs. Horton and Grunsky, and in. Messrs. Hayford and Folse. 
These results contrast with 35.6 in. obtained the authors the Symposium. 
The substantial excess the latter value estimated evaporation over the 
others impressive. 

Unquestionably, the mean annual rainfall and the mean annual run-off 
from the land area, 28.2 and respectively, obtained Colonel Pettis, 
may accepted reasonably reliable. The results for outflow from Lake 
Superior and changes storage the lake are determined methods which 
are reasonably precise and should correspondingly reliable. Consequently, 
Colonel Pettis’ further conclusion that underground flow more and losses 
are less than generally accepted seems rest fundamentally the validity 
the values adopted him for evaporation from the lake. 

The writer has participated rather extensive studies so-called land 
(or water) losses, the difference between rainfall and run-off, Eastern 
United States. (The term used herein includes all forms 
precipitation.) The studies were supervised Gordon Williams, Assoc. 
Am. Soc. E., under the general direction the 

They were limited streams the humid semi-arid regions the United 
States east the Rocky Mountains. The basic period for which the computa- 
tions were made the water year, year ending September 30. Only records 
which had continued without interruption for three years more were used. 
The basins were also selected with view to: (a) The avoidance complications 
due underground flow, indicated large springs other evidence 

Hydr. Engr., Geological Survey, Washington, 

Received the Secretary August 18, 1939. 


Losses the United Gordon Williams and others, Water-Supply Paper No. 
846, Geological Survey (publication pending). 
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material ground-water influence; (b) the avoidance artificial regulation 
storage; and (c) the availability sufficient number precipitation stations 
for the reliable determination the average precipitation the 
basins. The gives results for about 200 areas with aggregate period 
record more than 000 years, average about ten years per record, 


St. AND ADJOINING REGION 


MEAN ANNUAL 
TITIES, INCHES 


Drainage 
Period 
Basin above the indicated gaging station studied, 
square 
water years Pre- 
miles cipita- Run- Water 
tion off loss 
Streams Tributary Lake Michigan: 
Thornapple River near Caledonia, Mich. 1932-1934 773 32.0 9.3 22.7 
Muskegon River Newaygo, Mich..... 1932-1934 450 30.1 10.2 19.9 
Streams Tributary Lake Huron: 
Tittabawassee River Freeland, 1916-1920 530 29.7 9.3 20.4 
1932-1934 
Streams Tributary Lake Erie: 
River Rouge Detroit, Mich.......... 1932-1934 194 28.6 6.0 22.6 
Huron River Barton, Mich.......... 1915-1920 723 31.7 9.2 22.5 
(b) 
Chippewa River Basin: 
Jump River Sheldon, Wis........... 510 30.5 12.9 17.6 
Trempealeau River Basin: 
Trempealeau River Dodge, 1915-1919 633 29.5 8.3 21.2 
Black River Basin: 
Black River Neillsville, 1915-1934 756 31.1 9.6 21.5 
Crosse River Basin: 
Crosse River near West Salem, Wis. 1915-1934 412 30.3 10.0 20.4 
Wisconsin River Basin: 
10 Rib River at Rib Falls, Wis............ 1926-1934 309 30.1 12.4 17.7 
Yellow River Sprague, Wis.......... 1927-1934 420 28.8 6.3 22.5 
Kickapoo River Gays Mills, 1915-1933 ..... 31.6 9.3 22.3 
Rock River Basin: 
Sugar River near Brodhead, 1915-1934 529 32.5 9.2 23.3 
(c) Hupson River 
Red River Grand Forks, Dak......... 1882-1934 500 20.9 19.7 
99-199 
Red Lake River Crookston, Minn........ 320 19.5 1.6 17.9 


Table shows summary precipitation, run-off, and water losses for 
areas the western St. Lawrence River Basin and the adjoining region, they 
appear the report. The water-loss results are plotted Fig. order 
facilitate consideration relation their several geographical positions. The 
requirements the study for continuity and duration the records, freedom 
from complications caused artificial regulation stream flow, and adequacy 
the rainfall data caused the distribution selected basins sparse the 
region indicated. 


‘Mean Annual Water Loss, in Inches 
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used these studies, the run-off the surface stream flow measured 
gaging station, and the rainfall the computed average precipitation over the 
drainage area tributary the gaging station. The difference between these 
two factors, course, represents the evaporation and transpiration from the 
basin plus ground-water outflow from the basin that not measured the 
gaging station, minus ground-water accretion from sources outside the basin 
which appears the measured surface flow the gaging station. The common 
designation this difference water losses, representing the evaporation and 
transpiration, misnomer the extent that underground flow involved 
the ways indicated. Adjustments the water-loss results computed would 
subtractive for ground-water outflow that does not appear the gaging 
station and additive for ground-water accretion from outside the basin that 
appears the gaging station. Variations evaporation, transpiration, and 


a 
< 


Mean Annual Temperature, in Degrees Fahrenheit 


Fic. MEAN ANNUAL Loss AND MEAN ANNUAL TEMPERATURE FOR 
WITH MEAN ANNUAL PRECIPITATION Excess INCHES 


ground-water effects, both additive and subtractive, would reflected 
variations the computed difference between rainfall and run-off. 

The water losses show general uniformity under similar climatic conditions 
and tendency systematic transition value these conditions change. 
striking feature the progressive decrease the computed loss direction 
from south north. the Gulf States computed losses in. in. 
more are typical, whereas the latitude central New England and Wisconsin 
the losses seem range generally from in. 24in. the aforementioned 
studies data were computed for comparing the mean annual water loss and the 
mean annual temperature for each drainage areas selected, with view 
adequacy basic data and representative distribution. 

points plotted show the respective basin determinations and lines enveloping 
them. shows that for these basins the water losses vary general with 
the mean annual temperature. Unless one makes the unreasonable assumption 
that all these basins are losing large parts their water underground flow, 
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must conclude that Fig. shows approximately the losses evaporation 


and transpiration relation the mean annual temperature. 
The mean annual temperature the area tributary Lake Superior 
the order degrees. the relation indicated Fig. applicable, 


gPort Arthur 
William 


Two Harbors 


N 


would appear that the normal annual water loss evaporation and transpira- 
This compares with 16.6 in., obtained 


tion between in. and in. 
taking the difference between the rainfall and run-off results given Colonel 


Pettis, and differs radically from the 5.9 in. computed him. 
Colonel Pettis’ determination evaporation seems likely source 
discrepancy, and hoped that the work that factor may receive critical 


review. 
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DISCUSSIONS 


RECONSTRUCTION THE WALPOLE-BELLOWS 
FALLS ARCH BRIDGE 


Discussion 


between Bellows Falls, Vt., and Walpole, H.—a 540-ft arch span across the 
Connecticut River, built 1904—was ambitious project years ahead its 
time. When was planned, the Advisory Engineer, the late Snow, 
Am. Soc. E., and the Consulting Engineer, Worcester, Hon. Am. 
Soc. E., both recognized the advantage keeping the structure above 
possible ice damage; but the Bridge Committee, composed local business 
men from both sides the river, would not hear the extra expense. 

1906, year after the bridge was finished, flood damaged the bottom 
lateral bracing one end and the consulting engineer recommended that 
reinforced concrete slab installed brace the lower chord and protect the 
lower laterals, but there was authority hand see that this was done. 

1928, the dam was raised, making the pool higher; the canal 
was enlarged and the power-house expanded. was inevitable that sooner 
later this elevation the dam would cause destructive flood. 

this time some strong State agency could have bargained with the 
power company provide the necessary protection around the abutments 
partial offset for the benefit increased power production. usually 
the case, the local authorities were incapable handling these matters against 
power corporation, especially the bridge was interstate project. 
Had these precautions been taken, the heavy expense indicated the paper 
might have been avoided. Fortunately, the entire bridge was not lost and 
the State New Hampshire was able take charge and repair it. 

Offhand would seem that the cheapest method supporting the bridge 
during repairs would have been use falsework towers was done originally. 
far the writer knows the conditions were the same when the bridge 


paper Messrs. Langley and Edward Ducey was published April, 1939, 
Proceedings. Discussion this paper has appeared Proceedings, follows: September, 1939, 
French, Am. Soc. 

Roselle, 

Received the Secretary August 14, 1939. 


| 
‘ 
. 
} 
. 
‘ 
‘ 


1446 RIGHTS WALPOLE-BELLOWS FALLS BRIDGE Discussions 


was first erected, except that the water was deeper the pool. This would 
mean longer, and possibly few more, piles but would have been balanced 
shorter timbers the framing the towers. 

the work had been done after the spring rise, when the ice was gone 
and before freezing weather, there would have been very little risk, even 
summer flood had occurred. Assuming that the erection cost the bridge 
originally was one-half the cost the superstructure, $20 000, and taking 
into account the new repair material, difficult understand the lump-sum 
bid $120000. the ‘‘cost plan the State assumed the risk from 
operating hazards, and this had been done for the lump-sum plan this bid 
would probably have been much cheaper. Allowing for the changes the 
abutments would still leave difference $50 000 $75 000 explained. 

The writer hopes that the authors will give the costs the repairs detail; 
that is, the abutments, new repair steel, erection including cables and the new 
floor well information what was covered the lump-sum bid. 

All this does not the least detract from the high praise due the engineers 
for the skill whereby the span was shortened from 540 486 and the 
spring-line pins raised ft. The success their scheme the result bold 
design, careful planning, and constant care. 

The original bridge was erected upper New England, over capricious 
river, the depth winter because the delays the Committee awarding 
the contract. The bridge was repaired, also, the winter months, probably 
glory engineers that they are willing accept difficult and hazardous condi- 
tions forced circumstances beyond their control, without wasting time 


quibbling about them, and that they then proceed overcome these hazards 
successfully. 
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AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


DISCUSSIONS 


DESIGN OPEN-CHANNEL 
CONTROL SECTION 


Discussion 


Frank Assoc. Am. Soc. (by author 
seems have assigned himself difficult task and performed very ingeniously. 
The writer’s first impression this paper was that was too complicated for 
his understanding and possibly this case the first impression was the best. 
However, during leisure hour began trying follow the author’s reasoning 
and check the equations has derived. 

Becoming puzzled over Equation (6) wrote letter inquiry the 
author, who replied courteously: 


not think necessary call attention the other mathematical possibili- 
have practical significance our problem, for example 
when 


the writer has not found why must so. believes that 


and that 


Equation (5), obtained differentiating Equation (4), could read 


dA. 


which signifies the derivative f(A-) with respect Perhaps 
Equation (5) just valid; but Equation (47) seems more complete. Re- 


paper Kar! Kennison, Am. Soc. E., was published May, 1939, Proceedings. 
Discussion this paper has appeared Proceedings, follows: September, 1939, Messrs. Streeter, 
Emery Willes, and Robert Thomas. 


Braintree, Mass. 
Received the Secretary July 29, 1939. 
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which the derivative Therefore, placing the numerator 
tion (47) equal zero, the relation 

(46b). 
the propriety having factor 
the denominator which equals zero: 
the numerator must equal zero sat- 
isfy Equation (5) the denominator 
also zero seems that the value the 
fraction becomes which generally in- 
determinate, although not always. any event factors the 
nator must, necessity, equal zero, the denominator may have 
almost any value and still satisfy the requirement for maximum value 


Horace King, Am. Soc. E., describes the formula 


general criterion, applicable all channels, indicating the relation be- 
tween cross-section and discharge which must exist order that flow may occur 
critical Equation (48) resolves 


and identical with Equation (46b) being the top width 


the cross-section the stream the control section. 

The writer has made rough calculation the area corresponding 
flow 28.3 per sec, using values given for and Table 
this calculation was assumed that the sides the section Fig. were 
straight between the the sides are curved, the value 3.34 
ft, obtained for little too large. 

this value used Equation (48), the resulting numerical expression 
34.1 24.9, which does not fulfil the requirements 
the formula, the two sides should equal. the precise value 
found very close 3.00 (which probably the case), Equation (48) would 
check exactly and would proof the correctness Mr. Kennison’s results. 


10‘*Handbook Hydraulics,” Horace Williams King, McGraw-Hill Book Company, Inc., 
York and London, Third Ed. (1939), 373. 


c 
8 
: 


AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


DISCUSSIONS 


TENSION TESTS LARGE RIVETED JOINTS 


Discussion 


poses few problems and tests, and suggests few changes specifications 
the interpretation the present ones. 

rivet the fundamental element riveted joint. The 
two main variables rivets are the grip and diameter. suggested that 
there is, should be, some relation between small pin and rivet design. Fig. 
shows the relation between bearing, shear, and bending based simple 
the bending only half much assumed Fig. 20, still appears 
very important factor. The tests reported the paper indicate that 
decrease strength can expected the grip increases, but the tests 
have very limited variation grip. Fatigue tests the University 
Illinois support this conclusion. German fatigue tests have developed breaks 
bending. 

The writer suggests that the specifications long rivets made more 
conservative; for example, (long rivets) subjected calculated stress 
and having grip excess two diameters shall increased number 
least 16% per inch grip.” Tests isolate this problem should relatively 
simple. 

Fig. appears that, for thin plates, bearing important. These 
tests were made relatively thick plates that bearing for thin plates was 
not revealed. There appears considerable evidence that the allowable bearing 
stress can quite high. probable that the rivet heads confine and assist 
the joint distributing the load over larger area. 

End tests reported the paper did not indicate much need 
for greater end distances, but these tests were for relatively thick plates. One 
proposed specification calls for the end distance, edge distance, such 
give area equal the shearing area the rivet. appears reason- 


paper Raymond Davis, and Glenn Woodruff, Members, Am. Soc. E., and 
Harmer Davis, Assoc. Am. Soc. E., was published May, 1939, Proceedings. Discussion this 
has appeared Proceedings, follows: September, 1939, Messrs. Charles Goodrich, Frederick 

Shearwood, and Jonathan Jones. 
Associate Bridge Engr., State Div. Highways, Bridge Dept., Sacramento, Calif. 
Received the Secretary August 14, 1939. 
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1450 WINTER TESTS RIVETED JOINTS Discussions 
able demand that this area maintained between rivets prevent 
Proposed: ‘‘The pitch should not less than the end distance plus one-half 
the diameter the 

the “B” series, where full-row riveting was used, there 
indication that the maximum efficiency can obtained when the pitch 
two-thirds the gage. This point can possibly explained the 
tion the Lueders lines. the lines are 45°, then with the pitch two- 
thirds the gage, the maximum stress lines will fall between rivets and not 


Bearing = 4> 30000 Lb per Sq In. 
Shear =f,=15000 Lb per Sq In. 
Flexure =f,=30000 Lb per Sq In. 


fo) 


Bearing 


G=Total Grip, in 


through them. Photo-elastic studies should indicate the most efficient rivet 
pattern and the effect pitch and gage the stress the plates. 

Net Section.—These tests raised large question about net section. Under 
the heading ‘‘Conclusions from Tests: Net Section” the following statement 
made: experimental results bore consistent relation the specifica- 
tion requirements for determining net 

The writer concurs that the tests not agree with common practice and 
methods shown textbooks. However, net section, used the analysis 
these tests and common practice, appears contrary strict inter- 
pretation present specifications, which state that: net section shall 
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the least area which can obtained deducting from the gross sectional 
area, the area holes cut any straight zigzag section across the member, 
the foregoing specification permits the protection the net section 
omitting rivets the first row, then the specifications need revision. 

The grooving effect important when full-row riveting used. Why not 
take advantage this property and use the nominal diameter the rivet 
when the holes are subpunched and reamed? 


TABLE AREA AND NET AREA 


T' 
Specimen 
(see Figs. 
and 4) 
Protected area* Least areat Protected net area 


6.96 
6.93 
CCC2-3-1 6.94 
CCC2-3-2 


aon 


© 


OOD 


DCC2-3-1 
DCC2-3-2 
DCC2-4-1 
DCC2-5-2 
DCC2-6-2 
DCC5-1 

DCC5-2 

DCC7-1 

DCC7-2 


FSCA-1 


FSCB-1 
FSCB-2 


FSCD-1 
FSCD-2 


FSCE-1 
FSCE-2 


FSCF-1 
FSCF-2 


6. 


ONE 


oo wo 

an 


OO 


Protected net area the first row net area used the tests. 

The least net area based the nominal diameter the rivet, deducting one hole for each gage line. 

Table gives the ultimate strengths based the nominal diameter the 
one hole being deducted for each gage line, resulting the least net 
section obtainable. found that the least net section fair agreement 
the ultimate strengths based mill coupons. study the data from 


69.6 
69.6 
68.2 
68.2 
69.4 
CCC2-4-2 69.4 
70.0 
70.0 
CCC2-6-1 69.9 
CCC2-6-2 69.9 
67.9 
67.9 
CCC7-1 67.9 
67.7 
67.7 
68.1 
68.1 
11 68.1 ‘ 
68.1 
67.7 
67.7 
67.2 
67.2 
9.28 9.67 67.6 
9.29 9.67 67.6 
8.72 67.8 
8.18 8.72 67.8 
12.37 10.73 91.5 
10.81 91.5 
24.82 21.55 91.9 
91.9 
37.07 32.18 91.8 
24.76 21.51 90.4 
21.60 90.4 
36.15 31.40 92.5 
36.93 32.05 
36.87 32.00 90.3 
37.13 32.22 90.3 
‘ 
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the tests the Bureau Standards also gives more 
able results when based the least net section. 

One very desirable quality any specification basis upon which the 
ultimate strength can predicted. The present tests not have large 
riveted joints with full-row riveting and seems desirable perform these 
tests before the specifications limit the riveted joints 75% efficiency. set 
comparable butt-welded joints might also worthy consideration. 

Proposed: calculating the required net section riveted tension 
bers, net sections shall used all cases; and deducting rivet holes, the 
holes shall taken as: (a) in. larger than the nominal diameter the rivet 
for punched holes; the nominal diameter the rivet for reamed drilled 
holes. determining the net width, one hole shall deducted for each 
gage 

Transverse regard the findings transverse shear the 
rivets, the tests indicate many coarse square fractures with very little reduction 
area. there transverse shear the rivets, there must transverse 
tension the main plates and compression the cover plates. These stresses 
probably contributed the curling and buckling the plates. combined 
stress appears have caused the square breaks. The question might put: 
the effect yield, ultimate strengths and reduction area the 
lateral contraction 

The tests seem indicate the desirability using full-row riveting, short 
compact joints, and ductile materials. 

appears that co-ordinated program tests should undertaken and 
then uniform specification written. 


Am. Soc. (by the writer, the two 
most important conclusions from the tests given the authors are: 


(1) The ultimate strength large group driven rivets equals, approxi- 
mately, the sum the strengths the individual rivets; and, 

(2) Omitting alternate rivets from the outside row group rivets, 
increasing the transverse distance between the rivets, does not add greatly 
the effective net section riveted tension member. Conclusion (1) 
reassuring and highly gratifying. Conclusion (2) has been indicated other 
tests, and will the subject the present discussion. 


The specifications most used the design American steel bridges contain 
rules for determining the effective net section riveted tension members. 
These rules assign considerable importance the relation between pitch and 
gage some rivets are omitted from the outside row and not recognize the 
effect the transverse spacing the rivets upon the unit strength the 
section through the rivet holes. (The outside row rivets joint corre- 
sponds the outside row gusset plate the case riveted tension 
member connected gusset plate.) This true although there has been 


Investigation the Behavior and the Ultimate Strength Riveted Joints Under Load,” 
Gayhart, Transactions, Soc. Naval Architects and Marine Engrs., Vol. 34, 1926, 55. 


Research Prof., Structural Eng., Univ. Urbana, 
Received the Secretary September 18, 1032. 
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considerable quantity experimental data available for several years, indi- 
cating that the ratio pitch gage for staggered rivets minor factor and 
that the transverse spacing the rivets major factor affecting the strength 
riveted tension members. The tests reported this paper are value 
that they add the accumulating evidence that present specifications are 
error both these respects. 

The writer has searched American literature obtain experimental data 
used study the effective net section riveted tension members such 
are used for truss members steel bridges. With few exceptions, the tests 
included this study comply with the following requirements: (1) The width 
the specimens varied from 13.5 in. 22.5 in., inclusive; (2) the strength 
the plate material given tests standard control specimens was 
reported; (3) the specimens failed the plate; (4) the specimens had either 
Lin. l-in. rivets; and, (5) the joints were generally accord with what 
accepted good design practice. However, there were the following exceptions 
the requirements: One specimen was only in. wide; four specimens failed 
the rivets but plate stress that made the tests value. Forty-nine tests 
reported four investigators complied with these requirements, except 
noted, and have been used the study (they included eight tests made the 
Watertown twenty-one tests made for the four tests Mr. 
and, sixteen tests the University California and 
reported the paper being discussed). addition the foregoing, number 
tests specimens in. wide, made Mr. have been used the 
study the effect the ratio pitch gage upon the effective net section. 


(1) 


Plate 


(3) Tests; Holes Punched ® 


Ultimate Strength, in Kips 
ny 
> 
o 


Ultimate Strength, = 

0.50 1.00 1.50 2.00 per In. 
Ratio Pitch Gage, 


Preliminary the formulation rule for determining the effective net 
section riveted tension members, the writer made one study determine 
the effect the ratio pitch gage the strength joints from which some 
rivets were omitted from the outer row, and made second study deter- 
mine the effect the transverse distance between rivets upon the unit strength 
the net section. These studies are presented herein. 

The tests used the study determine the effect the ratio pitch 
gage upon the strength joint from which rivets were omitted from the outer 


Investigation the Behavior and the Ultimate Strength Riveted Joints Under 
Gayhart, Transactions, Soc. Naval Architects and Marine Engrs., Vol. 34, 1926, 55. 


1930 partial fulfillment the requirements for the degree Master Science. 
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row include those made Mr. Harris, some made for the Navy, and 
some made the University California. The tests made Mr. Harris, 
and the details the joint, are shown Fig. series tests riveted 
joints are reported: For one the holes the plates were punched (Curve 
and for the other the holes were drilled (Curve (2)). third series tests was 
made Mr. Harris for which the specimens were plates without joints, but 
with open holes the same pattern the rivets the specimens with riveted 
joints. The results these latter tests are shown Curve (1) 21, 
Tests were made duplicate for each series and the small circles represent the 
average the results tests duplicate specimens. The ordinates the 
diagrams represent the ultimate load and the abscissas represent the ratio the 
pitch the gage. The strength the plate material was determined tests 
standard control specimens for all Mr. Harris’ tests, and the ultimate 
strength for all specimens has been adjusted common plate strength 
000 noted that the plates with open holes tested 
Mr. Harris were stronger than the plates connected with rivets. The relation 
between the strength and the stagger, computed the American Railway 
Engineering Association specifications, given Curve (1), Fig. 21. 

Tests made for the Navy and used this study, 
include six each with the rivet patterns shown Table group 


TABLE 27.—INFLUENCE RELATION PitcH GAGE Upon 


1" Rivets 


3" 


a 


Description 


Item No. 


Me- High- Me- High- 
dium test dium test 
steel steel steel steel 


Strength control specimen, pounds per square inch...| 760 64975 320 
Effective net area, test, square 13.79 13.11 
Net Section American Railw Engineering ‘Association 

specifications, square inches. 12.31 12.38 13.58 13.80 
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six specimens included three with carbon-steel plates and three with low-alloy 
plates. For the specimens with carbon-steel plates, the efficiency the joints, 
the average three tests each instance, was 0.84 for Type and 0.81 for 
Type For the specimens with 
low-alloy plates the corresponding Rivets 
efficiencies were 0.82 and 0.80. That 
is, changing the pitch from 
the gage decreased the efficiency 
the joint only about per cent. 
(The term efficiency used herein 
the ratio the ultimate strength 
determined test the product 
the gross area the section and 
the unit strength the control 
specimen.) 

Two series tests made the 
ported the authors (Specimens Plate 
No. CCC2 and Fig. and 
Table 5), had the rivet pattern shown 
Fig. 22. The plate thickness 
was in. for Specimen CCC2 and in. for Specimen DCC2. The pitch, had 
for The values the efficiency the joints for the various values 
large consistent change the efficiency the joints. 


Net TENSION MEMBERS. 


CCC2 DCC2 


Description 
inches inches inches inches 
Ultimate load, thousands 
pounds per square 399. 388.0 28. 20. 691.0 |686.0 
Stress gross section, 
pounds per square inch. -| 51. 49.9 of 53.¢ 52. 52. 52.7 51.6 
Strength control specimen, 
Effective net section, test, 
Net section American 
Railway Engineering As- 
sociation specifications, 


7.17 7.13 10.97 11.62 11.57 
0.790 0.768 0.804 0.803 0.793 


The results the three independent agree one impor- 
tant respect: The strength riveted joints, with some rivets omitted from the 
outer row, not greatly affected the ratio the pitch the gage. 
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the study determine the effect the transverse distance between rivets 
the unit strength the net section the plates, only two types joints were 
included: (1) Joints with rivets omitted from the outer row; and (2), joints 
with some rivets omitted from the outer row but only those which had pitch 
equal to, greater than, the gage. They included tests from all four 

gations. The relation between the ratio, Strength control specimen’ and 
the transverse distance between rivets shown the small circles Fig. 


Strength Net Section 
Strength of Control Specimen 


io 
a 


Ratio, 


Transverse Distance Between Rivets, in Inches 


The numeral adjacent each circle indicates the number tests identical 
specimens averaged. The study includes twenty-eight tests and variation 
gage from in. 8in. The points marked and are for specimens that 
failed the rivets. apparent from Fig. that the unit strength the 
net section riveted tension member greatly affected the transverse 
distance between the rivets. 

view the discussion contained the previous paragraphs, appears 
that any rule for determining the effective net section riveted tension 
member should take into consideration the net section through the outer row 
rivets and also the transverse distance between rivets. Furthermore, 
would appear unnecessary consider the effect minor variations the 
ratio the pitch the gage for joints with rivets omitted from the outer row. 
With these ideas mind the following empirical rule for determining the 
effective net section has been developed. 

The effective net width (W,) riveted tension member given the 
equation 


Net Section 
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Section 


Section 


e 
Gross 


8 
Unit Strength of Net Section 
Unit Strength of Control Section 


1457 
0.76 
0.72 
Ratio, 
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which gross width the member; nominal diameter the rivet 
plus in.; number rivet holes deducted; and, transverse 


distance between rivets, two times the edge distance, whichever the 
greater. 


rivets are omitted from the outer row, and shall determined 
the rivets the outer row provided the pitch equals exceeds 0.80 the 
gage; otherwise the intervening rivets the second row and all the rivets 
the outer row will considered determining and 

Equation (4) composed two parts: One the actual net width, and the 
other factor compensate for the effect the transverse distance between 

Unit strength net section 
rivets upon the ratio, 


factor for various values are shown Fig. 24. 


The proposed rule for determining the effective net section riveted 
tension members has been applied the forty-nine tests the four investi- 
gations the beginning this discussion and the results are 
presented Table 29. this table, Column (5) contains values de- 
termined the manner described the preceding paragraphs. Column (6) 
gives the efficiency the joint, the efficiency being the ultimate load, 
determined test, divided the product the gross area and the unit 
strength the control specimen. Column (7) gives the effective net section 
given test, this being the ultimate strength, given the test, divided 
the unit strength the control specimen. Columns (8) and (9) give the 
computed values the net section, the former being the specifi- 
cations for railway bridges and the latter the proposed specifications. 
Columns (10) and (11) give the ratios the computed the test values the 
net section, the former being for the specifications and the latter for 
the proposed specifications. value greater than unity indicates that the 
computed value exceeds the value determined test. perusal the values 
Columns (10) and (11) reveals the agreement between the computed values 
and the values determined from tests. 

interest note that the largest value given Column (10) 
1.149, Line 21, and the smallest value 0.820, Line 15, indicating that the 
specifications gave value for the net section 15% too great for one 
specimen and 18% too small for another specimen. The largest value 
given Column (11) 1.061, Line 21, and the smallest value 0.928, 
Line 15, indicating that the maximum difference between the two values, one 
given the proposed specifications and the other the tests, only per cent. 
noted further that, general, the proposed specifications give value 
more nearly equal the value given tests than the specifi- 
cations. Moreover, the proposed specifications penalize the bad practice 
using extremely large transverse spacing the rivets and large edge distances 
cations, result high efficiency but which, according tests, resulted 


little, any, increase efficiency. 
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the proposed rule, defined the transverse spacing, two times the 
edge distance, whichever the greater. the transverse spacing varies, the 
widest space, two times the edge distance, used. 


RIVETS AND CARBON-STEEL PLATES 


No. 


ciency Railway 
inches joint Engi- Pro- 
neering posed 
Associa- specifi- 
tion cations 
specifi- 
cations 


(4) (5) (6) (7) (8) (9) (10) (11) 


(a) WATERTOWN ARSENAL TESTS 


eonnoon 
oounoou 
00 00 GO 
OR 
WW 
Ore © 
N 


Navy 


14.04 
13.79 

13.29 

13.13 

13.01 

13.72 

13.86 


NAAN 
Grower cr or orc 
ooooooo 


(c) Tests Harris 


3.00 0.812 12.20 10.00 11.32 
8.00 0.813 9.64 10.00 9.55 


1.010 
0.994 

0.973 

0.931 

1.061 

1.007 

1.023 

1.035 
0.947 
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all the transverse spaces are equal and the edge distance one-half 
the transverse spacing, the net section can expressed fraction the 
gross area and the value the fraction depends only upon the ratio 


The relation between the ratio Effective net section 


Gross section 


Fig. 24(b). apparent from this diagram that the efficiency joint does 


Thick- Diam- 
plate rivet 
| | | 4 
0.840 
0.823 
0.811 
0.809 
0.866 
0.851 
0.935 
0.910 
1.149 
1.032 
1.128 
1.119 
1.044 
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not increase appreciably with increase the transverse spacing beyond 

The authors have drawn the conclusion that the efficiency riveted 
tension member does not greatly exceed 75%, irrespective the transverse 
spacing the rivets. different rivet pattern had been used for the joints 
with large transverse spacing, believed that different conclusion would 
have resulted. All the specimens with large transverse spacing that were 
tested the authors had rivets the ends the outside row, and the edge 
distance was abnormally large for the weakest section the main plate. 
Thus, for Specimens CCC2 and DCC2, shown Fig. 22, the edge distance for 
the outside row rivets 6.5 in. Geometrically, this equivalent 
transverse rivet spacing in., whereas, for the interior rivets the outside 
row, the spacing only 8in. Moreover, the authors report another part 
the paper that, beyond the yield point, the unit longitudinal strain wide 
plate greater over the middle than the edges. This being true, the 
harmful effect placing the rivets the weak section near the mid-width 
accentuates the harmful effects wide edge distance. This statement 
supported the tests reported Table 29, which are discussed the following 
paragraphs. 

The only specimens for which the value the effective net section de- 
termined the proposed specification exceeds the value given the tests 
for those tests reported Lines 17, 21, 22, 23, and 24, Table 29. 

The details the specimen used the test reported Line are shown 
Table (Specimen 8510). The edge distance in. The proposed 
specifications give value for the effective net section 0.70% greater than the 
value given the test. 

The details the specimen used the test reported Line are also 
given Table (Specimen 9375). The edge distance 5in. The proposed 
specifications give value for the effective net section 1.2% greater than the 
value given the tests. 

The specimens used the tests reported Line 17, Table 29, had close 
spacing and small edge distance; but the fact that the computed value 
the effective net section exceeded 1.0% the value given the tests 
accounted for the fact that the failure the specimens was due the 
shearing the rivets and not failure the plate. All other specimens 
for which the computed value the effective net section exceeded the test 
value had the same rivet pattern the specimens Fig. 22. That is, the 
few instances which the computed value the effective net section exceeded 
the test value can accounted for the fact that the edge distance was 
large for the outer row rivets. 

The specimens shown Table have practically the same transverse 
spacing rivets, but Specimens DCC2 have two rivets near the middle the 
outer row, whereas Specimens Type have one rivet the center and one 
each end the outer The computed and test values the efficiency 
are equal each other and equal 0.80 for Specimens DCC2, whereas, for 
Specimens Type the computed efficiency 0.80 and the value determined 
from the test 0.85, the values being the averages from three tests each 
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instance. That is, although the computed value the effective net section 
for Specimens DCC2 reduced because the large edge distance, the computed 
value still greater, relative the test value for these specimens, than 
for the specimens Type which have small edge distance. 

The fact should noted, however, that the specimens used for the tests 
reported Lines and Table 31, had edge distances equal one-half 
the transverse spacing and yet they developed effective net width somewhat 
greater than the computed value. 


TABLE RIVETED TENSION MEMBERS FOR WHICH THE 
DISTANCE THE RIVETS THE Row EXCESSIVE 


2" Stee! 


Steel Rivets 


N 
“20. 27" x 
20.25" + 
SPECIMEN 8510 SPECIMEN 9375 
Joint Failed in Plate, Holes Drilled Holes Drilled 
Description 
Stress gross pounds per square inch 960 900 
Strength control specimen, pounds per square inch............ 550 700 
4 Effective net section, by test, in square inches .................+.- 8.32 8.41 
Net section American Railway Engineering Association specifi- 


result the analysis the data contained Table 29, the following 
rules for the design riveted tension members are submitted for the considera- 
tion the profession: 


Rule edge distance shall not exceed 2.5 for the outer row 


the foregoing equation (Rule 1): gross width the member; 
rivets plus in.; and transverse distance between rivets two times the 
edge distance, whichever the greater. 
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For joints with rivets omitted from the outer row, only rivets the outer 
row are considered determining and provided the pitch equals 
exceeds 0.80 the gage; otherwise all rivets the outer row and the inter- 
vening rivets the second row should considered determining and 


Upon STRENGTH RIVETED TENSION MEMBERS 


1" Rivets 


223" x2" 


4 @ 4" =16" 
Steel Rivets 
SPECIMEN DCC2 TYPE 


Carbon-Steel Plates Holes Holes Drilled 
Sub-Punched and Reamed 


SPECIMEN, 
Specimen DCC2 


Description 
Me- High- 

dium test 

steel steel 


Item No. 


inches inches 


691.0 686.0 669.0 902.0 113.0 
52.7 51.6 55.9 66.2 
65.6 65.1 65.7 81.6 


11.62} 


the preceding discussion has been assumed that the weakest section 
riveted tension member through the outside row rivets the gusset plate 
the ends the member. The writer knows published reports 
tests riveted tension members with stitch rivets, but the tests plates with 
holes Mr. Harris, reported Fig. 21, indicate open hole does not 
weaken plate much rivet transmitting stress. not reasonable 
suppose that stitch rivet weakens plate more than open 
would seem reasonable, therefore, that the foregoing rules design, which 
are applicable where the rivets transmit stress, give conservative results where 
the rivets not transmit stress but only serve stitch rivets. The 
probable exception this statement where stitch rivets are used connect 
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Stress gross section, kips per square inch.......... 
Strength control specimen, kips per square inch .... rati 
Effective net section, test, square inches 
Net section American Railway Engineering Association app 
specifications, square 10.97 
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two parallel pieces steel having widely differing stress-strain diagrams. 
Under this latter condition one piece, because its lower ductility, might 
before the other had been stressed its ultimate strength. This con- 
dition not likely occur member fabricated carbon steel struc- 
tural grade because the part the stress-strain diagram near the maximum 
stress nearly horizontal for considerable length. However, the behavior, 
under similar conditions, member fabricated from parallel pieces made 
steel which not uniform its physical properties somewhat uncertain. 


profession indebted the engineers the San Francisco-Oakland Bay 
Bridge who provided for this important series tests, and those who were 
responsible for making the tests and analyzing the results. has been advo- 
cated some consulting engineers, series well-planned tests, made 
connection with the important structures which are built, would soon yield 
valuable information concerning the behavior structures. 

The fact that under different stresses ‘‘the plates riveted joint tend 
elongate unequal important and should used the analysis 
many structural problems well that the riveted joint. allow- 
able unit stresses increase, structural design must placed more scientific 
basis. The structural designer must aware the many indeterminate 
factors that affect the action structure and must have some knowledge 
the influence those factors. 

order determine the effective net section and the effect end pitch, 
additional data are needed concerning the behavior plates with open holes 
under stress. The lines fracture for the specimens shown Figs. and 
all appear originate the regions high stress concentration. When 
condition plastic yielding reached, redistribution stress takes place, 
but initial fracture may occur the region where yielding first occurs. stress 
concentration the edges hole has any bearing the initial fracture, 
quite probable that the present ideas concerning net section need 
revised. The writer the opinion that, additional study and tests, 
rational analysis may result. the meantime, the authors’ 
appear represent reasonable procedure. 

view the uncertain action the riveted joint Stage II, the authors’ 
recommendations for assuming equal shear per rivet seem fully justified and 
additional refinement likely result from further tests. 


Associate Prof., Civ. Eng., Columbia Univ., New York, 
Received the Secretary September 18, 1939. 
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DISCUSSIONS 


LARGE CORE DRILLS AID CONSTRUCTION 
CHICKAMAUGA DAM 


Discussion 


Esa. (by engineers not understand that 
drills which will bore holes large enough for man enter can used any 
ground, from the loose water-bearing ground containing some boulders de- 
scribed the paper, the hardest rocks known the mining industry. 
the work described, was necessary use casing. Diameter was lost every 
time was placed, that hole which started diameter might 


the hole went through heavy ground and gouge which had sup- 
ported without losing size. This was accomplished satisfactorily using 
expanding steel form, behind which quick setting concrete was packed. These 
large drills are capable higher development and have much wider field 
possible usefulness than generally recognized. 


paper James Lewis, Jr., Assoc. Am. Soc. E., was published June, 1939, 
Proceedings. discussion printed Proceedings order that the views expressed may brought 
before all members for further discussion the paper. 


With Idaho Maryland Mines Corp., San Francisco, Calif. 
Received the Secretary August 14, 1939. 


Sinking with Shot Drill, Idaho Maryland Mine, Grass Valley, Newsom and 
Jackson, Bureau Mines Information Circular 6923. 
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DISCUSSIONS 


STRESS DISTRIBUTION AROUND TUNNEL 


Discussion 


RICHMOND, AND JACOB FELD 


Esq. (by stress distribution due 
gravity semi-infinite elastic body pierced circular hole interesting 
problem and the author’s treatment complete one. The application 
the results practise depend, course, upon the accuracy with which the 
assumptions are fulfilled. the shape the tunnel other than circular (as 
often is), mathematical solution practically impossible and recourse must 
had experimental methods. 

interesting experimental approach this problem was developed 
Mr. whose method consisted casting flat plate agar-agar 
horizontal position. network lines was placed the plate, and hole 
was cut, shape appropriately representing the tunnel profile. When the 
model was placed vertical position, the strain the material due its own 
weight could measured from the distortion the originally parallel lines, 
and hence the stresses could determined. For round tunnel profile Mr. 
Yamaguti obtained good agreement between his mathematical solution and the 
results his model test. 

Mr. Yamaguti’s method test reproduces the author’s Case with the 
appropriate value Poisson’s ratio for agar-agar. The author’s Case could 
also reproduced this method. would only necessary cast the 
agar-agar plate vertical position, that the stress distribution the un- 
pierced plate would approximate hydrostatic pressure. 

second experimental approach the problem the method photo- 
elasticity. the tunnel profile some distance below the surface, the stress 
distribution may obtained approximately the method superposition, 
stated the author. Assuming that the original stress distribution 


paper Raymond Mindlin, Assoc. Am. Soc. E., was published April, 1939, 
Proceedings. This discussion printed Proceedings order that the views expressed may brought 
before all members for further discussion the paper. 

Asst. Prof. Mech. Eng., Univ. British Columbia, Vancouver, Canada. 

Received the Secretary June 15, 1939. 


the Investigation the Stress Distribution Tunnel with Agar-agar Model Experiments,” 
Yamaguti, Proceedings, Third International Cong. for Applied Mechanics, Vol. II, 1930, pp. 


Je 
oll 
¢ 
. 
a 
‘ 
‘ 


1466 RICHMOND STRESS DISTRIBUTION Discussions 


represented Case II, for instance, the boundary stress around the profile 
just determined for uniform compression vertical direction magnitude. 


CO, = 


which the symbols used are those the author Fig. and the distance 
the hole from the free surface. The stress the boundary then obtained 
adding the two boundary stress distributions. For non-circular profile 
the boundary stresses are most easily determined photo-elastic tests. 

This procedure sufficiently accurate the profile some distance below 
the surface. For the case shown Fig. 4(b) (that is, for value Poisson’s 
ratio and 2.0), this approximate method gives uniform compressive 
stress around the boundary the tunnel 3.76 which density 
the material, and diameter the circular The maximum 
stress the author’s computation 4.26w The error the approxi- 
mate method thus slightly more than 10%, which perhaps accurate enough 
for practical purposes. Thus, the conclusion may drawn that the tunnel 
more than four radii, two diameters from the surface, the approximate 
solution obtained superposition sufficiently accurate. 

few years ago the writer (working under the direction Nadai the 
Westinghouse Research Laboratories) made some photo-elastic tests holes 
tunnel profile the suggestion Carstarphen, Am. Soc. The 
profiles tested were squares with various fillet radii the corners; the ratio 
boundary stress distribution was obtained photo-elastically when the bakelite 
plate which contained the holes was subjected pure tension parallel 
axis symmetry perpendicular the sides the holes. The stress concen- 
tration factors (that is, the ratio the maximum stress the boundary the 
hole the stress the plate long way from the hole) from these tests are 


TABLE 


THE SIDE THE 


Description 


Stress concentration factor for tension compression 
Maximum compressive stress boundary tunnel profile 


shown Table For the sake completeness, the case the circular hole 


tensile stress field, for included, the results for which can 


cor 
Then the boundary stress determined, for stress right angles, magnitude 
in’ 
rat 
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th: 
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obtained computation.” interesting note that minimum stress 
factor seems result for value about Thestress con- 
centration factor for square hole with sharp corners would theoretically 
infinite. 

The stress distribution around the boundary tunnels similar cross- 
section was then computed from the boundary stresses determined the 
photo-elastic tests. The case solved the writer was the author’s Case 
which the principal stresses before the tunnel existed were given Equations 
(69) and results obtained this procedure for value Poisson’s 


ratio are shown Item Table These results should sufficiently 


accurate, provided that the tunnel least two diameters below the surface, 
that the material which the tunnel driven behaves elastically, and that the 
original principal stresses are assumed. 


this paper—namely, circular tunnel uniform elastic material—Mr. 
Mindlin’s solution the most complete available. For use practical cases, 
however, the following items must taken into account: 


internal stress the rock its natural condition, which stress may 
relieved the local vicinity the tunnel the removal the tunnel rock. 

percentage error resulting from the use the method and for- 
mulas this paper when rock conditions cannot assumed 
seldom that tunnel driven through uniform granite, which probably the 
nearest approach the material assumed this paper. The more usual cases 
are tunnels driven through laminar rock, such shales, schists, limestones. 
extremely doubtful whether the stress determinations resulting from the 
formulas this paper can applied case where blocky limestone 
encountered, or, take different example, where schist thin layers, lying 
edgewise the face the tunnel, encountered. 

effect stresses around the tunnel temporary bracing 
wedged against the exposed rock during driving, the effect the rock stresses 
the tunnel lined with concrete and grouted under pressure; also the effect 
various grouting pressures balancing over-balancing rock stresses. 


taking the foregoing suggestions into consideration, might pertinent 
consider the assumption made under the heading the Problem”’ 
that, the circular boundary, the normal stress and the shearing stress are 
equal zero. Such assumption not correct when rock encountered 
which requires bracing prevent movement. seems the writer that 
the three cases mentioned covering the initial states the rock, although 
they cover great range possibilities, such range chiefly the variation 
the Poisson ratio Case II. Case that hydrostatic pressure which 


Cons. Engr., New York, 
Received the Secretary August 1939. 
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may applicable tremendous distances below the earth’s surface rocks 
that are under such pressure that they are practically liquid, and Case 
would only apply where tunnel driven through blocks rock with open ver- 
tical seams, which case could not considered isotropic. With reference 
Case and its dependence upon the Poisson ratio the material, might well 
emphasized that such ratio may vary for rocks, depending upon the amount 
the stress. However, stated Mr. Mindlin’s paper, such variation has 
little effect the stresses themselves, except when the tunnel unusually close 
the open surface. for such condition (namely, shallow tunnels) that 
careful stress analysis most necessary. 


bast 
sect 
the 
Flo 
con 
thr 
des 

elin 

thr 
tior 

act 

Soe. 

p.9 


AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


DISCUSSIONS 


DESIGN HIGH-HEAD SIPHON SPILLWAY 


Discussion 


spillway design presented this paper may prove meritorious. Equation (4), 
based the assumption that free vortex flow prevails the bend the throat 
section siphon spillway simplifying step, and results doing away with 
the necessity selecting masked under the name 
Flow Coefficient. Therefore, his attempt rationalize the common procedure 
design for throat sections siphon spillways worthy commendation. 

Flow Assumption.—Certain questions and doubts arise, however, 
concerning the soundness the author’s assumption vortex flow siphon 
throat sections, and concerning the practical application his formula the 
design siphon spillways. attempt answer these questions and 
eliminate these doubts, the velocity distributions, they actually exist the 
throat sections three different siphon were compared, respectively, 


the velocity distributions computed from Mr. Rock’s formula (see Equa- 
tion (10)): 


loge 
The numerical calculations involved this comparison are given Table 


and the results are portrayed graphically Fig. Table first lists the 
actual operating conditions the three selected siphons. These data were 


paper Elmer Rock, Jun. Am. Soc. E., was published April, 1939, Proceedings. 
this paper has appeared Proceedings, follows: June, 1939, Stevens, Am. 


Prin. Engr., Corps Engrs., A., Bonneville, Ore. 
Engr. Office, Bonneville, Ore. 

Assoc. Engr., Corps Engrs., A., Bonneville, Ore. 
Received the Secretary July 25, 1939. 
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TABLE 3.—Comparison BETWEEN VELOCITIES COMPUTED THE 
(Equation (11)) anp VELOCITIES 


Description Average 


Discharge cubic feet per second 
Operating head feet 27. 


Dimensions Throat Section, Feet: 


Radius the: 
Mid-section, ray 
Summit, 


our NO 
—+ 


0.833 0.605 0.466 
Discharge per unit width cubic feet per 
second 78.0 55.3 44.4 
Computed Velocities the Throat Section, feet per 
second: 


Mid-section, 
Summit, 


(c) (SEE (12)) 


Hydraulic Heads, Feet: 
Elevation points above the crest, the: 
The ratio p/w,|| the: 
Velocity the: 
Actual Throat Velocities, Feet Per Second: 


OR 
oo 


Ons] 


Percentage Differences Equation (13): 
Crest [(Item —Item 24)/Item 24] 100 
Mid-section [(Item —Item 25)/Item 25] 100... 
Summit [(Item —Item 26)/Item 26] 100...... 


See Transactions, Am. Soc. E., Vol. (1934), Fig. 989. 
Assumed average values order start the radii common 
Transactions, Am. Soc. E., Vol. (1934), Column (8), Table 995. 
Assume the datum Elevation 734.0 equal the elevation the crest all three spillways. 


Transactions, Am. E., Vol. (1934), Columns (4), (7), and (10), Table 994. 
Velocity head computed the difference between the energy gradient and the hydraulic gradient. 
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used Mr. Rock’s formula compute the velocities the crest, mid-section, 
and summit each siphon throat section. second set similar velocities 


——— Values Computed by Equation (11) 
__ Values Computed by Test Data (Total Head Minus 


Radius, 


Fic. 


these points was then computed subtracting the actually measured hydraulic 
gradient from the energy gradient. That is, 


(12) 
Comparison was finally made between the two sets velocities percentage 


difference order obtain index the validity Mr. Rock’s formula. 


V Rock V actual 


Percentage difference 


V actual 


will noted from Fig. that all three cases, the actual velocity values, 
derived from the test data, are greater the crest than the summit, and that 
the velocity values computed the author’s formula (Equation (11)) 
conform fairly well with the so-called actual values. When the number 
ways which experimental and personal error may have entered into the test 
data are considered, the comparison Fig. believed surprisingly good; 
and, furthermore, this comparison appears strong indication that the 
author’s assumption vortex flow justified. addition the foregoing 
comparison, the practical application Mr. Rock’s formula was tested 
checking the throat sections several other existing siphon spillways. was 
found that the computed and actual values checked quite satisfactorily. 
Further author’s treatment the importance 
proper lower leg design high-head siphons considered good. The mathe- 
matics involved this phase, well those concerning the derivation his 
vortex flow formula, were found valid. The numerical solution the 
example siphon, given Mr. Rock, aptly illustrates his point that plot 
the hydraulic gradient does not result true negative pressures the crest and 
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summit the throat section. However, the author had included his 
the actual prototype measurements and model experiments, which states 
prove conclusively that the negative pressures the throat section obtained 
from plotting the hydraulic gradient are erroneous, the value this part 
his paper would increased materially. 

Conclusions.—The inherent value Mr. Rock’s method that has 
extended the current method throat section design (which assumes 
velocity distribution) additional brief, his method consists 
using the distribution velocities the throat basis for determining the 
maximum velocities that will obtained under given discharge. The 
limitation these maximum velocities correspond maximum allowable 
negative pressures then permits more logical approach for the design the 
section. possible that when such siphon, throat-velocity, 
distribution theory completely established, the next logical step would lead 
attempt increase the efficiency siphon operation redesigning the 
throat section into other than rectangular circular section. Such design, 
course, would have result gain efficiency operation sufficient 
overcome increased construction costs. 
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